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I.  INTRODUCTION 

A  proper  understanding  of  mixing  in  chemical  reactors  is  of  vital 
importance  in  many  industrial  applications.  Liquid  phase  mixing  plays 
a  critical  role  in  precipitation  processes  where  the  shape  and  size  of 
relatively  insoluble  precipitate  can  be  significantly  altered  bv 
changing  the  local  mixing  environment  near  the  feed  port  (Garside, 

1985).  Similarly,  the  degree  of  mixing  in  the  reactor  strongly 
influences  the  initiator  consumption  in  some  commercial  polymerization 
reactors  (Van  der  Mblen  et  al. ,  1982).  Mixing  is  also  recognized  as  a  kev 
factor  affecting  overall  performance  of  a  combustor — both  in  terms  of 
fuel  efficiency  and  formation  of  pollutants  (Pratt,  1982). 

The  classical  perfectly  mixed  reactor  model  is  one  of  the  most 
useful  models  in  chemical  reactor  design.  However,  the  model  is  based 
on  the  assumption  that  immediately  after  entering  the  reactor,  the  feed- 
stream  is  distributed  uniformly  at  the  molecular  level  throughout  the 
reactor  volume,  and  there  are  no  spatial  gradients  of  concentration 
or  temperature.  Furthermore,  whether  the  reactants  are  fed  premixed 
or  separately  is  of  no  consequence.  Thus,  the  model  completely 
ignores  the  details  of  the  mixing  process  which  precedes  the  chemical 
reaction.  For  fast  reactions,  where  the  mixing  rate  is  comparable  or 
slower  than  the  reaction  rate,  the  perfectly  mixed  reactor  model  is 
clearly  inadequate.  The  problem  of  modeling  mixing  and  fast  reactions 
is  an  old  one  and  has  spurred  publication  of  several  plausible 
models  in  the  past.  However,  it  is  difficult  for  a  practicing 
engineer  to  choose  among  these  diverse  models  with  any  degree  of 
certainty  as  very  few  have  been  validated  against  a  well-defined  set 
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of  experimental  data.  The  principal  aim  of  this  work  is  to  obtain  data 
to  discriminate  among  several  rival  models  by  experimentally  studying 
the  effect  of  mixing  on  the  selectivity  of  two  fast  competing 
reactions. 

The  thesis  is  organized  as  follows.  In  Chapter  II,  several 
popular  approaches  to  chemical  reactor  modeling  are  reviewed.  A  par¬ 
ticular  class  of  models  called  "mechanistic  models"  is  given  special 
emphasis.  Five  mechanistic  models  having  a  turbulence  analogy  with 
known  statistical  turbulence  theories  are  briefly  described. 

Chapter  III  is  devoted  to  the  theoretical  foundations  of  the 
present  work.  The  turbulence  analogy  of  five  mechanistic  models  is 
evaluated  by  using  published  single  second-order  reaction  data.  This  is 
followed  by  numerical  simulation  of  two  competing  reactions  in  a 
highly  segregated  plug  flow  reactor.  The  results  of  these  simulations 
provide  a  background  for  the  experimental  design. 

The  details  of  the  experimental  design  are  considered  in  Chapter 
IV  while  Chapter  V  is  devoted  to  measurement  techniques,  methods  of 
analysis  and  details  of  the  operating  procedures. 

Experimental  results  are  presented  in  Chapter  VI  along  with  a 
comparison  of  model  predictions  with  the  data  and  a  discussion  of  the 
results.  The  major  conclusions  of  the  work  are  highlighted  in  Chapter 
VII.  Finally,  possible  avenues  of  further  research  are  suggested  in 


Chapter  VIII. 


II.  LITERATURE  REVIEW 

A.  General  Review 

A  great  deal  of  work  has  been  done  on  mixing  from  a  chemical 
engineering  viewpoint,  starting  in  1953  with  Danckwerts'  pioneering 
papers.  In  his  first  paper,  Danckwerts  (1953a)  introduced  the  concept 
of  residence  time  distribution  (RTD) .  RTD  is  a  useful  characteriza¬ 
tion  of  large-scale  fluid  motion  which  has  come  to  be  known  as 
macromixing.  In  another  paper  (Danckwerts,  1953b),  he  considered  the 
characterization  of  "local  texture,"  subsequently  known  as  micro¬ 
mixing,  of  a  nonuniform  mixture.  He  defined  two  characterizing  quan¬ 
tities,  namely,  the  scale  and  intensity  of  segregation.  The  scale  of 
segregation  is  related  to  underlying  fluid  mechanics  and  is  analogous 
to  the  scale  of  turbulence  used  in  the  statistical  theory  of  turbulence 
On  the  other  hand,  the  intensity  of  segregation  depends  upon  the 
interaction  of  diffusion  fluxes  with  the  chemical  kinetics  in  the 
context  of  chemical  reactors.  Although  Danckwerts  showed  how  to  use 
a  diffusion  controlled  reaction  (e.g.,  acid-base  neutralization)  to 
determine  the  intensity  of  segregation  (Danckwerts,  1958),  his  method 
of  characterizing  micromixing  did  not  find  much  practical  use  because 
of  inherent  experimental  and  conceptual  limitations.  For  example,  it 
was  not  possible  to  predict  intensity  of  segregation  for  a  slower 
reaction  from  the  acid-base  measurements.  Yet,  the  concepts  of  scale 
and  intensity  of  segregation  have  provided  the  framework  for  nearly 
all  the  subsequent  work  in  this  area.  Most  of  the  mechanistic 
modeling  has  proceeded  along  two  distinct  directions.  In  one,  much 
emphasis  is  laid  on  the  size  of  the  segregation  scale  which,  in  turn, 
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influences  the  intensity  of  segregation  (e.g.,  slab-diffusion  tvpe 
models  to  be  described  later  on).  In  the  other,  no  attention  whatso¬ 
ever  is  paid  to  the  scale  of  segregation  and  yet,  based  on  a  proposed 
physical  picture  of  turbulent  mixing,  all  levels  of  segregation 
intensity  are  accounted  for  (e.g.,  environment  and  coalescence- 
redispersion  models  to  be  described  later  on) .  Bourne  (1983)  has 
chronicled  part  of  these  developments  from  a  historical  perspective. 

Substantial  progress  is  also  being  made  in  modeling  turbulent 
reacting  flows  based  directly  on  the  Navier-Stokes  and  component 
continuity  equations.  However,  the  flow  in  most  systems  of  practical 
interest  is  sufficiently  complex  that  accurate  modeling  is  a  formidable 
task  (Patterson,  1981).  Since  mechanistic  models  can  accommodate  many 
of  the  complexities  of  practically  relevant  chemical  reactors,  our 
discussion  is  limited  only  to  the  mechanistic  models. 

Villermaux  (1983)  provides  an  extensive  review  of  mechanistic 
models.  All  of  the  mechanistic  models  contain  at  least  one  "micro¬ 
mixing  parameter"  whose  magnitude  determines  the  level  of  turbulent 
micromixing  in  the  reactor.  The  micromixing  parameter  has  generally 
been  treated  as  an  empirical  parameter  available  for  correlating 
experimental  data.  Because  of  the  supposed  empirical  nature  of  the 
micromixing  parameter,  most  of  the  mechanistic  models  have  been  viewed 
as  lacking  a  physical  basis  and  have  thus  been  regarded  with  skepticism 
— particularly  in  the  context  of  such  important  applications  as 
reactor  scaleup.  The  overwhelming  variety  of  models  have  also  led 
some  to  believe  that  given  a  little  imagination,  anv  number  of  equally 
plausible  models  can  be  conjured  up  (Villermaux,  1983). 

In  a  recent  work,  Tarbell  and  Mehta  (1985)  have  shown  that  five 


mechanistic  models  (the  coalescence-redispersion  (CRD)  model  of  Curl 
(1963),  the  slab-diffusion  (SD)  model  of  Mao  and  Toor  (1970),  the 
interaction-by-exchange  with  the  mean  (IEM)  model  of  Villermaux  and 
Devillon  (1972),  the  three-environment  (3E)  model  of  Ritchie  and 
Tobgy  (1979),  and  the  four-environment  (4E)  model  of  Mehta  and  Tarbell 
(1983a)) are  analogous  to  theories  of  isotropic  turbulent  micromixing 
as  expressed  by  Corrsin  (1957,  1964)  and  Rosensweig  (1964,  1966),  The 
analogy  is  based  on  the  analytical  demonstration  that  the  decay  law 
for  the  variance  of  concentration  of  a  nonreactive  tracer  as  derived 
from  the  mechanistic  models  has  the  same  form  as  that  derived  from 
direct  turbulence  theory.  For  a  batch  (plug  flow)  mixer  the  variance 
decay  takes  the  form  (Tarbell  and  Mehta,  1985) 

2 

V:r  exp(-T/Tm)  ,  (1) 

a 

o 

where  I  is  the  intensity  of  segregation,  X  is  a  time  scale  for 
s  m 

turbulent  micromixing  which  is  characteristic  of  the  fine  structure  of 

the  turbulence,  and  x  is  the  time  (residence  time).  The  following 

estimates  of  x  have  been  provided: 
m 

2  1/3  1/2 

X  -  0.5[4(L  / z)  +  (v/e)  In  Nc  ]  ;  Nc  >>  1  (liquids)  (2) 

ms  be  be 

T  -  - (5/t)2/3(L  /e)1/3  ;  NL  <  1  (gases)  (3) 

(3-N“c) 

where  Lg  is  the  integral  length  scale  of  the  concentration  field,  c  is 
the  turbulent  kinetic  energy  dissipation  rate  per  unit  mass,  v  is  the 
fluid  kinematic  viscosity,  and  N^c  is  the  Schmidt  number.  All  of  these 
quantities  are  directly  measurable  or  readily  estimated. 
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The  relationships  of  the  micromixing  parameters  of  the  mechanistic 

models  to  the  time  scale  of  isotropic  turbulent  mixing  (t  )  as  required 

m 

by  the  turbulence  analogies  are  summarized  in  Table  1.  Figure  1 
summarizes  the  physical  picture  underlying  each  model  for  the  case  of 
a  plug  flow  mixer  fed  with  two  separate  streams  carrying  A  (dark)  and 
B  (light),  respectively-  The  left-hand  panel  in  Figure  1  represents 
the  initial  inhomogeneity  while  movement  to  the  right  shows  the  decay 
of  inhomogeneitv  at  successive  residence  times.  Each  model  is  briefly 
described  below. 

B.  Mechanistic  Models  with  a  Turbulence  Analogy 


The  CRD  model  was  originally  proposed  by  Curl  (1963)  to  describe 
dispersed  phase  droplet  mixing  and  concurrent  chemical  reaction.  In 
Curl's  view,  a  chemical  reactor  population  consists  of  a  large  number 
of  equal  size  droplets  each  having  uniform  concentration  and  behaving 
as  an  ideal  batch  reactor.  Mixing  takes  place  by  random  and  instan¬ 
taneous  pairwise  coalescence,  concentration  homogenization,  and 
redispersion  into  an  identical  pair  of  droplets.  In  the  turbulence 
analog  for  a  homogeneous  fluid,  we  might  view  Curl's  droplets  as 
(isotropic)  turbulent  eddies  and  the  coalescence  rate  as  a  measure  of 
the  rate  of  turbulent  micromixing.  CRD  is  basically  a  deterministic 
population  balance  model  and  the  governing  equation  is  an  integro- 
differential  equation  for  the  probability  density  function  (Curl,  1963; 
Evangelista  et  al.,  1969;  Kattan  and  Adler,  1972).  The  micromixing 
parameter  of  the  model,  I,  is  the  average  number  of  collisions  experi¬ 
enced  by  an  eddv  during  its  passage  through  the  mixer. 


Summary  of  Turbulence  Analogies 


Micromixing 

Parameter 


Analogy 


I  =  — 
T 

m 


6  =  /  27TDT 


R  =  - 


R  -  — 

s  I 
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Comments 


Valid  at  all  flow  rate  ratios 


Approximately  valid  at  equal 
flow  rates,  invalid  other- 


Valid  at  all  flow  rate  ratios 


Valid  at  all  flow  rate  ratios 


Approximately  valid  at  equal 
flow  rates,  more  accurate 
at  unequal  flow  ratios 


en tat  ion  of  the  micromixing  process  for  various  models:  (A)  CRD 
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1 .  Slab  Diffusion  (SD)  and  Related  Models 

The  SD  model  of  Mao  and  Toor  (1970)  was  developed  to  describe  a 
chemical  reactor  with  a  plug  flow  RTD  and  separate  feedstreams  of 
reactants  (A  and  B) .  In  the  Mao  and  Toor  picture,  the  reactor  fluid  is 
made  up  of  two  planar  slabs  with  a  common  interface  open  to  diffusion 
and  impermeable  outer  surfaces.  The  slabs  enter  the  reactor  completely 
segregated  (pure  A  in  one  and  pure  B  in  the  other)  and  they  proceed 
to  intermix  by  molecular  diffusion  and  undergo  chemical  reaction  as 
thev  move  through  the  reactor  at  the  mean  fluid  velocity.  The  mixing 
parameter  of  the  model  is  the  double  slab  thickness,  5.  In  its 
original  formulation  the  SD  model  applies  only  to  reactors  with  equal 
flow  rate  reactant  streams,  and  the  turbulence  analogy  (Table  1)  only 
applies  to  that  case. 

For  unequal  flow  rates,  if  flow  inequality  is  reflected  through 
unequal  slab  thicknesses,  the  analogy  does  not  hold.  On  the  other  hand, 
if  slab  size  is  assumed  to  be  solely  determined  by  '  al  fluid 
mechanics,  one  could  argue,  as  Baldvga  and  Bourne  (1984)  have  done, 
that  slab  thicknesses  should  be  independent  of  flow  rate  ratio  for 
dilute,  isothermal  systems.  However,  in  that  case  it  is  not  clear  how 
the  model  can  account  for  the  effect  of  feed  flow  rate  variation  on 
reactor  performance.  It  is  simplv  concluded  that  the  SD  model,  in  its 
original  formulation,  does  not  possess  a  turbulence  analogy  when  feed 
flow  rates  are  unequal. 

The  original  SD  model  is  limited  to  reactors  with  a  plug  flow  RTD 
although  buerden  and  Thoenes  (1972)  have  extended  the  model  to  include 
an  arbitrarv  RTD  (with  the  restriction  that  the  RTD  for  each  feedstream 
be  the  same).  In  the  past  few  years,  Ottino  et  al.  (1979)  have 
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developed  a  "lamellar"  mixing  model  which  provides  a  generalization 
of  the  SD  model  through  incorporation  of  local  fluid  deformations 
(shear,  extension,  etc.)  as  well  as  micromixing.  However,  these 
models  are  limited  by  the  fact  that  flow  elements  are  segregated  from 
their  surroundings,  and  thus  "age"  mixing  is  prohibited.  The 
"stretching  lamina”  model  of  Angst  et  al.  (1982)  is  essentially  a 
much  simpler  (vis-a-vis  Ottino  et  al.)  version  of  an  SD  model  with 
local  fluid  deformation.  However,  it  suffers  from  the  physically 
unrealistic  assumption  of  a  diffusion  barrier  to  one  of  the  reactants. 
The  model  of  Angst  et  al.  (1982)  does  not  possess  a  turbulence  analogy 
in  the  same  sense  as  referred  to  earlier. 

3 .  Interaction  by  Exchange  with  the  Mean  (IEM)  Model 

As  pointed  out  by  Villermaux  (1984),  the  basic  idea  of  the  IEM 
model  was  suggested  by  Burgers  as  early  as  1951  (Burgers,  1951). 
However,  the  model  was  independently  developed  in  its  present  form  by 
Villermaux  and  co-workers  (1972,  1975)  as  an  outgrowth  of  the  work  of 
Harada  et  al.  (1962)  and  was  originally  conceived  for  application  to 
stirred  reactors  only.  For  a  two-feedstream  perfectly  stirred  reactor, 
the  IEM  model  consists  of  two  eddies — one  associated  with  each  feed- 
stream — which  act  as  well-mixed  batch  reactors  having  mass  exchange 
with  a  mean  environment  whose  concentration  is  assumed  to  be  constant 
and  to  coincide  with  the  mean  concentration  leaving  the  reactor. 
Micromixing  is  described  by  a  first  order  exchange  process  between  the 
mean  environment  and  each  eddy  (the  eddies  do  not  interact  directly) 
with  concentration  difference  as  the  driving  force.  The  exchange  rate 
is  characterized  by  a  mass  transfer  coefficient,  h,  which  is  the 
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micromixing  parameter  of  the  model.  The  turbulence  analogy  for  the 
IEM  model  was  apparently  first  recognized  by  Costa  and  Trevissoi 
(1972)  in  the  context  of  perfectly  stirred  reactors.  In  this  work  the 
IEM  model  is  applied  to  plug  flow  reactors.  It  is  easy  to  show  that 
the  same  turbulence  analogy  applies  to  plug  flow  reactors  by  inserting 
a  plug  flow  residence  time  distribution  function  into  the  analysis  of 
Tarbell  and  Mehta  (1985;  Equations  (18)-(23)). 

4 .  The  Three-Environment  (3E)  Model 

The  3E  model  was  developed  by  Ritchie  and  Tobgy  (1979)  to 
describe  chemical  reactors  with  two  feedstreams.  It  is  an  extension 
of  "two-environment  model"  of  Ng  and  Rippin  (1965)  for  premixed  feed 
reactors  to  separate  feedstream  reactors.  In  the  3E  model,  two 
entering  environments  (E.E.'s),  one  for  each  reactant  feedstream  in  a 
two-reactant  system,  are  assumed  to  supply  a  single  leaving  environment 
(L.E.)  at  rates  proportional  to  their  respective  masses.  While  each  E.E. 
acts  as  a  totally  segregated  reactor  (no  age  or  species  mixing),  the 
L.E.  behaves  as  a  maximum  mixedness  reactor  (complete  age  and  species 
mixing).  As  pointed  out  by  Mehta  and  Tarbell  (1983a),  the  E.E.'s 
should  be  viewed  as  lumped  representations  of  turbulent  eddies  of  pure 
reactant  which  have  only  interacted  with  other  pure  reactant  eddies 
and  as  such,  they  model  the  extremes  of  the  turbulent  concentration 
spectrum.  The  L.E.  is  a  lumped  representation  of  eddies  in  the  inter¬ 
mediate  region  of  the  turbulent  concentration  spectrum.  Micromixing 
is  modeled  by  the  first  order  transfer  of  material  from  the  E.E.'s  to 
the  L.E.  with  a  transfer  coefficient  R  which  is  the  model's  micro¬ 


mixing  parameter. 
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The  3E  model  has  a  very  general  structure  which  can  accommodate 
an  arbitrary  RTD  for  each  feedstream  all  within  the  computationally 
efficient  format  of  ordinary  differential  equations.  Since  there  are 
only  three  effective  eddies  in  the  model,  the  computations  should  be 
modest  by  comparison  to  those  required  for  Monte  Carlo  simulations  with 
the  CRD  model  (300  eddies) . 

In  spite  of  many  positive  attributes,  the  3E  model  is  lacking  in 
an  important  structural  feature  which  the  CRD,  SD,  and  I EM  models 
possess — the  presence  of  reacting  regions  which  are  rich  in  each  of  the 
reactants  in  a  two-reactant  (A  and  B)  feedstream  reactor.  The  3E  model 
has  only  one  reacting  environment  with  reactant  concentrations  at  the 
mean  of  the  concentration  spectrum.  This  structural  feature  of  A  rich 
and  B  rich  regions  becomes  very  important  in  competing  reaction  systems 
where  the  stoichiometric  ratio  of  the  reactants  may  have  a  profound 
effect  on  the  selectivity.  This  point  is  discussed  further  in  a  later 
chapter . 

5.  The  Four-Environment  (4E)  Model 

To  overcome  the  limitation  of  only  a  single  reacting  environment 
in  the  3E  model,  Mehta  and  Tarbell  (1983a)  developed  the  4E  model  by 
introducing  a  separate  L.E.  for  each  feedstream  and  allowing  mutual 
interaction  between  the  L.E.'s.  The  reactant  stoichiometry  of  each 
L.E.  may  now  be  rich  in  the  reactant  fed  through  its  respective  E.E. 

The  4E  model  retains  the  basic  features  of  the  3E  model.  The 

E.E.'s  are  segregated  flow  reactors;  the  L.E.'s  are  maximum  mixedness 

reactors;  and  the  transfer  of  material  from  E.E.  to  L.E.  is  first  order 

in  the  mass  of  the  E.E.  with  a  transfer  coefficient  R  .  It  is  further 

s 
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assumed  that  the  reversible  transfer  between  L.E.'s  is  first  order  in 
the  mass  of  each  L.E.,  also  with  a  transfer  coefficient  R  .  The  two 
L.E.'s  of  the  4E  model  are  a  somewhat  finer  lumped  representation  of 
the  intermediate  region  of  the  turbulent  concentration  spectrum  than 
the  single  L.E.  of  the  3E  model.  But  more  importantly,  they  provide 
the  structural  feature  of  multiple  reacting  regions. 

The  mathematical  structure  of  the  4E  model  is  similar  to  that  of 
the  3E  model  with  an  additional  set  of  material  balances  required  for 
the  second  L.E.  (see  Mehta  and  Tarbell  (1983a)  for  complete  details). 
Thus,  the  computational  requirements  of  the  two  models  are  similar. 


III.  BACKGROUND  FOR  THE  PRESENT  WORK 


The  existence  of  a  turbulence  analogy  for  the  models  described 
earlier  blunts  the  criticism  that  all  such  models  are  useful  onlv  for 
data  correlation.  A  priori  estimation  of  the  micromixing  parameter 
is  possible  now.  Tarbell  and  Mehta  (1985)  and  more  recently,  Chang 
et  al.  (1985)  have  shown  that  all  five  mechanistic  models  which  have 
a  turbulence  analogy  are  well  supported  by  single  second-order 
reaccion  data  in  various  reactors,  and  that  the  micromixing  parameters 
required  to  fit  single  reaction  data  are  in  accord  with  Corrsin's 
expression  of  x  .  However,  these  models  are  not  all  equivalent. 

Chang  et  al.  (1985)  have  demonstrated  this  through  numerical  simula¬ 
tions  of  competing  reactions  in  a  PFR.  These  results  emphasize  the 
importance  of  certain  structural  features  of  the  mixing  models  which 
are  independent  of  the  existence  of  a  turbulence  analogy.  They  also 
point  towards  the  importance  of  obtaining  experimental  data  for 
competing  reactions  in  a  highly  segregated  PFR  to  discriminate  among 
rival  mixing  models.  The  remainder  of  this  chapter  is  based  on  the 
study  by  Chang  et  al.  (1985)  . 

A.  Experimental  Evaluation  of  the  Turbulence  Analogies 
for  Single  Reactions 

First,  the  simulation  of  representative  experiments  of  Vassilatos 

and  Toor  (1965)  are  considered  using  the  five  models  possessing  a 

turbulence  analogy.  In  these  often  quoted  experiments,  single, 

irreversible,  second-order  reactions  (A+B  -*  P;  r,  =  -kC,CD)  were 

A  A  d 

studied  in  an  isothermal,  turbulent  flow,  tubular  reactor.  The  hydro¬ 
dynamics  were  maintained  constant  while  the  reaction  rate  was  varied 
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over  nine  orders  of  magnitude.  For  brevity,  simulation  of  only  three 
experiments  are  considered,  although  they  span  the  entire  range  of 
reaction  rates  (see  Table  2).  In  these  experiments,  the  two  aqueous 
reactant  solutions  were  separately  introduced  through  100  alternate 
jets,  and  the  reaction  proceeded  in  the  resulting  microscale  inhomo¬ 
geneous  mixture.  Since  no  RTD  measurements  were  reported,  it  has 
generally  been  assumed  in  models  of  this  reactor  that  an  ideal  plug 
flow  RTD  is  appropriate.  Since  the  flow  rate  of  each  stream  was  main¬ 
tained  constant,  the  feed  stoichiometry  was  varied  by  changing  the 
concentration  in  one  of  the  feed  streams. 

In  a  previous  work,  Mehta  and  Tarbell  (1983b)  simulated  Vassilatos 

and  Toor's  data  with  both  the  3E  and  4E  models  assuming  a  plug  flow  RTD. 

For  each  model  they  determined  the  value  of  the  parameter  Rg  which 

provided  the  best  fit  (mean  square  error  minimization)  to  the  entire 

data  set.  That  resulted  in  a  value  of  t  0.0082  sec  through  use  of 

m 

the  turbulence  analogy  (Table  1) .  This  value  of  x  is  used  to 

m 

calculate  the  micromixing  parameter  for  each  model  based  on  the 
turbulence  analogy  (Table  1) ,  to  simulate  the  three  cases  listed  in 
Table  2.  The  results  are  displayed  in  Figure  2.  Numerical  methods 
are  discussed  in  the  next  section. 

For  the  "slow"  reactions  (Figure  2a) ,  the  data  and  all  of  the 
models  follow  the  maximum  mixedness  prediction  closely.  Both  the 
"moderately  fast"  and  "very  fast"  reactions  are  fairlv  well  described 
by  all  of  the  models  although  IEM  provides  the  best  and  3E  the  worst 
fit  to  the  data.  Overall,  the  comparisons  are  quite  favorable  for 
all  models,  and  this  lends  significant  support  to  the  turbulence 
analogies  since  each  model  employs  the  same  value  of  T  . 
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Experimental  evaluation  of  the  turbulence  analogy  is  pursued 

further  bv  noting  that  McKelvev  et  al.  (1975)  have  determined  x  in  a 

m 

duplicate  of  Vassilatos  and  Toor's  reactor  bv  measuring  the  decay  of 

concentration  fluctuations  with  distance  downstream  of  the  injector 

in  a  nonreactive  svstem.  i  (as  defined  bv  Equation  (1))  was  observed 

m 

to  be  a  function  of  axial  position:  its  value  was  0.0042  sec  at  1.1 

cm,  and  0.015  sec  at  3.02  cm.  No  significant  reaction  took  place 

beyond  this  entrance  region  of  the  reactor  for  the  "moderately  fast" 

and  "very  fast"  reactions  (Figure  2).  Thus,  the  average  value  of  x^ 

measured  in  the  reaction  zone  was  0.0096  sec,  which  is  within  17 

percent  of  the  value  determined  independently  by  fitting  the  mean 

axial  concentration  profiles  (x  =  0.0082  sec  was  employed  in  the 

m 

simulations  of  Figure  2).  This  lends  support  to  the  interpretation 
of  micromixing  parameters  in  terms  of  turbulence  properties. 

Treleaven  and  Tobgy  (1973)  studied  single,  second  order,  liquid 
phase  chemical  reactions  in  a  tubular  reactor  with  one  reactant 
entering  through  a  single  central  jet  and  the  second  reactant  through 
the  annulus  surrounding  the  jet.  Individual  feedstream  RTD  data  had 
been  obtained  previously  (Treleaven  and  Tobgy,  1972a) .  A  range  of 
feed  stoichiometric  ratios  (0.1  -  15.0)  and  flow  rate  ratios  (1.0  - 
3.0)  were  considered  with  reactions  which  could  be  described  as 
"moderately  fast"  in  the  sense  of  Vassilatos  and  Toor's  experiments. 

Treleaven  and  Tobgy  (1972b)  used  experimental  RTD  information 
and  a  CRD  model  to  simulate  the  experiments,  while  Ritchie  and  Tobgy 
(1979)  used  the  same  data  in  3E  model  simulations,  and  also  reported 
results  of  simulations  with  the  SD  model  (Ritchie  and  Tobgy,  1978). 

For  the  cases  which  were  simulated  bv  all  three  models,  the  best  fit 
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values  of  t  determined  bv  the  turbulence  analogies  (Table  1)  are 
m 

nearlv  identical — differing  by  at  most  8  percent.  These  observations 

provide  additional  support  for  the  turbulence  analogies. 

It  should  also  be  pointed  out  that  Mehta  and  Tarbell  (1983b) 

accurately  modeled  the  extensive  set  of  perfectly  stirred  reactor  data 

reported  by  Plasari  (1976)  with  the  3E  and  4E  models  and  found  values 

of  t  which  were  within  10  percent  of  estimates  based  on  Equation  (2). 
m 

The  comparisons  in  this  section  indicate  that  the  five  mechan¬ 
istic  models  and  their  turbulence  analogies  are  quite  strongly  supported 
by  experimental  data  involving  a  single  chemical  reaction.  This  raises 
an  important  question:  Is  it  possible  to  clearly  discriminate  among 
the  models? 

It  has  been  generally  recognized  that  rapid  competing  reactions 
of  the  type 

kl 

A  +  B  •+  R 

k2 

R  +  B  ->  S 

are  particularly  sensitive  to  mixing  as  reflected  in  the  selectivity 

(2Cc/[Cn  +  2Cc ] ) .  However,  Ritchie  and  Tobgy  (1974)  have  shown  via 
b  K  b 

numerical  simulations  that  the  CRD  and  IEM  models  yield  nearly  the 
same  selectivity  predictions  for  such  competing  reactions  in  a  perfectly 
stirred  reactor  when  the  micromixing  parameters  of  the  two  models  are 
related  by:  4x/I  =  1/h,  as  prescribed  by  the  turbulence  analogies 
(Table  1)  and  I  >  4  (we  will  return  to  this  point  in  a  later  section) . 
Prompted  by  this  and  other  model  comparisons,  Villermaux  (1983) 
suggested  that  all  such  mechanistic  models  are  approximately  equivalent 
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Indeed,  from  analysis  of  single  reaction  data,  one  would  tend  to 
concur.  However,  in  the  next  section  it  will  be  demonstrated  that  the 
models  are  not  equivalent  through  simulations  of  a  plug  flow  reactor 
with  competing  reactions.  Unmixed  feedstream  plug  flow  reactors  are 
intrinsically  more  sensitive  to  micromixing  than  stirred  tanks  beer  ’.se 
both  reactant  concentration  and  reactant  segregation  are  elevated  in 
the  entrance  region  of  the  reactor.  Backmixing  in  stirred  tanks  tends 
to  reduce  both  of  these  effects. 

B.  Numerical  Simulations  for  Competing  Reactions 

For  simulations,  a  tubular  reactor  with  two  separate  feedstreams, 

each  characterized  by  the  same  plug  flow  RTD,  by  individual  flow  rates 

and  Qfi,  and  by  unmixed  reactant  concentrations  and  is 

considered.  For  the  competing  reaction  scheme  the  kinetics  of  the 

diazo-coupling  reactions  reported  by  Bourne  and  co-workers  (1981)  is 

3  3 

employed.  The  kinetic  rate  constants  are:  k^  =  7.3  x  10  m  /mol/s; 
k^/k,,  =  2086.  Through  several  preliminary  simulations  it  was  found 
that  a  clear  discrimination  of  the  mixing  models  is  apparent  when  the 
micromixing  time  constant  is  set  at  =  0.315  s,  and  thus  all  of  the 
results  to  be  presented  rre  based  on  this  value.  This  is  some  38 
times  larger  than  that  d  5cribing  Vassilatos  and  Toor's  reactor,  and 
is  thus  characteristic  of  a  highly  segregated  reactor. 

1 .  Numerical  Techniques 

All  of  the  simulations  were  carried  out  on  an  IBM  3033  computer. 
The  CRD  model  was  solved  by  Monte  Carlo  methods  (Spielman  and 
Levenspiel,  1965;  Treleaven  and  Tobgy,  1972)  using  a  random  number 
generator  to  select  pairs  of  eddies  for  instantaneous  mixing  and 


redispersion.  Each  eddy  in  the  ensemble  was  treated  as  a  batch 
reactor  with  appropriate  material  balances  accounting  for  the  reac¬ 
tion  kinetics.  Between  collisions  the  eddy  material  balances  were 
integrated  by  simple  first  order  explicit  predictor-corrector  methods. 
Three  hundred  eddies  were  maintained  in  the  ensemble  and  no  signifi¬ 
cant  changes  in  predictions  were  observed  when  the  number  of  eddies 
was  increased  to  500.  Typical  computation  times  for  CRD  simulations 
varied  between  500  and  5000  CPU  seconds  while  all  of  the  remaining 
models  required  no  more  than  15  seconds.  The  3E,  4E,  and  IEM  models 
required  numerical  integration  of  systems  of  nonlinear  ordinary 
differential  equations  which  was  carried  out  with  standard  IMSL  routine, 
DV0GER  based  on  Gear's  method.  The  partial  differential  equations  of 
the  SD  model  were  solved  by  a  fully  implicit  finite  difference  method 
with  100  uniformly  spaced  grid  points  to  cover  the  double  slabs.  The 
3F  and  4E  model  equations  were  available  in  Mehta  and  Tarbell  (1983a). 
The  IEM  and  SD  equations  were  simple  extensions  of  equations  available 
in  Tarbell  and  Mehta  (1985),  incorporating  appropriate  reaction  rate 
terms  and  separate  material  balances  for  each  species. 

The  computer  codes  were  extensively  checked  for  accuracy  and 

convergence  (see  Appendix  A.1-A.5  for  computer  codes).  In  particular, 

the  predictions  of  all  five  models  were  compared  with  available 

analytical  expressions  for  two  asymptotic  cases:  (1)  diffusion  of  an 

inert  tracer  (k^  =  k9  =  0;  =  0),  (2)  competing  reactions  under 

maximum  mixedness  conditions  (t  -*■  0) .  In  the  first  case,  the  predicted 

m 

value  of  the  tracer  concentration  variance  was  compared  with  the  exact 
analytical  expression  presented  in  Tarbell  and  Mehta  (1985),  and  the 
agreement  was  always  found  within  0.1  percent  on  a  relative  basis.  For 


Che  maximum  mixed  condition,  verv  low  values  of  x  were  emploved  (x  = 

m  '  m 

0.0001  sec,  0.0003  sec)  and  the  predictions  were  observed  to  be 

independent  of  X  .  Next,  for  a  sequence  of  values  of  C,  (average 
m  A 

concentration)  predicted  bv  each  model,  the  associated  values  of  CD 

K 

and  CL  were  compared  with  values  of  CL  and  CL  calculated  from  analytical 

D  K  D 

expressions  available  in  a  standard  text  (Levenspiel,  1962): 


.  CA  . ,  CA  k2/kl 

-R  -  l-k2/kl 


and 


Once  again,  all  he  comparisons  were  within  0.1  percent  on  a  relative 
basis . 


2.  Simulation  Results  and  Discussion 

Five  different  cases  were  simulated  to  provide  a  basis  for 
comparison.  The  input  parameters  and  predicted  selectivities  (Xg) 
are  summarized  in  Table  3.  X<,  represents  the  fraction  of  B  consumed 

to  produce  S  when  3  is  completely  depleted.  Representative  concen¬ 
tration  profiles  of  products  R  and  S  along  the  reactor  length  are 
plotted  in  Figures  3  to  5.  Figure  3  is  the  base  case  in  which  flow 
rates  and  concentrations  of  each  reactant  feedstre3m  are  equal.  Figure 
9  shows  the  effects  of  a  large  variation  in  the  flow  rate  ratio 

(0,/0n  =  100),  and  Figure  5  shows  the  consequence  of  varying  the  feed 
'A  'D 

concentration  ratio  ( C . „ / C =  2). 

AU  BU 

The  S  concentration  profiles  (Figures  3-5)  and  selectivities  (X^, 
Table  3)  verv  clearlv  reveal  large  differences  in  the  predictions  of 


rvyyyj  ^  *•-  >  -  f  '  r  ^  ^  && 


CO 

CD 


XI 

G 

£-* 


uut. 

r-H 

in 

00 

• 

1 

1 

• 

SIh- 

CM 

1 

1 

CO 

wi 

r— 4 

r— | 

04 

»—  j 

1 

G 

X 

<?■  1  £2 

00 

<r 

04 

m 

C“N 

r— i 

<u 

1  hv; 

Cl  QJ 

CM 

d 

X) 

00 

oo 

TJ 

ail 

r-i 

*— i 

CM 

o 

O  1 

s 

~i 

0) 

> 

**H 

y~V 

r— 1 

n 

XI 

un 

<3\ 

r^» 

cn 

tz, 

VC 

1 

1 

o 

00 

IO 

C/3 

CM 

1 

1 

CM 

>, 

CM 

. 

. 

. 

+ 

O 

o 

o 

"U 

0) 

a: 

IO 

r-H 

00 

m 

o 

v'— ' 

CM 

m 

on 

CM 

•l—t 

\ 

w 

CM 

CM 

o 

CM 

i — < 

TJ  j 

CO 

1— 1 

• 

• 

• 

. 

, 

<u 

IO 

C 

o 

o 

o 

o 

CM 

o  1 

II 

w 

G 

CO 

vO 

x> 

vO 

CM 

X 

w 

CM 

CM 

CM 

CM 

i 

*-» 

on 

C 

c 

o 

O 

i 

O 

. 

• 

. 

XJ  | 

o 

C 

o 

o 

u 

G 

rj 

& 

r-H 

o 

<r 

sD 

O'- 

ON 

g 

W 

<r 

o- 

CO 

O 

04 

c 

<r 

T— < 

r—i 

»-H 

O 

r— t 

• 

• 

• 

• 

. 

ex  ! 

! 

o 

o 

o 

o 

o 

to 

3 

•— 1 

CJ 

CO 

<T 

CM 

C^j 

o 

! 

G 

<r 

0^ 

CM 

i—i 

* 

Q£ 

T— 1 

CM 

i- ^ 

U 

♦ 

• 

• 

• 

• 

c 

o 

o 

o 

O 

o 

»(“* 

w 

£ 

ca 

0 

O' 

r—t 

o 

o 

r— 

_ _ ( 

•r* 

\ 

r—H 

o 

xJ 

< 

r—i 

'J 

ri 

i) 

^  i 

cn 

04 

on 

OC 

1 _ 

1^ 

1 

c 

O 

O 

•r* 

o  ^ 

r*n 

c—4 

r-H 

, _ 1 

XJ 

K 

>4 

V 

X 

X 

O  t 

O  v~' 

sO 

r— ' 

on 

nj 

m 

d 

* 

oo 

o 

. 

. 

e 

o 

• 

on 

o  ! 

u  , 

u*  1 

o  ! 

O 

on 

1 

on 

i 

oo 

on 

1 

m 

1 

■> 

_ 

i 

»— • 

v—' 

u 

O 

>— i 

f— < 

■ — . 

— 1 

— 1 

— H 

<  ^ 

X 

\s 

X 

y 

> 

u  — ' 

vC 

no 

oo 

04 

vC 

•H 

. 

on 

o 

• 

# 

XJ 

•  • 

• — - 

• — • 

y  : 

3 

. 

* 

o 

u 

in 

un 

in 

in 

— - i 

r— 1 

r— < 

-n 

_ , 

w  cn 

on 

on 

0^> 

on 

' — ' 

• 

, 

• 

. 

. 

O 

'  , 

O 

O 

° 

0) 

! 

i 

« 

r— i 

CM 

on 

<T 

un 

l 

i 

i 

x 

>. 


X 

X 


G 

u 


c 

•*H 

X 


CC 


2 

-5 


‘O 

T3 

OJ 

no 

cn 

i 

Q 

C 

u- 

_  o* 

•"• 

O 

u 

o 

o 

XJ 

X3 

E 

• 

Cn 

G 

X 

03 

C 

e 

0 

'j 

"O 

XJ 

•r-< 

on 

r— i 

OJ 

• 

> 

v£> 

> 

•r-1 

XJ 

cn 

o 

■— ■ 

03 

£ 

Cfi 

3 

c 

o; 

1* 

Cn 

o 

-J 

c 

i> 

U1 

'J 

c~ 

o* 

> 

— * 

*• 

J 

>^ 

11 

> 

c£: 

• — 

o 

•r-< 

r> 

r— 

U— 

03 

3  1 

- < 

— 

G  = 

'•— ' 

d  G 

X 

- 1 

1) 

G 

03 

U 

Z  I* 

•<H 

G 

•■—4  _G 

0> 

XJ  U 

c3 

"H 

G 

G 

U  TO 

H 

XC 

i-i  CD 

ojD 

V- 

C  G 

o 

G 

1  •<“* 

£ 

XJ 

U  C3 

u 

c  e 

> 

0 

0  CD 

<3 

W 

03 

XI 

w 

vv 

>V 

VSjS 


ss| 

?i 


m 

>■«' »< 


wv 


bjj^w. 


•\n  ,-,; 

•,*  N  .V 


isr jrjvjv'jvvvvrjv'jv'j* w  vv  ^  v  vn.vw.m-*  i 


Figure  4. 


Axial  concentration  profiles  for  competing  reactions  in  a 
PFR  predicted  by  five  models — Case  3:  (A)  profile  of 

intermediate  products;  (B)  profile  of  final  product 
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Axial  concencration  profiles  for  competing  reactions  in  a 
PFR  predicted  by  five  models— Case  5:  (A)  profile  of 
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Che  models.  For  example,  there  is  nearlv  a  tenfold  difference  in 
ijredictud  b  production  between  i_he  CRD  ana  3E  moueis  tor  case  3,  and  a 
fourfold  difference  for  the  same  case  when  the  CRD  and  IEM  models  are 
compared.  On  the  other  hand,  the  CRD  and  4E  models  give  quanti¬ 
tatively  similar  predictions,  and  the  same  may  be  said  of  the  IEM  and 
SD  models  (see  the  last  two  columns  of  Table  3) .  Although  not  shown 
in  Figure  5,  it  should  be  noted  for  case  5  (excess  of  A)  that  both  the 
maximum  mixedness  and  3E  models  predict  no  production  of  S  over  the 
entire  reactor  length. 

Of  course,  it  is  not  possible  to  decide  from  hypothetical  simu¬ 
lations  which  model,  if  any,  is  correct.  This  most  important  question 
must  ultimately  be  decided  by  comparisons  with  well  designed  experi¬ 
ments.  However,  these  simulations  do  seem  to  indicate  that  the  3E 
model  is  inadequate  for  selectivity  estimation  as  its  predictions  are 
very  insensitive  to  the  parameter  variations  we  have  considered  and  are 
always  close  to  maximum  mixedness  values.  This  finding  is  consistent 
with  the  conclusions  of  Mehta  and  Tarbell  (1983b)  who  simulated  a 
methane  combustor  with  the  3E  and  4e  models.  The  failure  of  the  3E 
model  is  apparently  a  consequence  of  reactions  taking  place  in  a  single 
concentration  environment  (L.E.)  which  is  at  the  mean  concentration. 

All  of  the  other  models  contain  at  least  one  rich  reacting  region  for 
each  of  the  reactants.  This  all  points  to  the  importance  of  structural 
features  of  mixing  models. 

Because  of  the  observation  that  the  CRD  and  4E  models  give  similar 
predictions  which  are  significantly  different  from  those  of  the  IEM  and 
SD  models,  we  naturally  wonder  if  there  are  prominent  structural 
differences  between  these  pairs  of  models  which  may  be  responsible  for 
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the  observations.  As  mentioned  previously,  all  four  models  possess 
the  strut! rural  Feature  of  "rich  ''cuctinf,  regions"  for  each  reactant. 
However,  only  the  CRD  and  4E  models  contain  unmixed  (unreacted) 
regions,  whiTe  the  IEM  and  SD  models  do  not  (see  Figure  1).  The 
presence  or  absence  of  "unmixed  regions"  in  the  reaction  mixture  may  be 
important  in  selectivity  predictions,  particularlv  in  highly  segregated 
reactors.  In  fact,  in  simulations  of  the  parallel-consecutive  reaction 
scheme  in  a  CSTR  by  both  the  IEM  and  CRD  models,  Ritchie  and  Tobgy 
( 1 9 7 A )  found  that  in  highly  segregated  reactors  (I  <  4),  the  predic¬ 
tions  of  the  models  differed  substantially,  whereas  the  comparisons 
were  quite  favorable  in  less  segregated  reactors  (I  >  4) .  The  simulated 
plug  flow  reactor  above  is  characterized  by  a  maximum  value  of  I  ~  6 
at  the  end  of  the  reaction  zone,  and  this  value  of  course  decreases  to 
zero  at  the  reactor  entrance.  Thus  the  simulated  reactor  is  highly 
segregated  over  most  of  its  length. 

The  present  experimental  study  derives  its  motivation  from  these 
results.  The  details  of  the  experiments  are  considered  next. 


IV.  EXPERIMENTAL  DESIGN 


A.  General  Discussion 

The  goal  of  this  experimental  study  is  to  obtain  reliable  data 
from  a  highly  segregated  turbulent  flow  reactor  such  that  the  observed 
mixing  effects  are  significantly  larger  than  the  errors  in  the 
measurements.  Furthermore,  we  want  to  discriminate  among  rival  models 
by  comparison  of  model  predictions  with  the  experimental  data.  Thus, 
proper  choice  of  reactor  and  reaction  system  is  very  critical. 

Turbulence  is  a  complex  flow  phenomenon  and  its  interaction  wiLh 
chemistry  of  competing  reactions  makes  it  even  more  complex.  In  this 
context  it  is  desirable  to  use  as  simple  a  flow  field  as  possible  in 
our  reactor.  Conceptually,  the  simplest  turbulent  flow  is  a  completelv 
isotropic  one  in  which  all  turbulence  properties  are  invarian  to 
axis  rotation,  reflection,  and  translation.  Such  a  fluid  implies 
spatial  homogeneity,  but  it  would  be  decaying  with  time  because  its 
definition  excludes  existence  of  any  turbulence  generation  within  the 
field.  A  field  which  is  nearly  isotropic  can  be  produced  downstream 
from  a  periodic  grid  placed  perpendicular  to  the  mean  flow  in  a  pipe. 

A  similar  one-dimensional  isotropic  field  can  also  be  produced  down¬ 
stream  from  a  multitube  mixing  head  placed  in  a  pipe.  The  resultant 
turbulence  field  is  homogeneous  in  any  plane  parallel  to  the  grid  or 
the  mixing  head,  except  in  the  wall  boundary  layers,  but  decavs  in  the 
downstream  direction.  Thus,  the  turbulence  field  would  be  statisticallv 
stationary  in  space  but  inhomogeneous  in  the  streamwise  direction. 

This  is  essentially  a  plug  flow  with  fairly  well  defined  turbulence 
characteristics.  It  is  important  to  note  that  a  PFR  with  separate 
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feedstreams  is  intrinsically  more  sensitive  to  mixing  than  a  perfectly 
stirred  reactor  as  the  absence  of  any  backmixing  in  a  PFR  allows  much 
higher  reactant  concentration  and  segregation  in  the  entrance  region. 
Thus,  a  turbulent  PFR  appears  to  be  an  appropriate  choice  for  our 
purpose. 

Both  screens  (Keeler  et  al.,  1965)  and  multitube  mixing  heads 
(Vassilatos  and  Toor,  1965;  Mao  and  Toor,  1972;  Torrest  and  Ranz, 

1970;  Fisher,  1974)  have  been  employed  in  the  past  to  generate  a 
characteristic  turbulent  field  in  a  uniform  mean  flow  field.  Both  of 
them  provide  similar  flow  conditions  downstream  from  their  backfaces. 
The  flow  emerges  as  jets  which  coalesce  some  distance  downstream  of 
the  turbulence  generator.  The  region  up  to  the  coalescence  plane  is 
somewhat  ill-defined.  Previous  studies  of  very  fast  reactions  in 
multitube  mixing  head  reactors  have  shown  that  the  conversion  of 
limiting  reactant  is  60  percent  or  higher  at  the  jet  coalescence  plane 
(Vassilatos  and  Toor,  1965;  Mao  and  Toor,  1972;  Fisher,  1974) .  This 
is  a  significant  shortcoming  of  a  multitube  mixing  head.  A  more 
satisfactory  arrangement,  with  screens. is  found  in  the  flow  reactor 
of  Keeler  et  al.  (1965).  It  consists  of  a  specially  designed  "hair¬ 
brush"  injector  and  an  ordinary  wire  mesh  at  the  reactor  entrance 
(see  Figure  6).  Reactant  A  enters  the  reactor  through  the  main  flow 
and  B  enters  isokineticallv  through  the  "hairbrush"  injector.  The 
isokinetic  injection  implies  a  large  difference  of  flow  rates  between 
two  streams.  As  a  result,  the  low  flow  rate  stream  (i.e.,  B)  is 
proportionately  more  concentrated  to  render  the  two  streams  stoichio¬ 
metric.  The  injector  provides  a  uniform  coarse  scale  distribution  ot 
reactants  over  the  cross-section  with  only  a  little  micromixing 
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(hence,  little  reaction).  The  screen  downstream  of  the  injector  trips 
the  turbulence  and  thereby  induces  micromixing  and  reaction.  Keeler's 
experiments  showed  that  the  extent  of  conversion  at  the  jet  coalescence 
plane  was  merely  2  to  5  percent  (Keeler  et  al.,  1965).  A  similar 
arrangement  of  "hairbrush"  injector  and  screen  is  employed  in  the 
present  work. 

Many  single  phase  competing  reactions  of  the  type 

A  +  B  +  R  (6) 

R  +  B  S  (7) 

have  been  reportedly  found  to  be  sensitive  to  micromixing.  Bourne 
(1982)  has  compiled  a  list  of  such  reactions.  All  the  reactions  occur 
in  a  single  liquid  phase.  Among  these,  the  azo-coupling  of  1-naphthol 
(A)  with  diazotised  sulfanilic  acid  (B)  in  dilute  aqueous  solution  at 
room  temperature  is  considered  the  most  suitable.  In  this  scheme 
(Bourne  et  al.,  1981),  the  first  coupling  leads  to  formation  of 
4-(4-sulfophenvlazo)-l-naphthol  (R) ,  brief lv  described  as  monoazo  dye, 
and  the  subsequent  coupling  of  R  with  B  produces  2 , 4-bis (4-sulf ophenylazo) 
-1-naphthol  (S),  briefly,  diazo  dye.  Both  R  and  S  absorb  light  in  the 
visible  part  of  the  spectrum  (400  nm-600  nm)  and  their  concentrations 
can  be  measured  spectrophotometr icallv .  The  reaction  system  has  several 
attractive  features: 

(1)  The  reactions  are  irreversible  and  the  kinetics  are  rela¬ 
tively  well  understood. 

(2)  ,\mong  the  candidate  competing  reactions,  these  two  reactions 
have  the  fastest  kinetics  at  room  temperature.  The 


characteristic  reaction  time  scale  for  the  first  coupling 
is  0.5  msec  and  for  the  second  one  is  1  sec.  Computer 
simulations  have  shown  that  for  a  given  mixing  time,  this 
system  of  reactions  yields  the  largest  difference  among  the 
predictions  of  rival  models  as  compared  to  other  reported 
reaction  schemes.  This  fact  makes  these  reactions  particu¬ 
larly  suitable  for  our  purpose. 

(3)  The  reactions  are  conducted  in  an  aqueous  medium  (pH  =  10.00). 
This  is  convenient,  safe  and  inexpensive. 

_3 

(4)  The  low  reagent  concentrations  (0. 1-0.1  mol  m  )  are 
economically  attractive. 


B.  General  Description 


A  schematic  diagram  of  the  flow  system  used  in  this  study  is  shown 
3 

in  Figure  7.  A  0.21  m  polyethylene  tank  serves  as  the  feed  tank 
for  the  reactant  A.  The  solution  is  pumped  by  a  746  W  sealless. 


magnetic  drive  centrifugal  pump  (Model  TE-7R-MD,  March  Mfg.,  Inc.,  IL) 
through  a  rotameter  (Fischer  &  Porter  Co.,  PA,  precision  bore  flowrator 
tube  No.  B6-35-10/77)  to  a  flow  distributor.  The  passage  of  the  flow 
through  this  distributor  leads  to  the  development  of  a  nearly  flat 
velocity  profile  and  the  subsequent  passage  through  a  contraction 
nozzle  further  reduces  the  existing  turbulence.  The  flow  then  passes 
through  the  "hairbrush"  injector,  the  turbulence  producing  screen, 
and  enters  the  reactor.  The  reactor  is  simply  a  t).45  m  section  or 
5.08  cm  I.D.  Pvrex  pipe  which  terminates  in  a  stainless  steel  overflow 


chamber . 


Detector 


The  reactant  B  solution  is  pumped  from  a  2x10  m  laboratory  glass 
beaker  bv  a  5.4x10  ^  W  laboratory  metering  pump  with  a  micrometer  flow 
adjuster  (Fluid  Metering,  Inc.,  NY,  Model  RP-D)  into  the  plenum 
chamber  of  the  "hairbrush"  injector.  The  solution  then  enters  the 
main  flow  through  the  "hairbrush"  (hypodermic)  needles  in  the  form  of 
axial  jets.  The  flowrate  of  B  is  usually  l/100th  of  the  flow  rate  of 
A  but  feed  B  is  proportionately  more  concentrated  than  feed  A. 

The  needle  valve  controlling  the  naphthol  flow  rate  and  the 
three-way  valve  before  the  injector,  are  both  made  from  stainless 
steel.  All  other  auxiliary  parts  such  as  piping,  bends,  tees,  and 
unions  are  made  from  PVC.  The  piping  for  the  naphthol  stream  is  2.54  cm 
I.D.  PVC  and  that  for  the  diazotised  sulfanilic  acid  is  0.64  cm  I.D. 

PVC.  The  mechanical  vibrations  from  the  pumps  are  reduced  by  putting 
styrofoam  paddings  underneath  the  pumps  and  by  connecting  the  pump 
outlet  to  the  remaining  flow  loop  by  a  Tygon  tubing  of  the  same 
diameter.  Flexible  Tygon  tubing  further  damps  the  vibrations. 

Copper-Constantan  thermocouples  are  placed  at  appropriate  loca¬ 
tions  and  temperatures  of  the  two  feedstreams  are  monitored  with  a 
digital  thermometer  (Omega  Engineering,  Inc.,  CT,  Model  2176A). 

C.  Flow  Distributor  and  Contraction  Section 
A  honevcomb  type  flow  distributor  is  used  in  this  work  where  the 
cells  have  a  large  length-to-diameter  ratio  (L/D  =  33).  Such  honey¬ 
combs  have  been  used  successfully  in  water  tunnels  to  reduce  turbulence 
to  acceptable  levels  (Lumlev  and  McMohan,  1967).  The  honeycomb  has 
two  effects:  it  reduces  the  level  of  the  existing  turbulence,  and  it 
creates  additional  turbulence  of  its  own.  By  evenlv  distributing  the 


flow,  the  honeycomb  also  yields  a  nearly  flat  mean  velocity  profile. 

For  economy  and  convenience,  the  packing  of  plastic  drinking  straws 
(6  mm  I . D . )  into  a  honeycomb-like  matrix  is  used.  Such  a  "poor  man's 
honeycomb"  has  been  successfully  employed  in  low-speed  wind  tunnel 
research  (Loerhke  and  Nagib,  1972).  The  straws  are  packed  into  a 
15  cm  I.D.  Pyrex  pipe  and  held  in  place  by  adhesive  on  the  straw-wall. 
The  estimated  solidity  (fraction  of  pipe  area  which  is  not  open  to  the 
flow)  of  this  honeycomb  is  0.33.  For  any  turbulence  manipulator,  such 
as  this  honeycomb,  it  is  important  that  the  solidity  be  maintained 
below  a  critical  value  (usually  0.43),  otherwise  jets  emanating  from 
the  manipulator  coalesce  to  form  a  central  jet  which  is  highly 
undesirable  (Tan-Atichat  et  al.,  1982).  A  settling  chamber  about  0.5 
diameter  long  follows  the  honeycomb.  This  allows  the  mean  profiles 
emanating  from  individual  straws  to  merge  and  settle  to  a  more  uniform 
overall  profile.  This  section  is  known  to  improve  the  flow  at  the 
entrance  to  the  test  section  (Pope,  1947). 

Next,  the  flow  enters  the  contraction  section.  It  has  been  shown 
that  the  fluctuating  component  of  velocity  varies  inversely  as  the 
square  of  the  contraction  ratio  (Pope,  1947).  The  contraction  ratio  of 
nine  used  in  this  research  conforms  to  standard  water  tunnel  practice. 
The  shape  of  the  contraction  zone  is  not  critical.  In  general,  anv 
smooth  curve  is  satisfactory  (Pope,  1947;  Steele,  1951).  The  section 
is  about  one  upstream  diameter  long  and  tapered  very  graduailv  at  the 
downstream  end  so  that  the  flow  has  space  to  even  out. 


D.  Injector  Head  and  Turbulence  Producing  Screens 


The  "hairbrush"  injector  is  a  duplicate  of  the  one  used 
by  Keeler  (1965) .  It  consists  primarily  of  a  matrix  of  18  gauge 
stainless  steel  hypodermic  needles  (total  of  37  needles)  arranged  on  a 
9/32  in.  square  lattice  as  shown  in  Figure  8.  Diazotised  sulfanilic 
acid  enters  the  circumscribing  plenum  chamber  (10  cm  O.D. ,  5  cm  I.D.) 
which  distributes  it  among  its  seven  3.18  mm  I.D.  tubes.  These  tubes, 
in  turn,  further  distribute  the  flow  among  the  needles  they  support. 

For  easv  fabrication,  the  plenum  chamber  is  constructed  from  brass. 
However,  brass  reacts  with  diazotised  sulfanilic  acid  at  experimental 
concentration.  Consequently,  the  interior  of  the  plenum  chamber  is 
coated  with  protective  varnish  ordinarily  available  at  a  local  art 
supply  store  (Univar,  F.  Weber  Co.,  PA).  The  thin  protective  coating 
prevents  any  measurable  reaction  between  the  B  stream  and  the  plenum 
chamber  wall  (see  Figure  9  for  a  photographic  view  of  the  injector). 

Keeler  (1965)  has  reported  that  the  variation  in  injection  rates 
from  needle  to  needle  is  less  than  0.5  percent.  Cairns  and  Prausnitz 
(1960)  have  shown  by  calculations  that  at  a  distance  of  about  13  mm 
downstream  from  the  injector  there  is  only  9  percent  "ripple"  in  the 
radial  concentration  profile  if  an  inert  dye  is  injected  through  the 
needles.  Thus,  this  system  is  an  attractive  device  to  establish  a 
nearly  plane  concentration  front.  The  feed  tubes  and  the  needles  do 
generate  turbulence  but  it  is  modified  by  the  turbulence  generator 
screen  downstream  of  the  needles.  In  the  present  work,  a  screen  is 
placed  either  upstream  or  downstream  of  the  injector,  producing  differ¬ 
ent  mixing  conditions  depending  on  the  screen  placement. 


The  screens  are  ordinary  wire  mesh  screens  (Newark  Wire  Cloth  Co., 


NJ) .  They  are  chosen  to  generate  different  turbulence  scales  in  the 
flow.  The  screens  are  placed  about  7  mm  downstream  or  upstream  of  the 
needle  tips.  A  Lucite  piece  with  a  6.4  mm  recess  provides  the  seat  for 
the  screen  (see  Figure  10) .  The  screen  is  actually  sandwiched  between  two 
carefully  cut  rubber  gaskets.  The  Lucite  piece  and  the  gaskets  have 
the  same  I.D.  as  the  reactor  and  the  assembly  rests  on  another  rubber 
gasket  placed  on  the  downstream  face  of  the  injector.  The  relevant 
screen  dimensions  are  listed  in  Table  4. 

Table  4 

Screen  Characteristics 
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10.  Turbulence  p 


E.  Probe  Mount 


The  probe  mount  is  a  50  cm  long  vernier  slide  which  is  attached 
to  a  steel  base  plate  (see  Figure  11).  The  groove  in  the  base  plate 
permits  movement  of  the  slide  in  a  horizontal  plane.  Since  the  groove 
can  be  rotated  around  the  center,  the  probe  mount  essentially  allows 
three-dimensional  movement,  with  an  axial  positioning  accuracy  of 
0.25  mm. 


F.  Pulsed  Ultrasonic  Doppler  Velocimeter 
In  this  study  a  well  characterized  turbulent  flow  field  is  desired 
in  the  reactor.  Measurements  of  local  mean  and  fluctuating  velocity 
at  several  locations  are  necessary  for  such  a  characterization.  A 
pulsed  ultrasonic  Doppler  velocimeter  (PUDV)  is  used  for  these  measure¬ 
ments. 

The  PUDV  operates  on  the  Doppler  principle.  When  sound  waves  are 
reflected  off  a  moving  object,  there  is  a  shift  in  their  frequency. 

This  Doppler  shift  is  related  to  the  velocity  of  the  object  in  the 
following  manner: 

2uf  cos0 

Af  ,  =  - 2 -  (8) 

d  u 

s 

where  u  is  the  velocity  of  the  object,  Af^  is  the  Doppler  shift  in 
'  frequency,  ug  is  the  speed  of  sound  in  the  medium,  3  is  the  acute  angle 

I  between  the  sound  beam  and  the  velocitv  vector  and  f  is  the  source 

,  o 

■  frequency. 

i 

i 

■  The  acoustic  transducer  consists  of  a  single  unfocused  piezo¬ 


electric  crystal  (3  mm  diameter)  which  covers  the  tip  of  a  tungsten 
epoxy  stem  (Etalon  Corp.,  IN).  The  crvstal  transmits  10  >tUz  ultrasound 
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in  0.4  ysec  pulses  with  a  pulse  repetition  frequency  of  20  kHz.  The 
echoes  returning  from  scattering  particles  in  the  fluid  are  quadrature 
phase  detected  and  sampled  bv  0.3  psec  pulses  with  variable  delay 
between  3  and  51  ysec.  This  delay  provides  range-gating  between  2  and 
39  mm.  The  sample  volume  is  located  3.5  diameter  way  from  the  trans¬ 
ducer  face.  Thus,  the  sample  flow  field  is  expected  to  be  free  from 
transducer  induced  disturbances.  The  quadrature  audio  signal  is 
filtered  at  the  pulse  repetition  frequency,  and  low  frequency 
components  below  about  200  Hz  are  removed  prior  to  processing  by  a 
directional  zero  crossing  counter  which  converts  Doppler  frequency 
shift  into  voltage.  This  voltage  proportional  to  velocity  signal  is 
ultimately  smoothed  by  a  one-pole  filter  at  120  Hz.  The  maximum 
velocity  which  can  be  resolved  by  the  PUDV  without  aliasing  error  is 
75  cm/sec.  The  system  was  designed  and  built  by  Dr.  C.  J.  Hartley  of 
Baylor  Medical  College,  TX. 

The  sample  volume  is  approximately  cylindrical,  0.45  mm  in  the 
direction  of  transmission  and  3  mm  in  the  transverse  direction  with 
the  equivalent  volume  of  a  sphere  having  1 . 8  mm  diameter.  It  is 
important  to  consider  how  the  turbulence  scales  might  compare  x^ith 
these  dimensions.  One  would  expect  the  energy  containing  eddies,  which 
make  the  predominant  contribution  to  the  root-mean-square  velocity 


fluctuations  (u  ) ,  to  have  a  length  scale  (Z  )  of  order  of  the  size 
rms  e 


of  the  dominant  turbulence  generator  (6.5  mm)  which  is  discussed  later 
on.  Thus,  the  sample  volume  of  the  velocimeter  would  appear  to  be 
sufficiently  small  in  the  direction  of  transmission  to  resolve  the 
most  important  scales  of  the  turbulent.-.  The  frequence-  of  fluctuations 


in  the  energy  containing  eddies  (w  )  should  have  an  order  of  magnitude 


given  bv  (u  /£  ).  For  the  maximum  u  measured  within  the  reactor 
rms  e  rms 

( ~0 -  8  cm/sec),  is  of  order  of  2  Hz.  Thus,  the  frequency  response 

of  the  velocimeter  output  (120  Hz)  should  be  more  than  adequate  for  the 
large  scales. 

The  smallest  scales  in  the  flow  may  be  characterized  by  the 
Kolmogorov  microscales  of  length  (ri)  and  frequency  (co)  as  follows: 

n  =  (v  !z)  ,  co  =  (e/v)  ,  .  .  .  (9) 

where  v  is  the  fluid  kinematic  viscosity  and  e  is  the  kinetic  energy 

3 

dissipation  rate  per  unit  mass  (e  =  u  /£,  ) .  Taking  Z  and  u  as 

rr.;  s  s  6  rm  s 

2 

before  and  noting  V  =  .008  cm  /sec,  we  find  n  ~  0.3  mm  and  to  -  10  Hz. 
Thus,  the  frequency  response  of  the  velocimeter  output  (120  Hz)  is 
also  adequate  for  these  scales,  but  the  sample  volume  is  an  order  of 
magnitude  too  large. 

Scattering  of  ultrasound  is  provided  by  fine  air-bubble  mist  found 
in  the  ordinary  tap  water  (mean  bubble  diameter  x  60  u)  .  The  mean  rise 
velocity  of  the  bubbles  is  3  percent  of  mean  flow  velocity  in  the 
reactor.  Therefore,  the  bubbles  are  expected  to  easily  follow  the 
fluctuating  velocities  in  the  reactor.  Jorgensen  et  al.  (1973) 
describe  the  PUDV  and  Garbini  et  al.  (1982)  discuss  the  use  of  this 
ultrasonic  technique  in  measurement  of  fluid  turbulence.  The  veloci¬ 
meter  yields  an  output  signal  of  0.25  volts/KHz  shift  in  frequency. 

This  analog  signal  is  passed  through  a  low  pass  filter  (Model  3750. 
Krohn-Hite  Corp.,  MA)  to  an  integrated  data  acquisition  and  data 
processing  microcomputer  system  (Model  DATA  6000,  Data  Precision,  MA) . 
The  microcomputer  provides  mean  and  r.m.s.  values  of  the  detected 


signal  from  which  the  mean  and  r.m.s.  velocities  are  calculated  (.see 
Appendix  B  for  a  sample  calculation).  The  results  are  reproducible 
within  5  percent  for  the  mean  and  15  percent  for  the  r.m.s.  voltaces. 
with  ±  mV  accuracy. 


G.  Spectrophotometer  System 

The  primary  goal  in  this  work  is  to  measure  local  concentrations 
of  R  (monoazo  Hye)  and  S  (diazc  dye)  at  several  locations  in  the 
reactor.  Since  both  R  and  S  absorb  visible  light,  absorption  spectro¬ 
photometry  (colorimetry)  is  employed  to  measure  these  concentrations. 
This  technique  is  based  on  the  Beer-Lambert  law.  According  to  this 

law,  when  a  monochromatic  light  beam  of  constant  intensity  I  '  passes 

o 

through  a  small  volume  containing  an  absorbing  specie  i  with  concen¬ 
tration  C  and  over  a  length  £,  the  resultant  attenuated  light  intensit 
I  is  given  by 


l'  =  IA  exp(-eXC.Z)  (.10) 

o  11 

E^  is  the  molar  absorptivity  of  the  specie  i  at  the  wavelength  \ . 

For  our  purposes,  the  Beer-Lambart  law  is  assumed  to  hold  for  a 
stochastic  concentration  C  as  well,  where 


C.  =  C.  +  C’. 

ill 


The  overbar  denotes  the  time-mean  vaJuc  while  prime  denotes  the 
fluctuating  part.  Then,  the  attenuated  light  intensity,  I  ,  is  also 


stochastic  in  nature.  Now,  if  only  R  and  S  absorb  in  the  visible 
range,  and  do  so  independently  at  the  prevailing  concentration  levels 
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(of  order  10  “  to  10  mol/m  ) ,  one  can  write 


tA  -  IA  exP(-eAV  -  4cs2) 


(12) 


Alternatively, 

A 


In  (^)  £r"Cr  +  es^Cs 


(13) 


Taking  time-mean  values  of  all  the  stochastic  variables,  one 
arrives  at 


-r  IA  ,  _  ,  __ 

A  =  £n  (-T-)  =  eR(CR  +  es£Cs 


(14) 


where  AA  is  the  mean  absorbence  of  the  medium.  In  general, 


AA  =  £n  (IAef/IA)  =  4£CR  +  es^CS  (15) 

■\ 

where  IA  _  is  the  intensity  of  light  in  a  reference  medium.  This 
ret  _  • 

linear  relationship  of  A  with  C  and  C,  is  very  critical  to  measure- 

R  5 

ments.  A  standard  linear  regression  can  be  used  on  AA  values  at 

several  wavelengths  (but  at  the  same  spatial  location)  to  determine 
—  —  )  \ 

Cn  and  Cc ,  provided  £nl  and  Zcl  are  known.  The  latter  quantities  can 

Kb  Kb 

be  determined  separately  by  calibration  of  a  sample  of  known  concen¬ 
tration. 

It  turns  out  that  R  and  S  are  highly  absorbing  in  the  range  400- 

600  nm  (i.e.,  eA  and  z'),  have  relatively  large  magnitudes).  This 
R  S 

indicates  that  the  path  length  £  must  be  about  1  mm  if  Equation  (14) 
is  to  be  used  to  measure  concentration  levels  of  10  to  10  mol/m 
without  saturating  the  detector  (AA  should  lie  between  0. I  and  1.0). 
Thus,  an  intrusive  fiber-optic  probe  is  needed  to  measure  absorbance 
locally . 
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In  the  past,  fiber-optic  probes  have  been  successfully  used  to 
measure  concentration  fluctuations  in  chemically  nonreacti^g  (Lee  and 
Brodkev,  1963;  Xve  and  Brodkey,  1967)  and  reacting  (Mahalingam  and 
Chevray,  1983)  flows.  For  greater  precision,  all  these  studies  have 
employed  their  own  designs  of  optical  and  electronic  measurement  system. 
Since  measurement  of  mean  concentrations  is  of  primary  interest,  which 
is  considerably  less  demanding  than  measuring  fluctuating  concentra¬ 
tions,  an  off-the-shelf  colorimeter  (Brinkman,  Model  PC-801)  is  used 
in  conjunction  with  a  specially  designed  fiber-optic  probe  for  the 
absorbance  measurements.  The  design  of  the  probe  is  based  on  the 
requirements  that  the  probe  introduce  only  a  minimum  of  flow  disturbance 
in  the  measurement  volume  and  that  it  provide  an  acceptable  spatial 
resolution.  These  and  other  fabrication  requirements  have  led  us  to 
the  fiber-optic  probe  shown  schematically  in  Figure  12.  The  probe  was 
fabricated  according  to  our  specifications  by  Fiberoptics  Fabrications, 
Inc.,  MA . 

The  transmitter  fiber  bundle  has  a  1  mm  O.D.  and  is  illuminated 
by  the  colorimeter  light  source  at  one  end.  The  receiver  fiber  bundle, 
which  aiso  has  a  1  mm  O.D. ,  is  separated  by  1  mm  from  the  end  of  the 
transmitter  bundle  (see  bottom  of  Figure  12)  and  carries  the  attenuated 
signal  back  to  the  colorimeter.  The  signal  passes  through  an  inter¬ 
ference  filter  mounted  on  a  filter  wheel  before  impinging  on  the  silicon 
photodetector.  The  filter  wheel  has  six  different  filters  mounted  on 
it  which  allows  selection  of  six  different  wavelengths  (250,  -<0,  -.90, 
520,  325,  and  570  nm) ,  with  a  20  nm  half-peak  bandwidth  for  each  filter. 
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intensity,  which  is,  after  passage  through  a  variety  of  electronic  cir¬ 
cuits  (somewhat  modified  bv  us),  available  as  an  output  voltage. 

A 

The  output  voltage  is  directly  proportional  to  tn(lVl  ),  which  is 
different  from  =  [Jln(I^/l^)  ]  .  The  former  involves  first  taking  the 
time  average  of  a  signal  and  then  taking  the  natural  logarithm  of  it 
while  it  is  exactly  the  other  way  around  for  the  latter.  Fortunately, 
even  very  conservative  estimates  show  that  the  error  introduced  by 
equating  these  two  quantities,  i.e., 

i.n(IA/lX)  AX  (16) 

o 

is  less  than  1  percent  in  the  experimental  range  (see  Appendix  C) . 

Thus,  we  are  able  to  get  an  output  signal  (V“)  which  is  directly 
proportional  to  the  absorbance  (AX) . 

The  Brinkman  colorimeter  is  chosen  because  it  has  several 
attractive  operating  characteristics: 

(1)  It  supplies  phase  shifted,  AC  modulated  light  to  the  probe 
tip.  At  the  detector  end,  it  uses  an  amplifier  with  an 
electronic  chopper  that  is  synchronized  with  the  light  source, 
but  phase-shifted  bv  90°;  as  a  result,  ambient  light  does  not 
affect  the  instrument  reading. 

(2)  It  has  electronic  circuitry  that  ensures  excellent  line 
stability.  Other  laboratory  appliances  such  as  stirring 
motors  and  pumps,  plugged  into  the  same  circuit,  do  not  affect 
the  stability  of  the  colorimeter  output. 

(3)  It  exhibits  good  linearitv  in  the  sense  of  Beer-Lamburt  law 
under  the  experimental  conditions. 


-  -  -  - 


It  is  important  to  note  that  the  Brinkman  colorimeter  is  a  single 
beam  spectrophotometer.  Since  the  electronic  chopper  at  the  detector 
end  is  electrically  synchronized  with  the  light  source,  1^  is  an 
arbitrary  internal  reference  which  is  always  a  constant.  For  an 
absorbance  value  with  respect  to  another  reference  medium  (1^  ^),  two 
separate  measurements  are  required.  The  basic  method  of  analysis  is  as 
follows.  For  a  given  solution  containing  R  and  S, 
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where  k  is  the  proportionality  constant. 
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Similarly,  for  a  reference 


solution, 
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Subtracting  (17)  from  (16),  we  get 
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where  k  =  k  and  k  =  k  cA£.  Both  k„  and  k,,  are  constants  and  are 

R  p  R  p  S  k 

independently  determined  bv  calibration  of  the  probe  with  known  solu¬ 
tions  of  R  and  S.  Thus,  AVA  values  are  needed,  each  of  which  involve 
two  measurements  (\A  and  V“  ,)  at  each  location,  to  determine  C  and 
C;,.  A  standard  linear  regression  routine  is  used  to  determine  the 

values  of  C  and  C  which  will  minimize  the  square  of  the  difference 
R  S 

between  the  measured  A\A  at  five  wavelengths  and  the  calculated  values 

of  AVA  (based  on  (L  and  C„) .  These  values  of  Cn  and  C„  are  termed  the 
R  S  R  s> 

"experimental  values."  In  general,  the  standard  deviations  of  these 
regressed  values  are  10  and  17  percent  for  CR  and  Cg,  respectively. 

The  gap  size  of  the  fiber-optic  probe  is  small  (1  mm)  and  it 
would  seem  that  the  boundary  layer  along  the  exposed  fiber-faces  may 
affect  the  local  flow  field  in  the  gap.  Based  on  the  fact  that  Lee 
and  Brodkev  (1963)  did  not  find  any  such  effects  in  concentration 
fluctuation  measurements  with  a  similar  probe  (1  mm  tip  diameter,  1  mm 
gap),  it  is  expected  that  the  boundary  layer  growth  has  negligible 
effect  on  the  mean  concentration  measurements  of  this  studv. 


V .  MEASUREMENTS 


The  experimental  work  was  divided  into  two  parts:  Flow  Field 
Characterization,  and  Mixing-Reaction  Studies. 

A.  Flow  Field  Characterization 

1 .  Comments 

The  flow  field  characterization  involved  measurement  of  the  mean 
and  r.m.s.  velocities  at  several  locations  in  the  reactor.  These 
experiments  were  conducted  separately  from  the  mixing-reaction  studies 
but  retained  the  same  hydrodynamic  conditions  (i.e.,  the  same  pipe 
Reynolds  number) . 

Two  points  deserve  comments  here. 

(1)  It  is  generally  necessary  to  seed  the  flow  with  scattering 
particles  when  a  PUDV  is  used.  Polystyrene  microspheres  at 
0.5  wt%  concentration  perform  satisfactorily  in  such  cases. 
Nevertheless,  it  was  difficult  to  maintain  a  uniform  sus¬ 
pension  of  microspheres  in  the  feed  tank.  This  was  mainly 

3 

due  to  the  large  hold-up  (0.05  m  )  and  the  inability  of 
the  stirrer  to  sweep  the  particles  from  the  floor  of  the  feed 
tank.  It  was  found,  per  chance,  that  the  fine  mist  of  air 
bubbles  in  ordinary  tap  water  obviated  the  need  for  any 
further  seeding  of  the  flow.  Therefore,  ordinary  tap  water 
was  used  as  the  working  fluid.  The  validity  of  the  experi¬ 
mental  procedure  is  considered  in  the  next  section. 

(2)  The  PUDV  and  the  data  acquisition  system  initially  were  found 
to  electrically  interact  with  each  other  resulting  in  a  highly 
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noisy  and  nonreproducible  signal.  A  specially  designed  low 
pass  filter,  which  filtered  both  on  the  "hot"  and  the  ground 
sides  of  the  connecting  cable  at  120  Hz,  along  with  a  large 
length  of  the  cable,  was  found  effective  in  reducing  the 
interaction  to  an  acceptable  level.  The  measured  values  con¬ 
tained  this  small  amount  of  electronic  noise.  Furthermore, 
the  Doppler  signal  in  this  technique  is  always  contaminated 
with  a  broad-band  random  noise  called  the  "Doppler  ambiguity." 
This  phenomenon  is  similar  to  that  occurring  in  laser 
anemometry,  but  its  effect  is  very  significant  for  this  tech¬ 
nique  because  of  a  much  larger  sample  volume.  The  random 
noise  is  caused  by  the  random  transit  of  scattering  particles 
through  the  finite  sample  volume.  There  is  also  noise 
associated  with  buoyancy  of  bubbles  and  the  bubble  size 
distribution.  All  except  the  Doppler  ambiguity  were  considered 
in  the  subsequent  calculation  cf  velocities  in  the  reactor. 

2 .  Validation  of  the  Experimental  Method 

To  establish  the  validity  of  our  experimental  procedure,  a  test 
case  of  fully  developed  turbulent  flow  in  a  pipe  was  chosen.  The 
principal  advantage  of  this  configuration  was  that  a  large  body  of 
experimental  results  using  hot-film  and  laser-Doppler  anemometry  was 
available  in  the  literature.  Comparison  of  our  data  with  the  reported 
data  provided  a  convenient  check  on  our  method.  Tine  mean  and  r.m.s. 
velocities  as  well  as  the  turbulence  spectra  at  the  centerline  or  the 
pipe  were  compared. 

A  schematic  diagram  of  the  test  flow  system  is  shown  in  Figure  13. 
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Ic  utilized  all  the  parts  of  the  original  flow  system  except  that  a  1.8 

cm  I.D.  pipe  was  used  in  place  of  the  reactor  assembly.  Tap  water  was 

3 

added  to  the  0.21  m  polyethylene  tank  continuously.  It  was  then 
pumped  through  the  rotameter  to  the  base  of  the  vertical  pipe, 
eventually  to  the  drain.  The  ultrasound  transducer  was  placed  at  an 
angle  0  =  43°  to  the  flow  direction.  The  transducer  emitted  a  3  mm 
wide  unfocused  beam  at  10  MHz,  which  had  a  range  of  2  to  35  mm  but 
began  diverging  after  20  mm.  For  experimental  purposes,  the  range  was 
adjusted  to  locate  the  measurement  volume  at  the  centerline.  This  was 
achieved  by  finding  the  location  where  the  PUDV  output  was  maximum. 

The  phasic  (raw)  signal  of  the  PUDV  was  passed  through  the  Krohn- 
Hite  filter  whose  cut-off  frequency  was  set  at  120  Hz,  and  also  through 
the  special  passive  filter  before  it  was  fed  to  the  microprocessor. 

The  sampling  period  was  2  ms  and  the  signal  was  continually  recorded 
and  processed.  The  data-acquisition  was  stopped  when  the  mean  and  the 
r.m.s.  values  of  the  entire  record  became  stable.  Typically  this 
corresponded  to  a  real  time  record  of  about  10  sec.  For  spectral 
measurements,  200  records  with  512  points  per  record  and  frequency 
smoothing  (Hanning;  were  ensemble  averaged  to  get  a  stable,  reproducible 
spectrum. 

The  microprocessor  gave  mean  and  r.m.s.  values  in  mV.  These  were 

appropriately  converted  to  proper  units  of  velocity  bv  using  Equation 

1 

(B.l)  of  Appendix  B.  The  power  spectrum  was  £iven  in  terms  of  mV“/Hz  as 
a  function  of  frequency  in  Hz.  This,  too,  was  appropriately  scaled  to 
compare  with  other  published  results.  Computational  details  are  given 
in  Appendix  D. 
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The  raw  data  are  presented  in  Table  5.  The  final  comparisons  are 
shown  in  Figures  14  and  15.  The  turbulence  intensity  measurements  fall 
in  the  same  range  as  the  reported  values  as  seen  in  Figure  14. 

Similarly,  spectral  measurements  at  pipe  Reynolds  number  of  5900  follow 
data  of  Garbini  et  al.  (1982)  very  closely.  Thus,  the  comparisons  were 
very  favorable.  They  provided  a  reasonable  basis  for  the  use  of  the 
present  method  in  investigating  the  fluid  mechanics  of  the  actual 
reactor. 

3 .  Experimental  Procedure 

Unlike  in  the  previous  verification  test,  for  characterization  of 
the  reactor  the  ultrasound  transducer  was  inserted  vertically  into  the 
flow  field  (8  =0°)  as  shown  previously  in  Figure  7.  This  allowed 
measurement  at  any  desired  location  in  the  flow  field  but  required 
special  precaution  in  the  data  collection.  In  this  configuration,  the 
bubble  mist  impinged  directly  on  to  the  transducer  face.  The  tiny 
bubbles  had  a  tendency  to  cling  to  the  transducer  face  and  after  a 
short  time  they  coalesced  to  form  a  much  bigger  bubble  which  distorted 
the  signal  considerably  and  caused  it  to  drift.  It  was  observed, 
however,  that  a  freshly  cleaned  transducer  face  gave  reproducible 
results  in  the  period  prior  to  the  drift.  Furthermore,  although  the 
transducer  was  inserted  into  the  flow,  the  actual  measurement  volume 
was  kept  about  L  cm  in  front  of  the  transducer  face.  Thus,  effectivelv, 
this  technique  did  not  disturb  the  flow  field  to  any  significant 
extent . 


An  experiment  was  conducted  in  the  following  sequence: 
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Table  5 

Centerline  Turbulence  Intensities  in  a  Straight  Pipe 
(Pulsed  Ultrasound  Doppler  Velocimeter  Data) 
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Figure  14.  Turbulence  intensity  measurements  [Hatched 

area — values  reported  in  literature  (Carbini 
et  al.,  1982)] 
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Formalized  turbulence  spectral  densities 
[( — )  ultrasound,  and  (-  -  -)  tint  t  i  In  data 
of  Far b ini  et  al.  (1982)] 


Figure  15. 


(1)  The  PUDV  and  the  low  pass  filter  were  turned  on  to  allow  for 
approximately  half  an  hour  warm-up  time. 

(2)  Ordinary  tap  water  and  another  stream  of  hot  water  were  con- 

3 

tinuouslv  added  to  the  0.21  m  tank  (A  stream).  The  height 

and  the  temperature  of  the  liquid  were  controlled  by  manual 

3 

adjustments.  Typically  the  tank  had  a  0.04  to  0.05  m 
hold-up  and  the  temperature  was  25-27°C. 

(3)  The  flow  rate  of  the  bubbly  water  through  the  rotameter  was 
adjusted  bv  the  needle  valve  and  the  globe  valve  in  the 
recycle  loop.  As  a  precaution,  the  rotameter  was  always 
calibrated  prior  to  each  run. 

(4)  Ordinary  distilled  water  was  used  as  the  working  fluid  for 
the  secondary  (B)  stream.  It  did  not  contain  any  air  bubbles. 
The  flow  rate  of  this  stream  was  adjusted  by  the  micrometer 
attached  to  the  metering  pump.  This  also  was  always  cali¬ 
brated  prior  to  each  run.  The  distilled  water  was  period¬ 
ically  replenished  throughout  the  run  to  ensure  a  continuous 
supply  to  the  "hairbrush"  injector. 

(5)  First,  the  PUDV  was  internally  calibrated.  The  output 
voltage  from  the  phasic  jack  was  adjusted  to  positive  and 
negative  0.610  volts  when  the  polarity  switch  on  the  front 
panel  of  the  instrument  was  kept  at  the  positive  and  negative 
positions,  respectively.  The  transducer  was  unplugged  during 
the  calibration. 

(6)  Next,  the  transducer  was  plugged  into  the  PUDV.  Since  the 
pulse-repet  it ion-f requencv  (PRF)  of  PUDV  was  19.53125  KHz. 
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the  calibration  for  the  range  control  was  0.3125  volts/cm. 

The  range  was  adjusted  accordingly  to  adjust  the  measurement 
volume  at  about  1  cm  in  front  of  the  probe. 

(7)  The  instrument  was  switched  from  the  calibration  (CAL)  to  the 
directional  (DIR)  mode,  the  polaritv  switch  to  the  "toward" 
position — since  the  flow  moved  towards  the  transducer — and 
the  filter  was  set  at  the  phasic  position.  The  sensitivity 
control  was  set  approximately  midway  so  that  the  background 
audio  noise  just  barely  triggered  the  decoder. 

(8)  The  ultrasound  probe  was  hand-held  in  air  to  measure  anv 
background  electronic  noise. 


the  desired  location. 


real  time  record.  By  this  time,  bubble  accumulation  on  the 
transducer  tip  began  to  distort  the  signal  and  hence  data 
collection  was  stopped.  The  measured  mean  and  r.m.s.  voltages 
were  recorded  from  the  microprocessor  screen  manually. 

(11)  The  transducer  face  was  wiped  gently  with  i  tissue  paper  and 
again  lowered  into  the  flow  field  for  another  measurement. 

The  entire  flow  loop  was  purged  frequently  with  a  high  flow  rate 
pulse  to  flush  out  most  of  the  bigger  air  bubbles  which  accumulated  in 
the  loop.  Otherwise,  such  big  bubbles  passed  through  the  measurement 
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volume,  making  that  measurement  unacceptable.  Nevertheless,  with 
proper  care  and  patience,  it  was  possible  to  obtain  reproducible  data. 

The  pulsed  ultrasound  technique  was  further  used  to  characterize 
the  scattering  bubble  mist  in  terms  of  a  mean  bubble  rise  velocitv. 
Since  the  rise  velocity  Ls  related  to  bubble  diameter,  it  is  also  a 
direct  measure  of  the  mean  bubble  diameter.  The  random  noise  asso¬ 
ciated  with  bubbles  was  also  obtained  from  the  same  measure- 
ne  n  t . 

At  the  end  of  each  run,  the  flow  was  suddenly  stopped  and  a 
freshly  wiped  transducer  was  introduced  into  the  fluid.  The  mean  and 
the  r.m.s.  voltages  measured  in  this  fashion  directly  corresponded  to 
the  two  quantities  of  interest — bubble  rise  velocity  and  the  random 
noise. 

4 .  Experimental  Conuit ions 

Mean  and  r.m.s.  velocities  were  measured  in  the  reactor  at 
several  axial  locations  on  the  centerline  and  also  a  few  radial  loca¬ 
tions.  each  measurement  was  repeated  at  least  twice  to  check  its 
renroducibilitv. 

Four  different  experimental  conditions  were  investigated  corre¬ 
sponding  to  -'our  different  placements  of  the  turbulence  generating 
screen  and  the  "hairbrush"  injector.  in  case  1,  the  honovcomb  section 
was  direct!-/  connected  to  the  reactor  pine  and  the  injector  and  the 
screen  were  not  u.-,ed  at  all.  Hi  is  allowed  measurement  of  details  o: 
t  he  !  low  from  the  downstream  section  <>f  the  honevcom.h .  ill  the  wav  to 
the  middle  >.  t  the  reactor  pipe.  In  the  second  case,  -niv  the  "h  t  ir- 
brush"  injector  was  introduced  and  no  screens  were  used.  I'll  is  was  to 
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1.  Material.' 


Distilled  water  at  room  temperature  was  used  as  the  working  iluid 
in  our  mixing-reaction  experiments.  Sodium  carbonate  and  sodium 
bicarbonate  ("Baker  analyzed,"  J.  T.  Baker  and  Co.,  XJ)  were  used  as 
buffer  salts.  Recrvstall ized  1-nanhthol  ("Certified,"  risher  Scientific 
('i.u,  XJ)  and  d  Lazoticed  sulfanilie  acid  ("puruas  p..i .  , "  Fiuka  Chemical 
Corn.,  XT) — reactants  A  and  B — were  both  reagent  grade  pure.  The 
latter  was  available  as  a  moist  caste  from  the  manufacturer  because  it 
is  explosive  when  dr’/.  It  was  always  stored  in  a  rei  r  iterator .  1  he 

actual  diazotised  sulfanilie  acid  content  oi  the  paste  was  supplied 


b6 

issess  the  effect  of  the  injector  on  the  turbulence  intensitv  Level. 

In  Case  3,  the  largest  screen  (M  )  './as  placed  just  I  mm  upstream  of  the 
injector,  to  evaluate  the  effect  of  this  configuration  on  the  turbu¬ 
lence.  Finally,  in  Case  4,  the  same  screen  was  placed  /  mm  downstream 
from  the  injector. 

The  presence  of  large  bubbles  in  the  system  and  their  tendency  to 
accumulate  at  the  screen  in  the  case  of  the  smaller  screens  (M2  and  M^) 
prevented  us  from  measuring  the  effect  ol  the  smaller  screens  on  the 
turbulence . 

The  raw  data  are  listed  in  Appendix  E,  along  with  the  data  for 
bubble  rise  velocity  and  the  random  noise.  The  latter  values  are  used 
to  correct  mean  and  r.m.s.  velocities.  The  corresponding  calculated 
values,  although  not  reported  numerically,  are  plotted  and  discussed 
in  the  next  chapter. 

3.  Mixing-Reaction  Studies 
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Che  manufacturer  as  actual  lot  analysis.  Typically  the  naste  contained 
about  10  percent  eater. 


2 .  Calibration  of  tire  Fiber-Ontie  Probe 

Our  technique  was  based  on  Equation  (19)  of  the  previous  chapter: 

av'  .  ^ 

The  determination  of  and  was  done  through  separate  calibrations 
of  standard  solutions  of  monoazo  dye  (R)  and  diazo  dye  (S). 

Due  to  its  use  as  an  indicator  in  microscopy,  R  was  commercially 
available  as  a  very  pure  powder  (’’Standard  Fluka,"  Fluka  Chemical  Corp., 
NY).  Thus,  preparing  a  scandard  solution  of  R  in  alkaline  water  of 
pH  =  10.00  was  a  simple  matter. 

A  standard  solution  of  S  was  prepared  by  reaction  of  a  standard 
solution  of  R  with  an  aqueous  solution  of  diazotized  sulfanilic  acid. 
Essentially,  the  second  azo  coupling  reaction  R  +  il  •-  S  was  carried  out 
in  a  beaker  under  similar  conditions  of  pH  and  ionic  strength  as  in 
reactor  and  was  allowed  to  go  to  completion.  By  adding  a  slight  excess 
of  P,  it  was  ensured  that  no  unreacted  R  remained  in  the  solution. 

The  preparation  of  primary  standards  in  the  above  manner  was  a 
departure  from  the  procedure  followed  by  Kozicki  (1980).  Kozicki’s 
method  synthesized  pure  R  by  dissolving  L-naphthOl  in  ethanol  and  then 
reacting  it  with  a  concentrated  slurrv  of  diazotised  sulfanilic  acid 
at  room  temperature.  Since  R  was  available  in  very  high  nuritv 
commercially,  the  commercial  R  was  preferred.  Regarding  synthesis  .>: 
pure  S,  it  was  not  possible  to  separate  S  in  a  solid  lorm  I rom  the 


solution  when  Kozicki ' s  procedure  was  followed.  Similar  dill  mult"  was 
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encountered  bv  other  researchers  as  well  (I-i.  1984).  Thus,  we  were 
compelled  to  use  a  solution  of  S  for  the  primary  standard.  In  this 
context,  reacting  R  with  B  to  form  S  in  a  buffer  solution  was  a 
reasonable  way  to  prepare  the  standard. 

The  molar  absorption  coefficients  of  these  solutions  are  plotted 

as  a  function  of  wave  length  in  Figures  16  and  17.  For  comparison, 

values  obtained  by  Kozicki  (1980)  and  Li  (1984)  are  also  displayed  on 

the  same  plots.  In  most  cases,  values  for  this  study  are  20-25 

R 

percent  lower  than  Kozicki' s  values  while  Li's  values  are  3-8  percent 
lower  than  the  latter. 

Li  (1984)  used  R  synthesized  by  Kozicki' s  method  but  in  the 
present  study  commercially  available  high  purity  R  was  used.  Hence, 
the  difference  in  values  is  attributed  to  relative  purity  of  R  and 
errors  involved  in  making  solutions.  On  the  other  hand,  the  e’  values 
for  this  study  are  47-56  percent  lower  than  Kozicki 's  values  and  Li's 
values  also  are  10-18  percent  lower.  Kozicki  (1980)  reported  that  the 
elemental  analysis  of  R  conformed  within  0.8  percent  of  the  theoretical 
values  but  that  of  S  differed  by  as  much  as  17  percent.  Thus, 

Kozicki 's  S  either  had  significant  impurities  or  even  a  somewhat 
different  molecular  structure.  Since  the  method  of  making  primary- 
standard  of  S  in  the  present  work  mimics  the  reactor  conditions,  it  is 
considered  more  reliable. 

In  each  case,  five  different  solutions  of  varying  concentrations 
in  the  range  of  0.07  to  0.35  mol/m^  were  prepared  by  suitably  diluting 
the  standard  solution  with  an  alkaline  buffer  of  l.’n  ,CO  j/N'aHCO  ^  (pH  = 
10.00).  The  same  buffer  was  used  as  the  reference  solution.  The 
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colorimeter  output  was  measured  first  bv  dipping  the  fiber-optic  probe 
in  the  solution  and  in  the  reference  solution  immediately  thereafter. 

All  solutions  were  made  in  duplicate  to  improve  the  accuracy.  As 
such,  ten  AV“  values  were  obtained  at  each  wavelength,  two  values  at  each 
of  the  five  concentrations  of  R  (or  S) .  Standard  least  square  linear 
regression  of  these  ten  values  yielded  (or  k^)  (see  Figures  18  and 
19). 

3 .  Method  of  Analysis 

In  the  mixing-reaction  experiments,  AV“  was  measured  at  five 
different  wavelengths  at  each  location  in  the  reactor.  When  these 
values  were  substituted  in  Equation  (19),  along  with  calibration  con¬ 
stants  kjj  and  k^,  a  set  of  five  equations  involving  only  two  unknowns 
K.  S 

(C  and  C  J  was  obtained. 

R  S 

A  standard  linear  regression  package  (MINITAB)  was  used  to  find 

the  best  fit  values  of  C  and  C_.  These  best  fit  values  are  termed 

K  5 

the  "experimeuta1  ’'alues"  in  this  thesis. 

This  method  of  calculating  concentrations  of  R  and  S  was  tested 
with  sample  solutions  containing  known  concentrations  of  R  and  S.  The 
concentration  range  covered  was  typical  of  the  range  expected  in  most 
of  the  mixing-reaction  experiments.  A  comparison  of  actual  and  mea¬ 
sured  values  is  displayed  in  Table  b.  It  is  evident  that  in  all  but  one 
case,  the  measured  concentrations  are  within  12  percent  of  the  actual 
ones.  This  was  considered  satisfactory  in  view  of  the  Low  level  or 
concentrations . 

A  similar  comparison  is  presented  in  Table  7  for  a  set  of  solu¬ 
tions  where  S  is  present  in  very  low  concentration.  It  is  clear  that 


mi 


Data  fur  Evaluation  of  Mu t nod  of  Analysis  for  Solutions  with  Low 


the  method  may  yield  negative  values  of  Cc  in  some  cases  and  the 
relative  error  in  measuring  Cg  is  fairly  high.  This  seems  unsatis¬ 
factory.  However,  it  may  be  argued  that  since  S  is  present  in  verv 
low  concentration,  it  is  reasonable  to  ignore  the  presence  of  S  in 
these  solutions  and  use  AVA  values  to  determine  C  only.  The  regression 
equation  then  becomes 

AvX  =  kRC^  (20) 

cs=o 

The  last  two  columns  of  Table  7  show  that  CL  values,  calculated  thus, 

R 

are  within  15  percent  of  the  actual  C  values,  except  when  CL  <  0.08 

L\  R 

3 

mol/m  .  Generally,  this  method  gives  closer  agreement  with  the  actual 
value  of  CR  when  CR  is  at  least  15  times  higher  than  Cg .  As  will  be 
discussed  later,  this  approximate  method  was  used  to  measure  the 
concentration  of  R  in  runs  with  a  low  ratio  of  B  to  A  in  the  feed. 

1 .  Characterization  of  Feedstreams 

The  primary  feedstream  in  the  mixing-reaction  experiments  contained 
l-nanhthol  (A)  in  an  aqueous  buffer  solution  of  Na^CO^/NaHCO^  (pH  = 

10.00  ±  0.02) . 

The  same  amounts  of  distilled  water  and  L-naphtlol  were  used  for 

each  experiment.  Thus,  the  unmixed  feed  concentration  C  ^  was  always 

3 

about  0.3  mol/m  .  The  exact  value  was  determined  bv  measuring  the 
I'Y-absorbance  of  the  solution  at  245  nm.  A  standard  solution  of 
L-naphthol  in  an  identical  buffer  was  used  to  determine  L'Y-absor bance 
at  245  nm  as  a  function  of  concentration.  This  was  used  as  calibration 
for  the  feedstream  A.  A  UY-visible  spectrophotometer  (Model  8450A, 


Hewlett-Packard  Corp.,  CA)  was  used  for  this  purpose.  It  was  important 


to  use  micropipettes  to  prepare  solutions  of  appropriate  dilution, 

usually  one  part  solution  in  sixteen  parts  buffer,  to  get  precise 

results.  It  was  also  necessary  to  use  quartz  couvettes  in  all  these 

measurements.  Also,  the  spectrophotometer  was  allowed  to  warm  up  for 

at  least  a  period  of  20  minutes  to  ensure  reproducibility.  A  sample 

calculation  is  shown  in  Appendix  F. 

The  secondary  stream  was  a  relatively  concentrated  solution  (10- 

3 

30  mol/m  )  of  diazotised  sulfanilic  acid  in  water.  The  actual  lot 
analysis  supplied  by  the  manufacturer  was  used  to  calculate  the  amount 
of  paste  required  to  prepare  a  solution  of  required  concentration. 
Usually  it  took  about  10  minutes  of  stirring  with  a  stirrer  rod  to 
dissolve  the  entire  amount  of  acid  into  the  water.  The  resultant 
solution  was  yellowish  brown  in  color. 

5 .  Experimental  Procedure 

The  nature  of  our  experimental  design  was  such  that  at  least  two, 
preferably  three,  persons  were  required  to  carry  out  a  mixing-reaction 
experiment.  The  actual  duration  of  an  experiment  was  about  20  minutes 
but  altogether  the  experiment  required  24  hours  for  the  preparation. 
The  protocol  outlined  below  had  evolved  during  a  series  of  preliminary 
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experiments  conducted  primarily  to  establish  an  appropriate  procedure: 

(1)  The  fiber-optic  probe  was  calibrated  first.  Standard  solu¬ 
tions  of  R  and  S  were  prepared  and  the  calibration  was 
carried  out  in  the  manner  described  earlier. 

(2)  Next,  the  flow  svstem  was  flushed  with  fresh  tap  water 


several  times  to  remove  anv  adhering  dust  particles  1 rom  the 


m 


3® 


,  .JtJ 


system.  The  system  was  drained  at  the  end  of  the  flushing, 
a 

(3)  The  0.21  m  tank  was  filled  with  distilled  water.  Through 

a  hose  pipe  connected  to  an  available  distilled  water  supply 

3 

line.  The  flow  loop  had  a  hold  up  of  about  0.013  m  . 

By  flushing  the  system  once  with  a  little  distilled  water, 
most  of  the  remaining  tap  water  was  pushed  out. 

(4)  In  a  small  mortar,  1-naphthol  was  ground  with  a  pestle  to  a 
very  fine  powder.  About  9.7  g  of  that  powder  was  weighed  on  a 
weighing  paper. 

3 

(5)  Next,  the  stirrer  in  the  0.21  m  tank  was  turned  on  and 
1-naphthol  powder  was  slowly  added  to  the  swirling  content  of 
the  tank  through  a  small  hole  on  the  tank  top.  It  was 
important  to  set  a  proper  stirrer  speed  by  a  variac,  as  a 
higher  speed  might  only  make  the  swirling  faster  but  not 
disperse  the  powder  effectively.  The  stirring  was  continued 
overnight,  altogether  for  a  period  of  12  hours.  That  much 
time  was  sufficient  to  dissolve  the  entire  amount  of  1-naphthol 
added.  Such  a  long  time  was  necessitated  by  the  low  solubilitv 
of  1-naphthol  in  pure  water  at  room  temperature.  It  was 
possible  to  speed  up  the  dissolution  process  by  adding  the 
buffer  salts  beforehand  but  it  also  Led  to  slight  coloration 

of  uhe  fluid  after  a  few  hours.  The  coloration  was  attributed 
to  i-naphthol's  sensitivity  to  Light  in  an  alkaline  environ¬ 
ment,  and  was  deemed  undesirable  tor  our  purposes. 

(5)  When  the  aqueous  solution  of  1-naphthol  was  readv,  the  buffer 
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aits  were  added  to  it  in  appropriate  amounts  (typically  235 
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of  'Ta^CO-j  and  190  ±  5  g  of  NaHCO^)  .  The  salts  dissolved  verv 
rapidly  and  brought  pH  to  10.00  +  0.02.  That  also  corre¬ 
sponded  to  an  approximate  ionic  strength  of  10  g  ion/m^.  The 
pH  was  measured  with  a  standard  pH  meter  (Orion  digital 
ionalvzer/501,  Orion  Corp.,  MA) . 

(6)  Next,  the  colorimeter  and  the  low  pass  filter  were  turned  on. 
The  cut-off  frequency  was  set  at  0.25  Hz. 

(7)  The  1  hp  pump  was  turned  on  and  the  alkaline  naphthal  solution 
was  recycled  through  the  loop  back  to  the  feed  tank.  The 
recycling  dislodged  any  air  bubbles  in  the  honevcomb  and  also 
removed  air  pockets  at  the  corners  of  the  flow  loop.  It  also 
brought  the  entire  loop  to  a  thermal  steady  state. 

(8)  The  rotameter  for  the  mainstream  was  then  calibrated  with  a 
measuring  cylinder  and  a  stopwatch. 

(9)  Next,  the  mainflow  was  stopped  and  the  needle  valve  was 
closed  so  that  the  hydrostatic  head  at  the  tip  of  the  needles 
of  the  "hairbrush"  injector  was  maintained.  Under  those 
conditions  the  secondary  stream  metering  pump  was  calibrated 
by  pumping  distilled  water  against  the  head.  Once  again  a 
small  measuring  cylinder  and  a  stopwatch  were  used  in  the 
cal ibration. 

10)  Thirty-five  hundred  cc  of  distilled  water  was  added  to  a 

separate  5000  cc  stainless  steel  beaker.  A  precalculated 
amount  of  diazotised  suLfanilic  acid  paste  was  weighed  in  <i 
small  beaker  and  was  completely  transferred  to  the  stainless 
steel  beaker  by  rinsing.  A  glass  rod  was  used  as  a  stirrer 
to  facilitate  rapid  dissolution.  I’sualiv  it  took  10  minutes. 


That  was  the  feed  solution  B  which  was  transferred  to  the 
2000  cc  beaker  connected  to  the  injector.  The  stage  was  now 
set  for  an  experiment. 

(11)  One  person  was  required  to  position  the  fiber-optic  probe 

in  the  reactor.  The  same  person  also  monitored  the  tempera¬ 
tures,  maintained  the  flow  rate  of  the  mainstream  at  a  constant 
level,  added  the  feed  solution  B  to  the  feed  beaker  to  maintain 
a  continuous  supply  to  the  injector  and  collected  effluent 
samples  to  check  for  the  stability  of  the  operation. 

The  other  person  collected  and  recorded  the  experimental 
data.  That  involved  manipulation  of  several  keys  on  the 
microprocessor  front  panel  to  begin  and  end  data  collection 
for  a  certain  probe  location.  The  person  had  to  monitor  the 
integrated  average  of  the  raw  signal  to  get  meaningful  values 
of  VA. 

(12)  An  experiment  began  when  the  overflow  chamber  outlet  was 
connected  to  the  drain  and  the  mainstream  was  turned  on.  The 
flow  rate  was  set  at  a  desired  level  bv  adjusting  the  needle 
valve . 

(13)  The  secondary  stream  flow  rate  was  set  very  high  in  the 
beginning  to  purge  any  residual  air  from  the  feedline  and 
"hairbrush"  injector,  and  then  the  rate  was  uuiusted  to  the 
proper  level  by  turning  the  micrometer  screw. 

( 13 )  A  small  beaker  of  pure  A  solution  was  kept  at  hand  for  use 

as  the  reference  solution.  The  fiber-optic  probe  was  dipped 
into  it  to  measure  the  reference  voltage,  V'  at  each 


wavelength  A. 


(15)  The  probe  was  then  lowered  into  the  reactor  and  the  data 
collection  begun.  I-'hen  the  person  monitoring  the  micro¬ 
processor  indicated  that  the  signal  was  statistically  stable, 
the  probe  was  moved  to  the  next  location.  The  sampling  time 
for  the  microprocessor  was  50  psec.  Two  thousand  forty-eight 
points  per  record  were  used  and  typicallv  32  such  consecutive 
records  gave  a  stable  value  that  corresponded  to  about  three 
seconds  of  real  time.  For  each  X,  six  different  axial  stations 
were  used.  After  measuring  at  the  sixth  location,  the  probe 
was  again  dipped  into  the  reference  solution  to  check  for  any 

baseline  drift.  An  arithmetic  mean  of  the  two  -  values 

ref 

was  used  in  calculation  of  AV^ .  Thus,  for  each  X  a  total  of 
eight  voltages  were  measured.  The  amount  of  feed  A  solution 
available  permitted  such  measurements  at  five  different  wave¬ 
lengths. 

(16)  At  Che  end  of  the  run,  the  metering  pump  was  turned  off  first 
and  tap  water  was  added  to  the  0.21  m  lank  to  continually 
supply  fluid  to  the  reactor.  That  permitted  washing  out  of 
the  colored  products. 

(17)  The  feed-beaker  of  stream  B  was  washed  and  filled  with  tap 
water  and  the  injector  was  purged  with  a  high  flow  rate.  That 
removed  all  the  residual  B  from  the  injector  and  prevented  its 
reaction  with  the  injector  wall  bv  diffusion  through  the 
protective  coating. 

(18)  Three  effluent  samples  were  collected  during  the  run  and  were 
analyzed  at  the  end  of  the  experiment.  The  procedure  for 


analvsis  ot  effluent  samples  was  the  same  as  before.  The 
probe  was  first  dipped  into  the  sample  and  then  into  the 
reference  solution;  the  voltage  difference  provided  the  needed 
aV'‘  value. 


6 .  Experimental  Conditions 

Thirty-seven  mixing-reaction  experiments  were  carried  out  over  a 
period  of  three  months.  The  selected  experimental  conditions  for 
successful  experiments  are  summarized  in  Table  8.  In  one  set  of  runs, 
the  effect  of  feed  stoichiometry  S  (ratio  of  B  and  A)  on  the  product 
distribution  was  studied  at  a  fixed  mixing  intensity.  Three  different 
stoichiometries  (3  =  0.39,  1.00,  and  1.54)  were  used.  In  another  set 
of  runs,  at  a  fixed  £  of  1.51,  the  mixing  intensity  was  varied  by 
altering  the  mesh  size  and  mesh  placement. 

In  general,  all  the  experimental  conditions  were  repeated  at  least 
three  times  to  check  for  the  reproducibility.  In  some  cases,  more  than 
three  runs  were  made  to  obtain  data  from  the  region  very  close  to  the 
injector  to  the  farthest  end  of  the  reactor.  A  six-inch  section  of 
reactor  was  added  to  or  removed  from  the  reactor  length  for  that  purpose 
It  was  connected  to  the  rest  of  the  reactor  by  a  short  two-inch  section 
of  hose  pipe  and  clamps. 

The  raw  data  thus  collected  are  compiled  in  Appendix  G.  The  analy¬ 
sis  of  these  data  and  the  consequent  discussion  is  considered  in  the 
next  chapter. 
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VI.  RESULTS  AND  DISCUSSION 


A.  Fluid  Mechanics 

The  turbulent  reactive  flow  considered  in  this  study  is  a  special 
case  where  reaction  takes  place  in  a  dilute  liquid  medium.  The  dilute 
nature  of  the  reactive  flow  field  provides  a  clear  separation  between 
fluid  mechanics  and  chemistry.  The  flow  is  isothermal  and  the  viscosity 
of  the  fluid  also  remains  unchanged  during  the  reaction.  As  a  result, 
the  velocity  field  drives  the  concentration  field  but  the  latter  does 
not  affect  the  former.  The  measurement  of  velocity  is  thus  useful  not 
only  to  characterize  the  flow  field  but  also  to  infer  something  about 
the  micromixing. 

The  turbulence  measurements  were  conducted  at  a  pipe  Reynolds 
number  of  3400,  which  was  the  same  as  in  the  mixing-reaction  experi¬ 
ments.  The  measured  axial  mean  velocity  at  the  reactor  centerline 
is  plotted  as  a  function  of  axial  distance  in  Figure  20(A),  and  the 
corresponding  turbulence  intensity  in  Figure  20(B) .  The  zero  on  the 
abscissa  of  these  figures  corresponds  to  the  plane  of  the  injector 
needles . 

The  mean  velocities  plotted  in  Figure  20(A)  incorporate  the 
correction  due  to  the  rise  velocity  of  bubbles.  Generally,  this 
correction  is  less  than  3  percent.  The  rise  velocity  of  hubbies 
corresponds  to  a  mean  bubble  diameter  (based  on  Stoke' s  law)  of  about 
60  um,  which  is  well  within  the  range  of  particle  sizes  generally  used 
in  PUDV  measurements.  The  axia.1  mean  centerline  velocity  profile 
shows  how  dramatically  the  flow  speeds  up  during  its  passage  through 
the  contraction  nozzle.  The  axial  profile  becomes  relatively  flat  as 


Che  flow  enters  the  reactor  pipe  although  the  velocity  continues  to 
increase  gradually  over  the  reactor  length. 

The  mean  centerline  velocity  increases  sevenfold  in  going  from  the 
settling  chamber  to  the  reactor  entrance  while  the  area  of  cross-section 
for  the  flow  decreases  by  a  factor  of  nine.  Material  balance  requires 
that  the  cross-sectional  average  velocity  also  increase  ninefold. 
However,  the  measured  velocities  are  centerline  velocities  and  are 
maximum  velocities  at  a  given  axial  location.  Therefore,  a  ratio  of 
seven  implies  that  the  radial  profiles  of  mean  axial  velocity  at  the 
two  locations — settling  chamber  and  reactor  entrance — are  not  completely 
flat.  This  may  have  resulted  from  growth  of  the  wall  boundary  layer. 
This  growth  is  further  reflected  in  the  gradual  increase  of  the  mean 
axial  velocity  over  the  reactor  length.  The  growing  boundary  layer 
acts  as  an  invisible  contraction  section  and  accelerates  the  flow  in 
the  pipe  core.  Figure  20(A)  reveals  a  16  percent  increase  in  the 
velocity  from  entrance  of  the  reactor  to  the  exit.  This  applies  to  all 
the  cases  studied.  Consequently,  it  is  important  to  determine  by  how 
much  the  flow  deviates  from  its  desired  one-dimensional  character. 

Velocity  measurements  were  made  at  the  reactor  exit  at  several 
radial  locations  (see  Appendix  E) .  Owing  to  the  limitations  of  the 
probe  design,  it  was  possible  to  access  only  the  central  3  cm  core  of 
the  reactor.  The  radial  profile  of  the  mean  axial  velocity  was  indeed 
found  to  be  relatively  flat.  The  maximum  variation  of  the  velocity  was 
less  than  3  percent  from  the  mean  axial  centerline  velocity.  These 
results  suggest  that  the  mean  velocity  field  was  uniform  in  the  central 
core  over  the  entire  length  of  the  reactor.  Since  the  boundary  layer 
growth,  which  causes  the  acceleration  of  the  core,  is  confined  to  the 
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unaccessible  annulus  between  the  core  and  the  pipe  wall,  the  boundary 


layer  thickness  can  only  be  estimated.  Assuming  a  flat  mean  velocity 


profile  in  the  central  core  and  a  linear  mean  velocity  profile  in  the 


boundary  layer  with  no  slip  at  the  wall,  it  is  possible  to  calculate 


the  boundary  layer  thickness  which  will  cause  a  16  percent  rise  in  the 


core  velocity  at  the  exit.  The  calculation  involves  a  simple  mass 


balance  of  the  fluid  between  the  reactor  entrance  and  the  exit.  Such  a 


calculation  yields  a  boundary  layer  thickness  of  0.37  cm.  The  actual 


boundary  layer  would  be  thinner  than  this.  Since  the  boundary  layer  is 


thin  and  the  measured  radial  profile  of  the  mean  axial  velocity  is 


flat,  it  is  concluded  that  the  reactor  provides  a  nearly  uniform  mean 


axial  velocity  field  along  its  length. 


The  axial  turbulence  intensity  profile  depicted  in  Figure  20(B) 


can  also  be  used  to  make  inferences  about  the  homogeneity  of  the 


turbulent  flow  field.  Note  that  the  plotted  turbulence  intensities 


incorporate  corrections  for  the  electrical  noise  as  well  as  the  random 


noise  associated  with  bubbles.  The  corrections,  which  involve  simply 


subtracting  the  noise  levels  from  the  measured  signals,  imply  that  the 


noise  processes  are  uncorrelated  with  the  true  turbulence.  No  attempt 


is  made  here  to  justify  this  assumption  since  the  turbulence  intensity 


values  are  not  used  in  any  critical  computations,  but  mainly  as  an 


indication  of  the  homogeneity  of  the  turbulence.  The  measured  and 


corrected  values  may  serve  as  upper  and  lower  bounds  on  the  actual 


turbulence  intensitv.  It  should  be  noted  that  no  such  corrections  were 


applied  to  the  fully  developed  pipe  flow  data  (cf.  Section  V.A.2), 


since  the  bubble  sizes  were  not  measured  in  those  experiments.  Con¬ 


sidering  the  fact  that  the  signal  level  was  three  to  six  times  higher 
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in  those  experiments  compared  to  the  reactor  measurements,  the  effect 
of  the  noise  is  expected  to  be  reasonably  small  and  should  not  affect 
the  conclusions  of  Section  V.A.2. 

It  is  evident  that  the  turbulence  intensity  decavs  very  rapidly 
to  an  asymptotic  value  of  about  4  percent  (10  percent  if  uncorrected) 
within  5  cm  downstream  of  the  injector.  Nearly  the  same  asymptotic  level  is 
reached  in  all  the  cases  studied  and  this  includes  when  the  screen 
is  placed  upstream  or  downstream  of  the  injector.  It  is  interesting  to 
note  that  when  the  screen  is  placed  upstream  of  the  injector,  it  does 
not  seem  to  alter  the  intensity  profile  at  all  as  compared  to  the  case 
where  only  the  injector  (compare  circles  and  squares  in  Figure  20(A)). 
In  terms  of  micromixing,  this  indicates  that  placing  a  screen  in  the 
upstream  location  may  not  induce  any  additional  mixing  compared  to  that 
induced  by  the  injector  itself.  The  highest  turbulence  level  is 
measured  within  2  cm  downstream  of  the  screen  when  it  is  in  its 
downstream  location  (note  the  elevated  positions  of  triangles  in  Figure 
20(A)).  This  implies  that  this  configuration  may  lead  to  better 
micromixing  than  the  other  cases. 

In  any  case,  the  first  one-eighth  of  the  reactor  has  varying 
turbulence  level  while  the  last  seven-eighths  has  an  almost  constant 
turbulence  intensity.  Actually,  the  turbulence  intensity  may  well  be 
decaying  over  the  entire  length  of  the  reactor,  but  with  the  present 
measurement  technique,  the  weak  decay  is  not  measurable.  An  estimate 
of  the  rate  of  decay  of  velocity  fluctuations  which  supports  this  idea 
is  given  below. 

Most  experimental  data  on  flow  behind  grids  are  correlated  by  the 
following  equation  (Hinze,  1959) : 
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rms 


-)  =  J5  [x  -  x  ] 

M  o 


(21) 


In  this  equation,  U  and  u  are  the  mean  and  rms  velocities,  respec- 

rms 


tively;  m  is  a  constant  which  depends  on  mesh  size,  shape  and  geometry; 
x  is  the  axial  distance  downstream  of  grid,  M  is  the  mesh  length  and 


x^  is  the  virtual  origin  where  experimental  data  intersect  the  x  axis 


,2  . 


when  ( u/ )  is  extrapolated  to  zero.  For  the  square-mesh  grids  used 


in  the  present  study,  m  is  estimated  to  be  86  based  on  many  measurements 
for  square-mesh  grids  with  round  bars  (Hinze,  1959).  The  axial  distance 
at  which  the  asymptotic  turbulence  level  is  reached  (5  cm)  is  taken  as 


the  virtual  origin  (x  )  and  M  =  0.42  cm  (M, ) .  It  is  estimated  that  the 

o  1 


turbulence  intensity  would  vary  from  3.1  percent  at  x  =  10  cm  to  1.2 
percent  at  x  =  40  cm,  the  exit  of  the  reactor.  Clearly,  the  screen 
generates  a  relatively  low  level  of  turbulence  which  decays  very  slowly 
over  the  reactor  length. 

The  radial  profile  of  turbulence  intensity  at  the  reactor  exit 
reveals  a  reasonable  uniformity  (within  10  percent  at  0.8  cm/sec)  over 
the  cross-section. 

The  foregoing  discussion  indicates  that  the  reactor  can  be 
divided  into  two  zones.  The  first  zone  corresponds  to  the  initial  5  cm 
of  the  reactor  which  has  higher  but  varying  turbulence  level  and  hence 
better  but  changing  micromixing  conditions.  The  second  zone  corre¬ 
sponds  to  the  rest  of  the  reactor  which  is  reasonably  homogeneous  for 
both  the  mean  and  fluctuating  velocity  fields.  This  zone  has  uniform 
micromixing  conditions  and  is  better  suited  for  characterization  by  .i 


single  value  of  the  micromixing  parameter  x  . 


m 


m 
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The  first  zone  covers 


only  one-eighth  of  the  reactor  and  is  not  expected  to  influence  the 
reactions  much  except  when  the  screens  are  placed  downstream  of  the 
inj  ector . 

Finally,  it  is  useful  to  employ  turbulence  intensity  values  in 

estimating  the  micromixing  time  (t  ) .  T  is  given  by  Equation  (2) 

m  m 

where  e  may  be  estimated  as 


£  ~  A  u^  /£ 
rms  e 


(22) 


The  proportionality  constant  A  is  cr  order  unity  (Hinze,  1959).  in 
the  above  equation  is  the  characteristic  integral  length  scale  of  the 
velocity  field  and  corresponds  to  a  characteristic  transverse  dimen¬ 
sion  of  the  flow  field,  that  of  the  dominant  turbulence  generator. 

There  are  three  turbulence  generators  in  the  present  system:  honeycomb 
injector  and  screen.  Fortunately,  all  of  them  have  similar  character¬ 
istic  dimensions:  straw  diameter  =0.6  cm,  injector  tube  spacing  = 

0.65  cm  and  the  mesh  spacing  (M^)  =  0.42  cm.  Any  one  of  the  above 

dimensions  is  adequate  for  the  purpose  of  estimating  t  .  Z  =0.65  is 

m  e 

chosen  in  the  following  estimate. 

A  reasonable  estimate  of  L  ,  the  integral  length  scale  of  the 
concentration  field,  is  0.8  cm,  which  is  the  spacing  of  the  injector 
needles.  It  is  based  on  the  fact  that  the  secondary  stream  jets,  which 
carry  the  reactant  B,  expand  and  fill  the  reactor  cross-section  at 
about  2  cm  downstream  from  the  needles.  Therefore,  the  expanded  jet 
diameter,  which  equals  the  spacing  of  the  injector  needles,  is  an 
appropriate  characteristic  dimension  of  the  concentration  field. 


Similar  arguments  are  used  in  literature  to  estimate  Lg  (McKelvev  et  al 
1975).  Equation  (2)  then  becomes 


90 


L~Z  1/3  v£  1/2 

x  *  0. 5 [4(— f— +  (— --)  In  NV  ]  (23) 

m  J  J  ■' 

u  u  L 

rins  rms 

2 

For  the  above  estimates  of  Z  and  L  and  V  =  0.008  cm  / s  (kinematic 

e  s 

viscosity  of  water  at  25°C)  and  a  typical  liquid  phase  Schmidt  number 

(Nr  )  of  700,  two  estimates  of  T  are  obtained: 

S  c  Ifl 

T  ^  2.20  sec  (u  =  0.8  cm/s  (uncorrected)) 
m  rms 

and 

T  ~  8.63  sec  (u  =  0.23  cm/s  (corrected)), 
m  rms 

Thus,  fluid  mechanical  measurements  indicate  that  the  micromixing  time 
in  the  reactor  may  lie  between  two  and  nine  seconds.  This  magnitude 
is  quite  striking  because  if  indeed  the  mixing  time  is  in  this  range, 
the  reactor  would  be  the  most  segregated  reactor  ever  to  be  used  in 
any  mixing-reaction  experiment.  (For  comparison,  the  reactor  of 
Vassilatos  and  Toor  (1965)  had  a  mixing  time  of  0.0082  sec.) 

B.  Mixing  and  Reaction 


1 .  Preliminaries 

Before  the  results  of  the  mixing-reaction  experiments  are  dis¬ 
cussed,  several  preliminary  points  are  considered.  An  important 
element  of  any  experimental  mixing-reaction  study,  such  as  the  present 
one,  is  to  have  reliable  kinetic  information.  This  refers  to  reaction 
rate  expressions  and  the  associated  rate  constants.  Generally,  a 
stopped-flow  apparatus  is  required  to  study  the  kinetics  of  such  fast 
reactions  as  employed  in  the  present  study.  Unavailability  ot  sucu  an 


apparatus  has  compelled  us  to  use  the  kinetic  information  reported  in 
the  literature  (Kozicki,  1980;  Bourne  et  al.,  1981).  Both  the 
reactions  (6)  and  (7)  are  reportedlv  irreversible  and  second  order  in 
concentration  dependence: 


"rA  =  klCACB 


rS  k2CRCB 


Bourne  et  al.  (1981)  report  the  following  values  for  the  rate  con¬ 


stants  at  pH  =  10.00  and  298  K: 


k^  =  7.3  x  10"*  m^/mol/s 


kn  =  3.5m  /'col/ s 


Hie  temperature  dependence  of  k0  is  reported  bv  Kozicki  (1980) : 


=  1.93  x  10'  exp(-5989/T)  m"*/mol/s 


It  is  instructive  to  see  how  the  second  azo-coupling  reaction 
(Equation  (7))  was  studied  by  these  workers.  Equimolar  solutions  of 
R  and  B  were  injected  into  the  stopped-flow  apparatus  and  absorbance 
was  measured  at  \  =  620  nm  as  a  function  of  time.  Absorption  due  to  R 
was  neglected  at  that  wave  Length  and  the  absorbance- L ime  data  was 
used  to  calculate  k7  (Kozicki,  1980). 

This  method  has  some  obvious  limitations.  Estimated  values  or 

9 

C  and  c  at  620  nm  are  126  and  850  mVmol,  respectively.  ClearLv. 
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absorption  due  to  R  mav  be  small  but  not  negligible.  Since  absorbance 
was  attributed  entirely  to  S,  the  reaction  rate  was  overestimated  by 
some  unknown  amount,  resulting  in  a  higher  than  actual  value  of  k£. 
Recently,  Baldyga  and  Bourne  (1984)  reported  another  expression  for  k2: 

k2  =  1.11  x  10^  exp(-3957/T)  m^/mol/s  (27) 

Baldyga  and  Bourne  (1984)  gave  no  details  of  their  experimental 
method,  but  since  the  new  value  of  k^  is  half  of  the  old  value  of  k2 
at  298  K,  Equation  (27)  appears  to  be  an  improvement  over  Equation 

(26) .  Nevertheless,  there  is  some  ambiguity  associated  with  the  value 
of  k2* 

The  question  of  the  accuracy  of  k2  acquires  further  significance 
in  light  of  variations  in  values  of  £R  and  Eg  as  measured  in  the  present 
work  and  reported  by  others  (see  Figures  16  and  17).  Figure  17 
suggests  that  if  indeed  Kozicki's  values  of  Eg  were  higher  than  the 
actual,  his  measured  Cg  values  were  lower  than  the  actual  Cg  values. 
Consequently,  k2  given  by  Equation  (26)  may  be  lower  than  the  actual 
k2*  On  the  other  hand,  errors  introduced  by  a  higher  Eg  value  and 
neglecting  the  absorption  due  to  R  may  cancel  each  other  out.  In  that 
event,  k2  given  by  Equation  (26)  may  be  more  reliable  than  Equation 

(27)  even  though  the  latter  is  the  most  recent  value.  In  the  present 

work,  the  values  of  e  and  e  used  in  the  data  analysis  are  different 

K  b 

(and  lower)  than  those  of  Bourne  et  al .  (1981).  Therefore,  k2  given  by 
Equation  (26)  is  considered  more  accurate  and  is  emphasized  in  fitting 
models  with  the  data.  However,  the  sensitivity  of  predictions  to 
variations  in  k2  are  considered. 
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The  first  rate  constant,  k^,  is  three  orders  of  magnitude  higher 

than  k^.  Therefore,  the  reported  value  of  k^  at  298  K  can  be  used 

equally  well  for  all  experimental  runs.  The  model  predictions  are 

insensitive  to  minor  variation  in  k^. 

Next,  it  is  important  to  check  the  internal  consistency  of  the 

experimental  concentration  measurements.  The  values  of  and  £^ 

R  b 

displayed  in  Figures  16  and  17  were  obtained  with  the  Hewlett-Packard 

UV-visible  spectrophotometer  and  not  the  fiber-optic  probe  used  in 

the  reactor  experiments.  In  these  experiments,  k^  and  kg  were 

obtained  through  calibration  of  the  fiber-optic  probe  prior  to  each 

set  of  experiments.  Recalling  the  fact  that  k^  =  k  and  k^  = 

RPR  b 

kpE^Jl,  it  follows  that  a  ratio  of  k^  (or  kg),  measured  by  the  probe, 
and  e  (or  £^) ,  measured  in  a  separate  spectrophotometer,  must  be  a 
constant  (k^£) ,  independent  of  compounds  as  well  as  the  wavelength. 

It  is  found  that  the  ratio  is  in  the  range  0.50-0.53  for  R  and  0.63- 
0.71  for  S  for  all  five  wavelengths.  The  relative  constancy  of  the 
ratio  for  R  or  S  is  encouraging.  The  difference  in  the  values  for  R 
and  S  may  stem  from  the  fact  that  the  optical  fiber  was  directly  in 
contact  with  the  solutions.  The  chemical  contact  may  have  affected 
(e.g.,  etched)  the  exposed  fiber  surface.  Nevertheless,  the  existence 
of  a  nearly  constant  ratio  lends  further  support  to  the  concentration 
measurement  technique. 

The  final  preliminary  point  concerns  models  to  be  used  for 
comparison  with  experimental  results.  In  Chapter  III,  where  simulation 
of  the  present  competing  reactions  was  considered,  it  was  shown  that 
of  the  five  mechanistic  models  with  a  turbulence  analogy,  the  3E  model 
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is  grossly  inadequate  for  selectivity  estimation.  The  failure  of  3E 
was  attributed  to  the  lack  of  "rich"  and  "lean"  reacting  regions  in 
the  model  structure.  The  evaluation  of  the  remaining  four  models  can 


be  attempted  by  comparison  of  model  predictions  with  the  experimental 
data.  Unfortunately,  for  the  present  experimental  conditions — which 
involve  highly  unequal  reactant  flow  rates — the  SD  model  does  not 
possess  a  turbulence  analogy,  and  the  CRD  model  required  a  staggering 
amount  of  computer  time  to  achieve  meaningful  results.  Owing  to  the 
high  inequality  of  flow  rates,  the  concentration  profiles  predicted 
by  the  CRD  model  were  very  sensitive  to  the  number  of  eddies  chosen 
for  the  simulation,  although  the  exit  concentrations  were  fairly  stable. 
For  meaningful  predictions  of  the  profiles,  either  the  number  of  eddies 
had  to  be  increased  by  an  order  of  magnitude  or  the  simulation  had  to 
be  conducted  several  times,  each  with  a  different  sequence  of  random 
numbers,  and  the  results  ensemble  averaged.  Either  of  these  alterna¬ 
tives  was  economically  prohibitive.  However,  the  previous  simulations 
(Chapter  III)  did  show  that  the  CRD  and  4E  models  give  similar  predic¬ 
tions  and  that  the  I EM  and  SD  models  also  give  similar  predictions. 

In  the  following,  therefore,  the  4E  and  IEM  models  are  used  to  compare 
model  predictions  with  the  experimental  data. 

2.  Experimental  Results 

The  results  of  the  mixing-reaction  experiments  are  displayed  in 
Figures  21-25.  The  ordinates  of  all  the  figures  are  concentrations 
normalized  by  the  (mixed)  feed  concentration  of  B,  C^q.  C^q  is 
calculated  by  CDr,  =  C„„Q  / (Q  +  Q_).  The  ordinates  of  profiles  for  S 
S  are  magnified  compared  to  R  except  in  Figure  21.  The  abscissa 
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Figure  21. 
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AXIAL  DISTANCE  ,  x  ,  cm 

Mean  concentration  profile  of  products  (6  =  0.4;  screen 
(Mi)  upstream  of  injector;  Tm  =  2.0  sec)  [4E,  4E':  4e 
model  predictions  for  k.2  given  by  Equations  (26)  and 
(27),  respectively;  MM:  maximum  mixedness  (PFR) 
predictions ] 


AXIAL  OISTANCE  ,  x  ,  cm 

Figure  22.  Mean  concentration  profile  of  products  (S  =  1.0; 

screen  (M^)  upstream  of  injector;  Tm  =  2.0  sec): 
(A)  profile  of  monoazo  dye,  and  (B)  profile  of 
diazo  dye 


AXIAL  DISTANCE,  x,  cm 


Mean  concentration  profile  of  products  (8  =  1.5, 
screens  (Mi,M2,M3)  upstream  of  injector;  xm  =  2.0 
(A)  profile  of  monoazo  dye,  and  (B)  profile  of  dia 
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is  the  axial  distance  from  the  injector  needle  tips.  The  squares 
represent  the  arithmetic  mean  of  experimental  data  at  each  location, 
and  the  "error  bars"  represent  the  range  over  which  the  experimental 
results  varied  in  successive  runs  (typically  three  runs).  The  "error 
bars"  are  not  based  on  statistical  measures  (e.g.,  standard  deviation) 
since  experiments  were  repeated  only  three  times.  The  smooth  curves 
are  the  predictions  of  various  models — maximum  mixedness  (classical 
PFR) ,  IEM  and  4E.  The  curves  labeled  as  4E'  are  4E  predictions  for  k2 
given  by  Equation  (27). 

Figures  21-25  reveal  that  the  "error  bars"  are  quite  large, 
implying  relatively  large  variations  in  the  measured  concentrations. 

The  "error  bars"  merit  some  discussion  before  the  model  discrimination 
is  pursued  in  detail.  The  magnitude  of  variations  is  similar  for  both 
R  and  S.  Many  sources  may  contribute  to  these  variations,  such  as: 
errors  involved  in  measuring  and  maintaining  the  feed  flowrates  (Q^ 
and  Q  )  at  a  constant  level;  errors  associated  with  preparing  feed 

D 

solutions  of  known  concentrations  (C  C__) ;  minor  variations  in  pH 

Ad  dU 

from  run  to  run;  errors  caused  by  a  fairly  low  signal  to  noise  ratio 
associated  with  the  present  measurement  system;  the  fact  that  the 
optical  fiber  was  exposed  directly  to  the  solutions  and  hence  subjected 
to  possible  undesirable  chemical  interactions;  and  that  the  measure¬ 
ments  could  be  made  only  at  five  wavelengths.  Apart  from  Qg,  which 

3 

is  highly  concentrated  (Cg^  ~  10-40  mol/m  ),  the  feed  conditions  were 
controlled  quite  precisely.  A  maximum  variation  in  Qg,  determined 
from  the  metering  pump  calibrations,  is  estimated  to  be 
±  0.06  ml/sec  and  can  cause  only  a  3  percent  variation  in  CgQ.  Azo 
coupling  is  an  electrophilic  aromatic  substitution  and  the 
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reactive  species  are  1-naphtholate  and  the  diazonium  ions,  not 
the  molecular  1-naphthol  and  diazotised  sulfanilic  acid.  The  concen¬ 
trations  of  the  reactive  ions  are  pH  dependent  and  so  are  the  rate 
constants  k^  and  k9  (which  are  based  on  total — ionized  and  nonionized — 

concentrations).  Nevertheless,  the  reaction  medium  was  adequately 

13  3 

buffered — 10  mol/m  and  10  mol /m  NaHCO^  (I  =  40  g  ion/m  )  as 

3 

compared  to  0.3  mol/m  1-naphthol — and  the  pH  was  found  to  be  constant 

(within  ±  0.02  of  10.00).  The  effect  of  operating  temperature  also 

is  expected  to  be  minimal.  Thus,  it  is  concluded  that  the  large 

"error  bars"  result  principally  from  inherent  limitations  of  the 

present  absorbance  measurement  system.  Of  course,  it  should  also  be 

emphasized  that  the  experimental  concentration  levels  are  very  low  and 

are  intrinsically  difficult  to  measure  accurately. 

The  input  parameters  for  the  models  were  the  feed  flow  rates 

(Qa  and  Qg) ,  the  unmixed  feed  concentrations  (Caq  and  Cbq) ,  the 

operating  temperature  (Ta  )  ,  the  mean  axial  velocity  (0)  and  the 

micromixing  time  (T  ) .  Since  each  experiment  was  repeated  at  least 

m 

three  times,  the  corresponding  values  of  the  input  parameters  were 
taken  as  arithmetic  averages  of  the  actual  input  conditions  for  all 
runs  at  the  same  nominal  operating  conditions.  A  constant  value  of 
U  =  7.6  cm/s  was  used  as  input  for  all  the  simulations.  Although  the 
value  of  k2  given  by  Equation  (26)  is  emphasized  in  the  discussion  of 
model  comparison  with  experiment,  predictions  corresponding  to  a 
given  by  Equation  (27)  are  also  plotted  (curves  4E')  to  show  the 
sensitivity  of  the  comparison  to  the  value  of  k2* 

The  only  adjustable  parameter  of  the  models  is  t  .  The  actual 
model  parameters — R<,  for  4E  and  h  for  IEM — are  related  to  Tm  through 
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the  turbulence  analogy  (cf .  Table  1).  The  models  may  be  used  to  fit 
the  experimental  data  by  varying  T  .  The  fluid  mechanical  results  in 
conjunction  with  the  turbulence  analogy  have  suggested  that  may  lie 
between  two  and  nine  seconds.  In  view  of  the  scatter  in  the  experi¬ 
mental  data,  an  "eyeball"  fit  was  used  to  determine  the  "best"  rather 
than  a  more  ■  phisticated  statistical  method. 

Ordinary  y,  both  models  should  be  used  to  fit  the  data  inde¬ 
pendently  anc  then  evaluated  for  their  relative  success.  In  the 
present  study,  however,  the  IEM  model  was  unable  to  simultaneously 
match  the  experimental  profiles  of  both  R  and  S,  regardless  of  the 
value  of  r  or  used, for  6  _>  1.  A  reasonable  match  of  the  IEM  model 
predictions  for  CR  always  led  to  very  poor  predictions  of  Cg.  In  fact, 
Cs  profiles  could  not  be  matched  at  all  unless  an  arbitrarily  higher 
value  of  k2  was  used.  For  6  =  0.4,  when  virtually  no  S  was  produced, 
reasonably  good  predictions  were  possible  with  the  IEM  model.  As  a 
consequence  of  the  overall  inability  of  the  IEM  model  to  make  reason¬ 
able  predictions,  the  question  of  model  evaluation  was  narrowed  down  to 
considering  the  ability  of  only  the  4E  model  to  predict  the  experimental 
data.  In  Figures  21-25,  the  predictions  of  the  IEM  and  maximum 
mixedness  models  are  displayed  only  to  highlight  the  differences 
between  them  and  the  4E  model  predictions.  In  addition,  the  4E  model 
predictions  for  k2  given  by  Equation  (27)  are  also  displayed  to  show 
that  conclusions  remain  unchanged  for  either  of  the  k2  values. 

Figure  21  shows  the  results  of  experiments  number  35-37.  These 
experiments  were  conducted  with  screen  M^  in  its  upstream  position. 

This  case  was  studied  earlier  from  a  fluid  mechanical  viewpoint  (Case  4) 
and  was  found  to  be  most  suited  for  model  evaluation  because  the 
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turbulence  field  was  most  homogeneous.  In  the  other  case  where  screen 

was  placed  downstream  from  the  injector,  local  inhomogeneity  in  the 

turbulence  field  very  close  to  the  screen  might  obscure  the  model 

comparison  somewhat.  The  relative  uniformity  of  turbulence  in  the 

former  case  is  preferable  for  an  unambiguous  evaluation  of  the  4E 

model.  In  experiments  35-37  the  stoichiometry,  0,  was  very  low  (0.4). 

Because  of  the  limited  B  present  and  the  fact  that  the  first  coupling 

(A  +  B  +  R)  is  three  orders  of  magnitude  faster  than  the  second 

coupling  (R  +  B  -»•  S) ,  the  concentration  of  S  is  expected  to  be  nearly 

zero.  This  is  essentially  a  single  reaction  case  and  we  do  not 

expect  to  be  able  to  discriminate  between  models  here  (recall  Chapter 

III) .  Because  the  S  concentration  was  expected  to  be  near  zero  in 

this  case,  the  experimentally  observed  absorbance  was  attributed 

entirely  to  R  and  the  approximate  method  of  concentration  calculation 

was  used  (cf .  Chapter  V.B.3).  The  reaction  results  displayed  in 

Figure  21  indicate  that  T  =2  sec  is  a  good  estimate  of  the  micro- 

m 

mixing  time.  This  value  of  t  is  in  accord  with  its  estimate  based  on 

m 

fluid  mechanical  results  (x  ~  2-9  sec).  With  T  =2  sec,  the  I EM 

m  m 

model  also  provides  a  reasonably  good  fit  to  the  data.  In  fact,  the 

fit  of  the  IEM  model  could  be  improved  if  a  somewhat  lower  value  of 

x  were  chosen.  As  mentioned  above,  this  data  set  can  not  be  used  for 
m 

model  discrimination;  however,  it  is  quite  useful  in  determining  the 

value  of  for  either.  Effectively,  this  operating  condition  (0  = 

0.4)  serves  as  a  fast  single  reaction  calibration  case  allowing  us  to 

fix  a  value  of  x  . 

m 

The  results  of  experiments  with  an  equimolar  feed  (0  =  1.0; 


experiments  8-13),  but  with  an  identical  screen- inj  ec  tor 


configuration  as  in  Figure  21  (B  =  0.4),  are  depicted  m  Figure  22. 
Although  there  is  scatter  in  the  experimental  data,  particularly  for 


S  close  to  the  injector  (due  to  very  low  values  of  Cg) ,  the  4E  model 

predictions  are  in  satisfactory  agreement  with  the  data  for  a  value 

of  t  =2  sec.  Note  that  the  hydrodynamic  conditions  in  this  case 
m 

(B  =  1.0)  are  identical  to  those  of  the  previous  case  (B  =  0.4)  and 

thus  T  should  be  the  same.  The  striking  failure  of  the  IEM  model 
m 

(x  -  2  sec)  is  also  evident  in  this  figure,  particularlv  in  the 
m 

prediction  of  S  concentration. 

The  data  for  an  even  higher  stoichiometric  ratio  (0  =  1.5),  but 
again  for  the  same  screen- inj ector  arrangement  is  shown  in  Figure  23 
(experiments  1 6- 2 3 ,  34).  Again,  the  4e  model  compares  quite 

favorably  with  the  experimental  data  for  a  value  of  =  2.  It  should 
be  pointed  out  that  experiments  22  and  23  involved  use  of  two  finer 
screens  (M2  and  Mg)  in  the  upstream  position,  and  yet  no  change  in  the 
R  and  S  concentration  profiles  were  observed  (see  Appendix  G) .  This  is 
consistent  with  our  findings  from  the  fluid  mechanical  measurements  that 
when  a  screen  was  placed  upstream  of  the  injector,  its  effect  or.  the 
turbulence  intensity  was  indistinguishable  from  that  of  the  injector 
itself.  The  IEM  model  grossly  overpredicts  the  concentration  of  R 
and  predicts  negligible  concentration  of  S,  which  is  in  stark  contrast 
with  the  experimental  measurements. 

In  view  of  the  three  cases  discussed  above,  it  is  clear  that  the 
effect  of  mixing  on  the  experimental  selectivity  (concentrations  of  R 
and  S)  is  very  large.  This  is  also  reflected  in  the  large  difference 
between  the  experimental  data  and  the  maximum  mixedness  (classical 
PFR)  predictions.  Both  models — 4E  and  IEM,  predict  a  large  effect  on 


the  concentration  of  R,  but  only  the  4E  model  is  able  to  fit  the  S 

concentration  data  well.  It  is  concluded  from  these  comparisons 

with  experiments  that  the  4E  model  is  able  to  simulate  the  effect  of 

varying  feed  stoichiometry  (8)  on  the  selectivity  (C  /(C  +2C  ))  under 

b  K  b 

constant  mixing  conditions  (constant  t  )  quite  well. 

m 

The  ability  of  the  4E  model  to  respond  to  changes  in  mixing  time 
at  a  fixed  feed  stoichiometry  (8  =  1.5)  was  tested  by  experiments  16- 
23,  24  (M.  in  upstream  position),  31-33  (M^  in  downstream  position), 
25-27  (M9  in  downstream  position)  and  28,  30  (M^  in  downstream  posi¬ 
tion).  There  was  no  observable  difference  between  the  results  for  the 
two  finer  screens  (M2  and  M^;  see  Appendix  G) .  Consequently,  those 
results  are  averaged  together  and  displayed  on  the  same  plots  (Figure 

24).  The  4E  model  predictions  for  T  =1.0  sec  are  in  reasonable 

m 

agreement  with  the  data  for  this  set  except  in  the  region  near  the 

injector  where  the  model  underpredicts  the  conversion.  For  the  case 

where  the  coarser  screen  (M^)  was  placed  in  its  downstream  position 

(Figure  25),  the  predictions  are  satisfactory  for  T  =0.5  sec.  It  is 

m 

also  clear  that  screens  in  their  downstream  location  induce  additional 
mixing  as  compared  to  upstream  position. 

It  is  important  to  recognize  the  effect  of  the  initial  zone  of  in¬ 
homogeneous  turbulence  (first  one-eighth  of  the  reactor)  on  the  experi¬ 
mental  concentrations.  The  first  reaction  (A  +  B  -►  R)  is  very  fast 
(characteristic  time  of  0.5  msec)  and  is  dominantly  mixing  limited. 

The  second  reaction  (R  +  B  -*■  S)  is  much  slower  (characteristic  time  of  1 
sec)  and  is  not  as  sensitive  to  the  reduced  mixing  time  (due  to  higher 
turbulence  intensities)  in  the  initial  zone  as  the  first  reaction.  This 
is  evident  in  the  proportionately  greater  variation  in  R  concentration 


106 


in  going  from  Figure  23  to  25  (t  =  2.0,  1.0  and  0.5  sec,  respectively) 
as  compared  to  S. 

In  view  of  the  foregoing  discussion,  it  is  concluded  that  the  4E 
model  is  a  more  reliable  model  of  mixing  and  chemical  reaction  than 
the  3E,  SD  and  IEM  models.  The  failure  of  the  IEM  model  is  simulating 
the  experimental  reactor  is  very  striking  and  may  be  attributed  to  its 
structural  deficiency  as  pointed  out  recently  by  Zwietering  (1984)  or 
earlier  in  this  thesis  (cf.  Chapter  III). 

Additional  comments  on  several  recent  suggestions  which  have 

appeared  in  the  literature  are  appropriate  in  light  of  the  results  of 

this  work.  It  has  been  proposed  that  in  a  competing  reaction  scheme 

similar  to  the  present  one,  that  the  final  concentrations  of  reaction 

products  (i.e.,  C  and  C  when  C„  =  0),  are  always  the  same  for  fixed 
Kb  d 

initial  conditions  (feed  flow  rates,  concentrations,  temperature,  pH) 
in  a  plug  flow  reactor,  regardless  of  the  axial  concentration  profiles 
(Brodkey  and  Lewalle,  1985).  This  suggestion  implies  that  any  reason¬ 
able  mixing  model  should  predict  the  same  final  selectivity,  regardless 
of  the  value  of  the  mixing  time  and  also  that  all  mixing  models  should 
yield  the  same  final  selectivity  as  long  as  the  qualifying  conditions 
of  the  above  proposal  are  met.  The  first  implication  is  in  conflict 
with  the  4E  model,  as  it  predicts  significant  differences  in  the  final 
selectivity  for  the  experimental  case  of  8  =  1.0  (case  shown  in  Figure 
22).  For  example,  when  =  2.0  sec,  CR/CBQ  =  0.57  and  Cq/CRn  =  0.21 


S  BO 


and  when  T  =0  sec,  0-/0.,,,  =  1.0  and  Cc/Cc  =  0.00.  The  second 
m  R  BO  SB 

implication  also  is  not  in  accord  with  the  present  study  as  large 

differences  are  observed  between  the  4E  and  IEM  model  predictions  at 

8  =  1.0.  For  example,  when  t  =2  sec,  at  long  times,  4E  predicts 

m 
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0.02.  However,  no  significant  differences  are  observed  between  the  two 
model  predictions  if  a  lower  value  of  t  ,  sav  the  one  appropriate  to 
Vassilatos's  reactor  (Vassilatos  and  Toor,  1965),  i.e.,  0.0082  sec,  is 
used.  In  such  a  case,  the  model  simulations  indicate  that  the  reactions 
go  to  completion  within  the  first  3  cm  of  the  reactor  length.  This 
again  demonstrates  the  importance  of  having  a  highly  segregated  experi¬ 
mental  reactor,  as  was  the  case  in  the  present  study,  to  discriminate 
among  mixing  models. 

A  surprising  feature  of  the  experimental  results  is  that  the  finer 

screens  (M»  and  M7)  seem  to  induce  less  mixing  (t  =1.0  sec)  as  cora- 
o  m 

pared  to  the  coarser  screen  (M  ;  r  =0.5  sec)  when  the  screens  are 

l  m 

placed  downstream  of  the  injector.  Recall  chat  they  do  not  affect  mixing 
when  placed  in  the  upstream  location.  Also,  both  the  fine  screens  appear 
to  induce  the  same  amount  of  micromixing.  This  is  contrary  to  one's 
intuition  that  the  finer  screens,  due  to  their  smaller  mesh  size,  would 
shed  ner  eddies  and  hence  lead  to  better  mixing  compared  to  the 

coar.  creen.  The  explanation  is  based  on  the  possibility  that  under 
the  operating  conditions,  the  finer  screens  may  not  shed  anv  eddies  at 
all.  It  should  be  noted  that  the  pipe  Reynolds  numbers  were  identical 
for  all  three  screens,  but  for  flow  behind  screens,  the  mesh  Reynolds 
number  (Re^)  is  more  appropriate  to  characterize  the  flow  regime. 

Support  for  the  above  explanation  comes  from  the  work  of  Keeler 
(1964).  Keeler's  experimental  system  was  similar  to  the  present  one. 

It  was  a  tubular  reactor  with  similar  screens  and  "hairbrush"  injector. 
Keeler  carried  out  experiments  to  determine  the  flow  rate  in  the 
reactor-screen  system  at  which  regular  shedding  of  eddies  commenced. 
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From  his  results,  one  can  readily  estimate  the  onset  of  eddy  shedding 
in  terms  of  a  critical  Reynolds  number  based  on  the  wire  diameter 
(Re^^) ,  provided  the  solidity  of  the  screen  is  known.  The  Re^  2nd 
the  actual  screen  Reynolds  numbers  (Re^)  are  listed  in  Table  9. 


Table  9 

Operating  Characteristics  of  the  Screens 


Screen 

M1 

m2 

M3 

Solidity 

0.43 

0.39 

0.60 

ReCr 

* 

50 

53 

35 

Red 

96 

43 

30 

Red  =  Ud/v. 


Clearly,  the  finer  screens  do  not  exceed  the  critical  Reynolds  number 

and  thus  may  not  shed  any  eddies  under  the  operating  conditions.  As  a 

result  they  are  less  effective  in  generating  additional  turbulence. 

Nevertheless,  they  do  cause  some  extra  mixing,  perhaps  due  to  physical 

obstruction  of  the  flow  during  its  passage  through  the  screens,  and 

hence  the  associated  mixing  time  is  smaller  (t  =1.0  sec)  than  the 

m 

mixing  time  for  the  case  where  no  screens  were  placed  downstream  of  the 

injector  (t  =  2.0  sec), 
m 

A  final  point  in  this  discussion  is  related  to  the  overall 
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exit  concentrations  of  R  and  S.  Overall  material  balance  demands  that 


So  ■  S  +  S  +  2S  <28) 

For  all  the  feed  stoichiometries  used  in  this  study  (8  <  2),  we  should 

expect  C  =  0  at  the  exit.  Thus, 

B 

So  -  S +  2S  (29) 


The  values  of  C  ,  C_  and  Cc  at  the  reactor  exit  are  listed  in  Table  10 
and  the  percentage  error  in  the  material  balance  is  shown  in  the  last 
column  of  the  table.  The  material  balance  is  satisfied  within  10  per¬ 
cent  for  6  £  1.0.  However,  for  8  ~  1*5,  with  the  exception  of  experi¬ 
ment  16,  the  output  of  B  (CR  +  2C<.)  always  exceeds  the  input  (CBQ)  by 
a  significant  amount. 


It  is  clear  from  Equation  (29)  that  any  error  in  measuring  Cg  is 

amplified  by  a  factor  of  two  in  the  material  balance  calculations.  It 

is  reasonable  to  expect  errors  in  C0  and  C  values  whose  magnitudes 

K  b 

are  similar  to  those  implied  by  the  "error  bars"  (i.e.,  variation  by 

3  -  - 

0.02  to  0.06  mol/m  ).  For  experiment  number  30,  if  C  and  Cc  values 

K  b 

3 

were  decreased  by  0.06  mol/m  ,  the  overall  material  balance  would  be 
met  within  10  percent  rather  than  50  percent  as  reported.  It  is  also 
interesting  to  note  that  experiment  number  34  has  an  80  percent  error 
in  the  mass  balance  and  experiment  number  22  has  only  a  13  percent 
error,  and  yet  both  the  experiments  yielded  very  similar  concentration 
profiles  of  R  and  S  (see  Appendix  G) .  The  main  difference  between 
these  two  experiments  is  only  in  the  exit  concentration  of  S  (0.13 
mol/m  for  Run  #22  and  0.28  mol/m  for  Run  #34) .  Thus,  our  inability 
to  measure  Cc  very  accurately,  mav  lead  to  such  discrepancies  in  the 


material  balance  as  are  seen  in  experiments  19  to  34.  It  is 

felt  that  placing  a  glass  seal  at  the  probe  tip  to  protect  the  fiber 

ends  will  improve  the  accuracy  of  the  fiber-optic  probe  considerably 


VII.  SUMMARY  AND  CONCLUSIONS 


Several  important  points  have  emerged  from  the  present  work. 

First,  it  was  shown  in  Chapter  III  that  the  five  mechanistic  models 
of  mixing  and  chemical  reaction  having  a  turbulence  analogy — namely, 
the  CRD,  SD,  IEM,  3E  and  4E  models — are  not  equivalent  in  terms  of 
estimating  selectivity  of  competing  reactions  in  a  PFR.  The  non¬ 
equivalence  of  models  was  demonstrated  via  numerical  simulations, 
although  the  CRD  and  4E  models,  the  SD  and  IEM  models,  and  the  3E  and 
maximum  mixedness  models  were  found  to  yield  similar  predictions. 

Based  on  simulations,  a  highly  segregated  plug  flow  reactor  was 
designed  and  constructed  for  an  experimental  study  of  the  effect  of 
mixing  on  selectivity  of  two  competing  azo-coupling  reactions.  The 
fluid  mechanical  measurements  discussed  in  Chapter  VI. A  provided 
sufficient  evidence  to  conclude  that  the  constructed  plug  flow  reactor 
yielded  a  nearly  homogeneous  mean  axial  velocity  field  and  a  turbu¬ 
lence  level  which  was  homogeneous  over  7/8  of  the  reactor  length.  The 
mixing-reaction  data  collected  in  this  study  were  fairly  reproducible 
in  spite  of  the  low  concentration  levels  at  which  the  reactor  was 
operated.  The  data  provided  a  valuable  test  to  discriminate  between 
two  representative  models — 4E  and  IEM.  The  4E  model  was  found  to 
successfully  simulate  the  experimental  behavior  at  all  the  conditions 
studied  while  the  IEM  model  was  found  unable  to  match  the  experimental 
data.  The  failure  of  the  IEM  and,  as  shown  earlier,  3E  models,  and 
the  relative  success  of  the  4E  model,  point  to  the  importance  of 
certain  structural  features  of  mixing  models,  such  as  the  existence 
of  "rich"  and  "lean"  reacting  regions  and  presence  of  regions  of  total 
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segregation.  The  SD  model  is  inapplicable  to  unequal  flow  rates  and 


hence  can  not  be  used  to  simulate  the  present  experiments.  The  CRD 
model  may  be  able  to  compare  favorably  with  the  data  but  requires 


prohibitively  large  computer  time. 

It  is  also  concluded  that  a  discrimination  among  several  plausible 
mechanistic  models  of  mixing  and  chemical  reaction  is  possible  only 
under  conditions  of  imperfect  mixing  similar  to  ones  covered  in  the 
present  study.  The  present  experimental  data  provide  valuable  infor¬ 
mation  to  test  many  turbulence  closure  schemes  for  concentration 
correlations  as  well.  The  plug  flow  reactor  of  the  present  study  is  a 
very  useful  arrangement  to  further  study  the  interaction  of  turbulence, 
mixing  and  chemical  reaction  in  an  unambiguous  manner. 

Finally,  the  4E  model,  which  has  many  attractive  features,  like 
its  applicability  to  arbitrary  feed  flow  rates,  arbitrary  feed 
residence  time  distributions,  and  arbitrarily  complex  kinetics;  its 
mathematical  simplicity  (O.D.E.'s)  and  modest  computational  require¬ 
ments;  the  fact  that  its  model  parameter  can  be  readily  estimated 
a  priori  or  measured  from  fast  single  reaction  experiments,  is 
recommended  as  the  most  general  mechanistic  chemical  reactor  model  for 
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homogeneous  systems  now  available. 
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VIII.  FUTURE  WORK 


Many  industrial  reactors  are  inherently  limited  by  imperfect 
micromixing.  The  flow  field  of  such  reactors  are  so  complex  that  a 
mechanistic  model  like  the  4E  model  might  be  a  powerful  tool  in  a 
systematic  study  of  the  reactor  behavior  as  compared  to  more  tradi¬ 
tional  turbulence  models.  Based  on  the  encouraging  results  of  the 
present  study,  it  seems  reasonable  to  consider  analyzing  reported  data 
from  a  high  intensity  or  low  calorific  value  fuel  combustor,  or  a 
precipitation  reactor,  or  a  polymerization  reactor  within  the  frame¬ 
work  of  the  4E  model. 

Experimentally,  the  present  reactor  is  a  very  useful  system  for 
obtaining  reliable  mixing  and  reaction  data  under  relatively  uniform 
and  controlled  fluid  mechanics.  However,  the  fluid  mechanics  of  the 
present  reactor  should  be  better  characterized  by  more  accurate  methods 
such  as  hot-film  or  laser-doppler  anemometry.  Even  the  chemistry  side 
of  the  present  work  shows  need  for  more  accurate  kinetic  information. 
Moreover,  the  present  fiber-optic  system  has  demonstrated  its  ability 
to  provide  meaningful  mean  concentration  data.  Attempts  should  be 
made  to  improve  the  signal  to  noise  ratio  of  this  system  so  that  even 
fluctuating  components  of  the  concentrations  may  be  measured.  Such  a 
data  base  would  contribute  significantly  in  evaluation  of  many  avail¬ 
able  closure  schemes  in  turbulence  modeling. 
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From  a  modeling  point  of  view,  incorporating  a  variation  in  the 
mixing  time — a  critical  factor  in  practical  situations — would  greatly 
enhance  the  generality  of  the  4E  model. 
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The  computer  codes  given  in  this  appendix  were  used  in  all  the 
model  computations  fcr  competing  reactions  in  this  study.  The  main 
input  parameter  for  the  models  are: 

3 

QA,  QB  =  Feed  flow  rate  of  two  streams  (cm  /s) 

3 

CAO,  CBO  =  Unmixed  feed  concentration  (mol/cm  ) 

U  or  VEL  =  Mean  axial  velocity  in  the  reactor  (cm/s) 

TEMP  =  Mean  reactor  temperature  (°C) 

TAUM  -  Mixing  time  (s) 

These  values  are  either  read  as  data  or  directly  specified  in  the 
main  program.  The  rate  constants,  RK1  and  RK2,  also  are  either  read  as 
data  or  internally  specified  or  calculated  within  the  program.  All 
the  governing  equations  are  in  their  dimensionless  forms.  Other 
parameters  are  specified  as  needed  by  the  respective  numerical  schemes 
for  the  numerical  integration  purposes  or  for  obtaining  printouts  at 
appropriate  time  steps.  The  printed  values  are  normalized  values  of 
concentrations  relative  to  (mixed)  concentration  of  B,  C  . 

D  U 

The  model  equations  for  the  4E  and  IEM  models  are  as  follows: 

The  4e  Model 


* 

^  <C°*  -  CL*)  exp(-R  Tt*> 

dt  s  J1  Jl 


k  .  L*  L*  ,  * 

+  XT  ( C . ,  -  C  .  . )  ]  ±  Y  .  . 

jk  ji  'ji 

i 


j  ~  A ,  B ,  R ,  S 
i  =  1,2 
k  =  3-i 


(A.  1) 


where  notations  are  in  accordance  with  Mehta  and  Tarbell  (1983a). 


is  the  local  production  (+)  or  consumption  (-)  term  given  in  terms  of 


local  concentrations  and  rate  constants  (k.  ana  k.,) . 


APPENDIX  A.  1— Slab  Diffusion  (SD)  Model 


cc  - 

CC  SLAB  DIFFUSION  MODEL  FOR  COMPETING  REACTIONS  - 

cc - 

IMPLICIT  REAL*3(A-H.O-Z) 

COMMON /CONC/  CA ( 200 ) ,C8(200) ,CR(200) ,CAOLD(200) 

COMMON/TRIO/  B(200)  ,A(200)  . C ( 200 )  ,0(200) ,V(200) 

COMMON/PREV/  CAP(200)  ,C8P(200)  ,CRP(200) 

COMMON/GRID/  X ( 200 ) , CC ( 200 . 3 ) 

PI  =  3.1415926536000 

CC - 

CC  DEFINE  SYSTEM  PARAMETERS 

CC - 

T  AUM=  0.630 
QA  =  1 .DO 
QB  =  1 .DO 
CBO  =  0.2D-3 
CAO  =  0.2D-3 
U  =6.30 
RK 1  =  7.3D6 

RK2  =  3.5D3 
PA  =  QA  /  (QA+QB) 

PB  =  QB  /  (QA+QB) 

CBOBAR  =  CBO’PB 
CAOBAR  =  CAO«PA 
BETA  =  CBOBAR/CAOBAR 
GAMMA =  RK2/RK1 
TSCALE  =  2.D0*TAUM*PI**2 
DA  =  CB0*RK1 »TSCALE 
WR I TE ( 6 , 800 )  QA.QB.PA.PB 
WRITE (6. 100)  CAO, CBO, BETA. RK1 , RK2, GAMMA 
100  F0RMAT(/2X, 'CAO  =  ' . F8 . 3 . 2X , ' CBO  =  ' , F8 . 3 . 2X , ' BETA  =  ' . F8 . 3 . 2X 

*  ,'RK1  = ' , D 1 1 . 4 , 2X , ' RK2  =  ' , D 1 1 . 4 . 2X , ' GAMMA  =  '.Dll. 4/) 

WRITE ( 6 , 1 10)  TAUM.  DA 

110  F0RMAT(2X, 'TAUM  =  ' . 0 1 1 . 4 , 2X . ' DA  NUM  =  '.D11.4) 

800  FORMAT ( 2X , ' QA  =  ' , F 1 0 . 4 . 2X , ' QB  =  ' , F 1 0 . 4 , 2X . ' PA= ' , D 1 1 . 4 , 

*  2X , '  PB= ' ,D1 1 .4) 

CC  - 

CC  DEFINE  NUMERICAL  SOLUTION  PARAMETERS 
CC - 

C  -  DT  =  TIME  STEPS  .  NT=  *  OF  TIME  STEPS  TO  BE  RUN  - 

OT  =  0.001 
NT  =  400 
NXA  =  49 
NXB  =  49 

CC  -  NXA 1 =  »  OF  GRIDS  IN  PART  A  ,NAB1=  *  OF  GRIDS  IN  PART  B  - 

NXA 1  =  NXA+1 

NXB 1  =  NXB+1 
NX  =  NXA1+NX81 
OXA  =  PA / ( NXA+O . 5 ) 

DXB  =  PB/ ( NXB+O . 5 ) 

DXA2  =  DXA  *  *  2 
ALFA  =  DXA2/DT 
DX82  =  OXB*  *  2 
ALFB  =  0XB2/DT 
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0XL2  =  (  (DXA  +  OXB)  **2/4.00  +DXA**2  )  /2.D0 
DXH2  =  (  (DXA+DXB)**2/4.D0  +DXB**2  )  /2.D0 
ALFL  =  0XL2/DT 
ALFH  -  DXH2/DT 
00  1 4  I  =  1  ,  NXA  1 

14  X(I)  =  DXA *  (  I  -  1 ) 

DO  24  I  =  1  ,  NXB 1 

24  X(I  +  NXA1)  =  X(NXA1  )  <■  (  DXA  +  DXB  )  /  2  .  DO  +  DXB  *  (  I  -  1  ) 

X( 1 )  =  0.00 
X ( NX ) =  1.00 

WRITE (6. 180)  NXA 1 , NXB 1 , DXA , OXB , OT , NT 
180  FORMAT ( 2X ,  'NXA 1  =  '  . I  5 , 2X ,  ' NXB 1  =  '  . 1 5 , 2X .  'DXA=' . D 1  1 . 4 , 2X . 
*  ‘DXB= ' , 01 1 .4, 2X, ' DT= '  ,  0 1  1  .  4 , 2X  ,  ' NT = ' ,15/) 

MX  =  150 
TOR  =  1 . 00-3 

CC  - 

CC  INITIAL  CONDITIONS 

CC - 

DO  2  1=1 .NX 

I  F  ( I  .LE.NXA1 )  CA( I )  =  1  .00 
I F ( I .LE.NXA1 )  CB(I)=O.DO 
I F ( I .GT.NXA1 )  CA ( I ) =0 . DO 
I F ( I .GT.NXA1 )  CB(I )=1 .DO 
CRCI)  =  0.00 
2  CONTINUE 

WRITEC6.145)  (CA( I ) , 1=1 .NX) 

WRITEC6.145)  (CB( I ) , 1=1 , NX ) 

145  FORMAT ( 7D 1 1 .4) 

CC  START  THE  NUMERICAL  SOLUTION  ROUTINES 

CC  SEQUENCE  :  A  — >  R  — >  B  — >  CHECK  CONVERGENCE  OF  B 

CC  I  | 

CC  l-< - 1 

TIME  =  0.00 
DO  1  LT  =  1 , NT 
DO  17  1=1. NX 
C A P (  I  )  =  CA  (  I  ) 

CBP ( I )  =  CB  (  I  ) 

CRP( I )  =  CR(  I  ) 

17  CONTINUE 
NPRINT  =  Q 

IF(LT/10*10.EQ.LT)  NPRINT=1 
TIME  =  TIME+DT 
DO  5  NS  =  1 .MX 


CC - 

CC  SOLVE  FOR  CA 
CC  - 


CC 

0.C. 

AT  X=0 . DO 

A  (  1  )  = 

-(2.D0  +  ALFA  +  DA*CB( 1  ) *0XA  2 ) 

C  f  1  )  = 

2.  DO 

0(1)  = 

-ALFA*CAP( 1 ) 

CC - LOWER  PART 


00  15  1=2. NXA 


STBVTBVWV^C’^'.V^  ’”-■' 


B( I)  =  1 . DO 

A  (  I  )  =  -(2.00+ALFA+DA*C8(I)*DXA2) 

C(I)  =  1 .DO 

0(1)  =  -ALFA«CAP(I) 


CONTINUE 


SPECIAL  CASES  NEAR  THE  INTERFACE  LINE 
B(NXA1)  =  1.00 

A ( NX A  1  )  =  - ( 2 . 00+ ALF L+DA *CB ( NXA 1  }*DXL2) 

C ( NXA 1 )  =  1.00 

0(NXA1)  =  -A LFL* CAP ( NXA 1 ) 


B(NXA1+1)  =  1.00 

A ( NXA 1  +  1  )  =  -(2.00+ALFH+DA»CB(NXA1+1 )*DXH2) 

C( NXA 1+1 )  =  1 .DO 

D(NXAl-M)  =  -ALFH*CAP(NXA 1 + 1 ) 


-  UPPER  PART 


00  25  1=2, NXB 
B ( NXA 1  + 1  )  =  1 
A ( NXA 1 + I )  =  - 
C(NXA1+I)  =  1 
0 ( NXA 1  + 1  )  =  - 


=  -(2.00+ALFB+0A*CB(NXA1+I )*DXB2) 


=  -ALFB*CAP(NXA 1 +1 ) 


CONTINUE 


B.C.  AT  X=1 .00 


B ( NX )  =  2. DO 

A ( NX )  =  -(2.D0+ALFB+DA«CB(NX)*DXB2) 
D ( NX )  =  -ALFB»CAP(NX) 

CALL  TRIDAG( 1 , NX , B , A , C , 0 , V ) 

DO  3  1=1 .NX 
C  A  (  I  )  =  V  (  I  ) 


SOLVE  FOR  CR 


B.C.  AT  X=0 . DO 


A  (  1  )  =  -(2.D0+ALFA+OA*DXA2*GAMMA»C8(  1  )  ) 

C(  1  )  =  2.00 

0(1)  =  - ALFA  *CRP ( 1 )  -DXA2*0A*CAP( 1 ) *CBP(  1 ) /BETA 


LOWER  PART 


DO  16  I =2, NXA 
B ( I )  =  1.00 


A ( I )  =  -(2.00+ALFA+DA«OXA2*GAMMA*CB( I ) ) 

C(I)  =  1.00 

0(1)  =  -ALFA*CRP(I)  -DXA2*0A«CAP( I ) *CBP( I ) /BETA 


CONTINUE 


SPECIAL  CASES  NEAR  THE  INTERFACE  LINE 
B ( NXA  1  )  =  1.00 

A ( NX A  1  )  =  -( 2 . D0+ALFL+DA*DXL2»GAMMA *CB ( NXA  1  )  ) 

C ( NXA 1 )  =  1,00 

0 ( NXA 1  )  =  -ALFL*CRP(NXA1 )  -DXL2 *DA *  CAP ( NX A  1  ) *CBP ( NXA 1  ) /BETA 


B ( NXA 1  +  1  )  =  1.00 

A ( NXA 1  +  1  )  =  -(2.D0+ALFH+DA*DXH2*GAMMA*CB(NXA1+1 ) ) 

C ( NX A  1  +  1  )  =  1  .00 

0( NXA 1  +  1  )  =  -ALFH*CRP(NXA 1  +  1  ) -0XH2 +  DA * C AP ( NXA 1  +  1  )*CBP(NXA1  +  1  ) / BETA 


-  UPPER  PART 
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00  35  1=2. NX8 
S ( NXA I  +  I  )  =  I  .DO 
A(NXA 1  +  1 )  =  -( 2  . 


A(NXA1  +  I )  =  - ( 2 . D0+ALF8+DA*DXB2  * GAMMA  *CB  ( NXA 1 + I  ) ) 

C ( NXA 1 + I )  =  1.00 

D( NXA 1+1 )  =  -ALF8*CRP(NXA 1+1 ) -DA *DXB 2 *C AP ( NXA 1+1 ) »CBP(NXA1+I ) /BETA 


35  CONTINUE 


B.C.  AT  X=1 .DO 


8  (NX)  =  2.  DO 

A ( NX )  =  -(2.DO+ALFB+OA*DXB2*GAMMA*CB(NX)) 

D ( NX )  =  -ALFB*CRP(NX)  -DXB2 *DA *CAP ( NX ) »CB P ( NX ) /B ETA 
CALL  TRIDAGC 1 . NX , B , A . C . D , V ) 


DO  36  1=1 .NX 


36  CR ( I )  =  V ( I ) 


SOLVE  FOR  C8 


B.C.  AT  X=0 . DO 

A ( 1 )  =  -( 2 . D0+ALFA+0A*DXA2* (CA( 1 ) /8ETA+CR( 1 ) *GAMMA) ) 

C(  1  )  =  2. DO 

0(1)  =  -ALFA*C8P( 1 ) 


LOWER  PART 


DO  45  1=2. NXA 


B ( I )  =  1 .DO 

A ( I )  =  -( 2. D0+ALFA+DA+DXA2* (CA( I ) /BETA+CR( I ) *GAMMA) ) 

C(I )  =  1 .DO 

D(I)  =  -ALFA  *CBP ( I ) 


45  CONTINUE 

SPECIAL  CASES  NEAR  THE  INTERFACE  LINE 


B ( NXA 1 )  =  1 . DO 

A ( NXA 1 )  =  -(2.D0+ALFL+DA*DXL2*(CA(NXA1 ) /BETA+CR ( NXA 1 )*GAMMA) ) 
C ( NXA 1 )  =  1 . DO 
D ( NXA 1 )  =  -ALFL+CBP ( NXA 1 ) 


B ( NXA 1  +  1  ) 
A (NXA 1+1 ) 


-( 2 . D0+ALFH+DA*DXH2«(CA (NXA 1+ 1 ) /BETA 
+CR (NXA 1 + 1 ) *GAMMA )  ) 


C ( NXA 1 + 1 )  = 

D ( NX A  1  +  1  )  = 

CC  -  UPPER  PART 


=  -ALFH*CBP(NXA1+1 ) 


DO  46  1=2, NXB 
B (NXA 1+1 )  =  1 . 


A (NXA 1+1)  = 


-( 2 . D0+ALFB+0A*DXB2« (CA(NXA 1 +1 ) /BETA 
+CR (NXA 1 +1 ) ‘GAMMA )  ) 


C(NXA 1+1 ) 
D ( NXA 1 + I ) 


=  -ALFB*CBP(NXA1+I ) 


46  CONTINUE 


B.C.  AT  X= 1 . DO 


8  (  NX  )  =  2. DO 

A ( NX )  =  -( 2.DO+ALFB+DA+OXB2* (CA (NX) /BETA+CR (NX ) *GAMMA ) ) 

D ( NX )  =  -ALFB*CBP ( NX ) 

CALL  TR I DAG (1 ,NX,8,A,C,D,V) 


DO  47  I = 1 . NX 
47  CB( I)  =  V( I ) 
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CHECK  CONVERGENCE  FOR  THIS  TIME  STEP 


IF  (NS.EO. I )  GO  TO  9 
ERRL  =  0.00 
DO  11  1=1 .NX 

ERR  =  0A9S( (CA( I )-CAOLO( I ) )/CA( I )  ) 

IF(ERR .GT. ERRL)  ERRL  =  ERR 
1 1  CONTINUE 

170  FORMAT ( 15X ,  ' NS  =  '  , I  5 , 2X .  ' ERRL  =  '.012.5) 
IF ( ERRL. GT. TOR)  GO  TO  9 
GO  TO  12 
9  DO  56  1=1 .NX 
CAOLD(I)  =  C A (  I  ) 

56  CONTINUE 
5  CONTINUE 


INTEGRATION  TO  FIND  AVERAGE  CONCENTRATIONS 


12  I F ( NPR I  NT . NE .  1 )  GO  TO  1 
WR I TE ( 6 . 1 60 )  LT 
160  FORMAT (  /2X  , ' LT  =  '.15) 

CALL  ICSCCU(X, CA .NX.CC . 200 . IER) 

CALL  DCSQDU(X.CA, NX.CC, 200. 0.00, 1 . DO . Q . I ER ) 

Q 1  =  Q 

CAAVG  =  Q»CAO/CBOBAR 

CALL  ICSCCU(X,CR .NX.CC . 200 . IER) 

CALL  DCSQDU(X,CR, NX.CC. 200. 0.00. 1 .DO. Q, IER) 

CRAVG  =  Q*CBO/CBOBAR 

CALL  ICSCCUtX ,CB , NX.CC . 200 . IER) 

CALL  DCSQOU ( X . CB , NX , CC . 200 , 0 . DO , 1 . DO . Q . I ER ) 

C8AVG  =  Q*  CBO/ CBOBAR 

FA  =  1 . 00-(QA+QB)/QA*Q1 

CSAVG  =  CAOBAR/CBOBAR  -CAAVG-CRAVG 

RTIME  =  TIME’TSCALE 

Z  =  U  *RT I  ME 

WR I TE ( 6 . 150)  Z. RTIME. FA 

150  FORMAT ( 5X ,  ' Z  = '  . F9 . 4 . 2X ,  ' REAL  TIME  = '  . D 1 1 . 4 . 2X ,  *  FA  =',F8.5) 

:  OUTPUT  RESULTS  FOR  THIS  TIME  STEP 

WRITE (6. 1 40 )CA A VG. CBAVG. CRAVG. CSAVG. NS. ERRL 
140  FORMAT ( 5X , ' C A  =  ' , D 1 1 . 4 , 2X , ' CB  =  ' . D 1 1 . 4 , 2X , ' CR  ='.D11.4. 

*  2X , ' CS  = ' . 01 1 . 4 , 2X. 'NS  =  ' .I5.2X, 'ERR=‘ ,01 1 .5) 

IF(C8AVG.LE. 1 .00-6)  GO  TO  260 
1  CONTINUE 
260  STOP 


SUBROUTINE  TRIDAG  ( I F . L , B , A , C , 0 . V ) 
IMPLICIT  REAL*8(A-H,0-Z) 


SUBROUTINE  FOR  SOLVING  A  SVSTEM  OF  LINEAR  SIMULTANEOUS 
EQUATIONS  HAVING  A  TR 1 0 1  AGONAL  COEFFICIENT  MATRIX. 
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SUB-DIAGONAL.  DIAGONAL,  AND  SUPER -D I AGONA L  COEFFICIENTS 
ARE  STORED  IN  THE  ARRAYS  B , A ,  AND  C.  THE  COMPUTED 

SOLUTION  VECTOR  V(IF) . V(L)  IS  STORED  IN  THE  ARRAY  V. 

DIMENSION  BET  A ( 20 1 )  , GAMMA ( 20 1 ) 

DIMENSION  B ( 200 ) . A (200) . C(200) ,D(200) ,V(200) 


.  COMPUTE  INTERMEDIATE  ARRAYS  BETA  AND  GAMMA 

BETA (IF)  =  A ( I F ) 

GAMMA ( I F ) =D ( I F ) /BETA (IF) 

I FP 1  =  IF  +  1 
DO  1  1=  I FP 1 .  L 

BETA ( I ) =A ( I ) -  B( I )*C( 1-1 )/BETA( 1-1 ) 

1  GAMMA ( I )  =  ( D ( I )  —  S ( I ) * GAMMA (I-1))/BETA(I) 

.  COMPUTE  FINAL  SOLUTION  VECTOR  V  . 

V ( L )  =GAMMA(L) 

LAST  =L-IF 
DO  2  K= 1 . LAST 
I  =  L-K 

2  V(I)  =  GAMMA ( I )  -  C( I ) *V( 1+ 1 ) /BETA( I ) 


RETURN 


APPENDIX  A. 2 — Coalescence-Redispersion  (CRD)  Model 


CRD  MODEL  FOR  COMPETING  REACTIONS 


IMPLICIT  REAL*S(A-H.O-Z) 

DIMENSION  CA ( 505 ) , CB(505) ,CR(505) .CS(505) 

INTEGER  I R ( 2 ) 

DEFINE  SYSTEM  PARAMETERS 


TAUM=0 . 630D0 
QA  =  1 .DO 
QB  =  1 .DO 
CAO  =  0.2000-3 
CBO  =  0.2D-3 
RK1  =  7.3D6 
RK2  =  3.5D3 
VEL  =  6.30 

CAOBAR  =  CAO*QA/ (QA+QB) 

CBOBAR  =  CB0*Q8/( QA+QB ) 

BETA=C80BAR /CAOBAR 

GAMMA  =  RK2/RK1 

GAMMA  1  =  1  .DO/ ( I  .DO-GAMMA) 

WR I TE ( 6 . 100)  CAO, CBO, BETA, RK1 ,RK2, GAMMA 
100  F0RMAT(2X, 'CAO  =  ' , F9 . 5 , 2X . ' CBO  =  ' , Fg . 5 . 2X , ' BETA  =  ' , F8 . 3 . 2X 
«  , ' RK 1  =  ' ,01 1 .4.2X. 'RK2  =  ' . 0 1 1 . 4 , 2X , ' GAMMA  =  '.Dll. 4/) 

WR ITE(6,200)  VEL 
200  FORMAT ( 2X , ' VEL  =  ',F8.3/> 


DEFINE  NUMERICAL  PARAMETERS 

NED  =  300 
NED2=  150 
NED2 1 =  NED2+1 
FREQ  =  NED/TAUM 
DTAU  =  1. DO/FREQ 
INT  =  10 
DINT  =  DTAU/INT 
DINTH=  DINT 
NT  =  1400 

WR ITE(6,30Q)  NED, FREQ, INT, NT 
300  FORMAT ( 2X .  ' M  OF  EDDIES  =  '  , I  5 , 2X ,  ' FREQ 
*  . 1 5 , 2X .  ' #0F  TIME  STEP  =  ',15/J 

DEFINE  INITIAL  CONDITIONS  FOR  EACH  EDDY 


=  ’ , D1 1 . 4 , 2X, ' INT 


DO  10 

1=1 . NED2 

C  A  (  I  ) 

=  CAO 

CB  (  I  ) 

=  O.DO 

CR  (  I  ) 

=  0.00 

CS(I) 

=  O.DO 

10  CONTINUE 

00  20 

I  =NE02 1  , 

CA(I) 

=  0.00 

CB  C I  ) 

=  CBO 

v 
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CR { I )  =  0 . DO 
CS(I )  =  0 .DO 
20  CONTINUE 


CC 

CC 

CC 


START  THE  COALESENCE  .REDISPERSION  AND  REACTION  PROCEDURE 


T I  = 
DO  5 


0  .  DO 
LT  = 


1  .NT 


TI  =  TI+DTAU 


z 

VEL*TI 

GENERATE 

RANDOM  NUMBERS  FROM 

1  TO  NED 

2  CALL  TIME ( I TI ) 

DSEED  =  ITI *500 . DO 

CALL  GGUD (DSEED.NED, 2, IR) 

N 1  =  IRC  1 ) 

N2  =  I R  (  2 ) 

IFCN1.EQ.N2)  GO  TO  2 

500  FORMAT (216) 

CONCENTRATION  REDISTRIBUTION 

AFTER  REDISPERSION 

CA(N1 ) 

=  (CA(N1 )+CA(N2) )  / 

2.  DO 

CACN2) 

=  CA  (  N  1  ) 

CB(N) ) 

=  (CB(N1 )+CB(N2))  / 

2.  DO 

CBCN2) 

=  CBCN1  ) 

CRCN1 ) 

=  (CRCN1 )+CR(N2) )  / 

2.00 

CRCN2) 

=  CR ( N 1  ) 

CSCN1 ) 

=  (CSCN1 J+CSCN2) )  / 

2  .DO 

CSCN2) 

=  CSCN1  ) 

BATCH  REACTIONS  IN  THE  EDDIES 

DO  3  I 

= 1 .NED 

CAOLD 

=  CA(  I) 

CBOLD 

=  CBCI) 

CROLD 

=  CR (  I  ) 

CSOLD 

=  CSC  I ) 
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IF  (CAOL.D  .  EQ  .  0  .  DO  .  OR  .  CBOLD.  EO  .  0  .  DO  )  GO  TO  3 
DO  4  J=1 . INT 
DINT  =  DINTH 

IF(CA( I ) -RK1 *CA( I ) *CB( I ) *DINT. LE .0 . DO)  GO  TO  101 
KK  =  1 
GO  TO  103 
KK  =10 

DINT  =  DINTH/FLOAT(KK) 


103  00  50  JJ  =  1 , KK 
C AP=  C A ( I ) 

CBP=  CB ( I ) 

CRP=  CR ( I ) 

CSP=  CS(I) 

C A ( I )  =  CAP  -RK1 *CAP*CBP*OINT 
CAR  =  CA ( I ) / CAP 


'  *  m  % 
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CR (  I  )  =  (GAMMA  1 *( CAR» «GAMMA-CAR )  +CR P / C A P * C A R * • GAMMA ) « C A P 

CS(I)  =  CAP+CRP+CSP  -CA(I)-CR(I) 

CB { I )  =  CBP-  (CR(I)-CRP)  -2.D0*(CS(I)-CSP) 

C A ( I )  =  ( CA ( I ) +CAP )  /  2. DO 
CB ( I )  =  ( CB ( I ) +CBP )  /  2. DO 
CA ( I )  =  CAP  -RK1 «CA ( I ) *CB( I ) *DINT 

CAR  =  CA (  I  )  / CAP 

CR ( I )  =  ( GAMMA  1 *(CAR**GAMMA-CAR)  +CRP/ CA P » CAR » »GAMMA ) *C AP 
CS(I)  =  CAOLD+CROLD+CSOLD  -CA(I)-CR(I) 

C8(I)  =  C80LD-  (CR( I ) -CROLD)  -2 . DO* ( CS ( I ) -CSOLD ) 

50  CONTINUE 
4  CONTINUE 
3  CONTINUE 

IF(LT/20»20.NE. LT)  GO  TO  55 

CC - 

CC  CALCULATE  THE  AVERAGE  CONCENTRATIONS  FOR  ALL  EDDIES 

CC  - 

CAAVG  =  0 . DO 
CBAVG  =  O.DO 
CRAVG  =  O.DO 
CSAVG  =  O.DO 
DO  6  1=1 , NED 
CAAVG  =  CAAVG-CA(I) 

CBAVG  =  CBA VG+CB ( I ) 

CRAVG  =  CRAVG+CR ( I ) 

CSAVG  =  C  SA VG+CS ( I ) 

6  CONTINUE 

CAAVG  =  CAAVG  /NED 
FA  =  (CAOBAR -CAAVG)  /  CA08AR 
CBAVG  =  C8AVG  /NED  / CBOBAR 
CRAVG  =  CRAVG  /NED  /CBOBAR 
CSAVG  =  CSAVG  /NED  /CBOBAR 
CAAVG  =  CAAVG  /  CAO 


OELB  =  1 . DO-CBAVG-CRAVG-2 . 00*CSAVG 
WRITEC6.400)  Z.TI.FA 

400  F0RMAT(/2X,  '  Z  =  '  .F9.4.2X,  'REAL  T I  ME  =  '  , D 11 . 4 . 2X ,  ■ FA  =',F8.5) 
WRITE (6. 600)  C AAVG. CBAVG. CRAVG. CSAVG. DEL A , DELB 
600  FORMAT ( 2X , '  CA=  ' .01 1 .4.2X. 'CB=' , D 1 1 ,4,2X, '  CR=  ' , D 1 1 .4,2X, 'CS=' 
*  011 .4 , 2X. 'DELA= ' ,D11.4,2X.'0ELB='.011.4) 

55  CALL  TIMREMf ITJ) 

IF( ITJ .GE. 200)  GO  TO  5 
WRITEC 16. 140)  Z.TI 

140  FORMAT ( 2X , ' Z  =  ' . F9 . 4 , 2X . ' REAL  TIME= ' , D 1 1 . 4 ) 

WRITEC 16. 145)  ( C A ( I ) , 1=1 .NED) 

WRITEC  16. 145)  (CBC I )  ,  1  =  1  .NED) 

WRITEC  16,  145)  (CRCI  )  ,  1  =  1  .NED) 

WRITEC  16,  145)  (CSC  I )  . I  =  1  . NED) 

145  FORMAT ( 70 1 1.4) 

5  CONTINUE 
STOP 
END 


> 


APPENDIX  A. 3 — Three-Environment  (3E)  Model 
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c - 

C  THREE  ENVIRONMENT  MODEL  FOR  COMPETING  RE AT  I ONs 

c - 

IMPLICIT  REAL*8  (A-H.O-Z) 

EXTERNAL  DFUN 

DIMENSION  V(8 ,4)  ,  VMAX(4)  ,WK(84. 1)  ,ERR0R(4)  ,CEXIT(4) 
COMMON /COM  1 / R  ST , Q 1  .Q2.PHI  ,  C I  N  (  4  )  , DM  1  , DM2 . TAU . VLCAP 
READ  (5.10)  Q! . 02 . CAO , CBO . TEMP 
READ  (5,20)  TAUM 

READ  (5,30)  HMAX.H.HMIN.EPS.TSTOP 
READ  (5.31)  RK 1 , RK2 
31  FORMAT ( 4X , D9 . 2 . 5X , D9 . 2 ) 

U=7 . 6000 

TAU=24 . 0D0*2 . 54D0/U 
HI  N  =  H 

RST=TAU/TAUM 
PHI =Q2/Q1 

CB0BAR=CB0*Q2/ (Q1 +Q2) 

CAOBAR=CAO*Q1 / (Q1+Q2) 

BETA=CB08AR/CA0BAR 

WRITEI6.40) 

WR I TE ( 6 , 50 ) 

WR I TE ( 6 , 60 )  Q 1 .Q2.CA0BAR.CB0SAR, TEMP 
IF  (BETA. GT. 2.0)  GO  TO  200 
CA 1 0  =  (  1  .0+Q2/Q1 J/BETA 
CB20= ( 1 .0+Q1/Q2) 

DLCAOB- 1 . O/BETA 
DLCBOP - 1 .0 
COBAR=CBOBAR 
GO  TO  210 

200  CA 1 0= ( 1 . 0+Q2/Q 1 ) 

CB  20  =  BET  A  * ( 1  .0  +  Q1/Q2) 

DLCAOB= 1 . 0 
DLCBOB=BETA 
COBAR=CAOBAR 
210  CONTINUE 

DM  1 =  RK1 *COBAR*TAU 

DM2  =  RK  2  *DM 1 /RK 1 

TRXN 1 =T AU / DM  1 

TRXN2=1 . 0/RK2/CA0BAR 

V ( 1 . 1 ) =CA 1 0 

V( 1 . 2 ) =CB20 

V( 1 . 3 ) =0 . ODO 

V( 1 , 4 ) =0 . 000 

CA=CA 1 0«C0BAR 

C8=CB20*C0BAR 

WR I TE ( 6 , 70 )  R  ST . BETA 

WR I TE ( 6 , 75 )  TAU . TAUM , TRXN 1 . TRXN2 

WR I TE ( 6 . 80 )  HMAX. H.HMIN.EPS 

WR I TE ( 6 . 90 )  CA.CB 

1 0  FORMA T(4X . F6. 1 .6X.F6. 1 .6X.09.2.6X.D9.2.5X.F6.2) 

20  F0RMAT(7X , F9 . 3 ) 

30  FORMAT(4(5X,D10.3) .6X.F7.3) 

40  FORMAT ( IX. 1 20 ( '•')/) 


vv j v, 


W.v-^ W .V  'W  ■ 


-_'  -.  --  >. 


Lw  vw  tw  »“  ■ 
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50  FORMAT  (  40X  ,  '  SELECTIVITY  IN  PLUG  FLOW  R  E  ACTOR '/ 35X  , 

1  ' SOURNE-S  KINETICS  :  THREE  ENVIRONMENT  MODEL'/) 

60  FORMAT ( 1  OX ,'Q1  =  '.F5.1,'  CC/S '  ,5X ,  ' 02= '  . F6 . 2 ,  '  CC/S',5X. 

1  / 1  OX ,  ' CAOB AR  =  '  ,D10.3,  'MOL/L'  .  5X.  ' CBOBAR= '  ,010.3,  'MOL/L'  , 

2  .5X.-TEMP  (DEG  C)=',F6.2/) 

70  FORMAT (10a, 'T/TAUM=  ' , 0 1 0 . 3 . 5X , ' BETA=  '.DIO. 3/) 

75  FORMAT( 15X, 'TAU=' ,F7.3, '  S' ,5X, *TAUM=' ,010.3. '  S ' . 5X , 

*  ' TRXN 1 = ' . D 1 0 . 3 , '  3' ,5X, ' TRXN2= ' ,010.3, '  S'/) 

80  FORMAT ( 1  OX,  'HMAX= '  ,010.3,5X,'H=',010.3,5X.  ' HMIN= '  ,010.3, 

1  5X , ' EPS= ' .010.3/) 

90  FORMAT( 15X.  'CA=  ',09.2,'  MOL/L ', 5X ,' CB=  '  , D9 . 2  ,  '  MOL/L'/) 

WRITEC6.40) 

WR I TE ( 6 , 1 30 ) 

WRITEC6. 140) 

DO  100  1=1,4 
CIN( I)=V( 1 . I ) 

YMAX ( I ) = 1 . 000 
100  CONTINUE 

T=HIN/100.000 
N  =  4 
MTH  =  2 
MAXDER=6 
J  START  =  0 
H  =  HIN 

TKEEP=0 . 1 000-03 
TSTEP=0 . 002 
D I ST=0 . 000 

900  CALL  OVOGERCOFUN.Y.T.N.MTH.MAXOER, JSTART . H , HM I N . HMAX , 

*  EPS, VMAX, ERROR, WK, IER) 

IF(T.LT.TKEEP)  GO  TO  900 
TKEEP=TKEEP*TSTEP 
IF(TKEEP.GE.0.01 )  TSTEP=0.01 
IF(TKEEP.GE.O. 1 )  TSTEP=0 . 1 
IF  (RST*T.GT. 150. )  GO  TO  250 
THRI=DEXP(-RST*T) 

GO  TO  260 
250  THRI=0.0 
260  CONTINUE 
QS  =  Q1 

DO  110  L=1 .4 
I F ( L . EQ . 2 )  QS=02 

1  1 0  CEXIT( L)  =  (VLCAP*V(  1  , L ) +  THR I  * ( QS *C I N ( L ) ) ) / (01 +Q2) 

OLCA=CEX I T ( 1 ) 

IF  (BETA . LE. 2 . 000)  DLCA=DLCA«BETA 

DLCR=(  1  . 000/(  1  . 000-RK2/RK 1  ) )*( ( OLC A ) * « ( RK 2 / RK 1  ) -DLCA ) 

IF  (BETA . LE. 2 . GOO )  DLCR =DLCR / BETA 
DLCA  =  C  EX  I T (  1  ) 

OLCS=DLCAOB-DLCA-OLCR 
OLC8=OLCBOB-DLCR-2 . 000* OLC S 

CHECK  1  =  (OLC AOB  )-CEXIT(  1  )-CEXIT(3)-CEXIT(4) 

CHECK2  =  DLCB0B-CEXIT(2) -CEXIT(3) -  2 . 000*CEX I T ( 4 ) 

SELECT  =  2 . 0D0*CEXIT(4) / (CEXIT( 3) *2 . 0*CEXIT(4) ) 

AT  =  T*TAU 
01  ST  =  T*U*  TAU 


•;v\ 


v 


A 


■A^a 


y-:- 


i 


m 


*/  v"  Oi 

*j 

/  /  . 
AAAj 


C.'jCm  AjL 


WRITEC6, 120)AT.DIST, (CEXIT(JK) . JK=1 .4) . DLCB , DLCR . DLCS . SELECT 

*  , CHECK  1  .CHECK2 

IF  (CEXITC2) .LT.O. ID-2)  GO  TO  160 
IF  (T.LT.TSTOP)  GO  TO  900 
C 
C 

160  WRITE(6, 1 20 ) AT , D I  ST , (CEXIT(JK)  ,JK=1  ,4)  ,  DLCB , DLCR , DLCS . SELECT 

*  ,  CHECK  1  . CHECK2 
WRITE(6.40) 

SELECT=CEXIT(3) /CEXIT(4) 

C0NVA=  ( DLC AOB-CEXI T (1))/DLCA0B 
CONVB=  (DLCB0B-CEXIT(2) )/DLCB0B 
WRITE(6,170)  SELECT, CONVA . CONVB 
WR ITE (6,40) 

1 20  FORMAT ( 2X . F6 . 3 . 8 ( 4X . F6 . 3 ) . 3 ( 2X , 09 . 2 ) ) 

130  F0RMAT(32X, ' SCALED  CONCENTRATIONS ', 1 2X ,' MAX .  MIXED'/) 

140  FORMAT (3X, 'TIME(S) ' ,3X, 'DIST.(CM) ' ,2X,  ' CA ' .  9X.'CB'.  BX , 

*  'CR',  8X, 'CS 7X, 'CBM' . 7X. ' CRM' , 7X . ' CSM ', 7X . ' SELECT ' . 

*  7X, 'DELI ' ,6X, ' DEL2 ' ) 

170  FORMAT ( 10X,  'SELECTIVITY:  CR/CS  =  '  . D 1 0 . 3 .  1  OX .  ' CONVER S I ON  OF  A 

*  =', DIO. 3. 10X. 'CONVERSION  OF  B=',D10.3/) 

STOP 

END 

SUBROUTINE  OF UN ( VP , TP , M , OV . PW , I ND ) 

IMPLICIT  REAL'S  (A-H.O-Z) 

DIMENSION  VP (8, 4) ,PW(4,4) ,DV(4) ,RXN(4) 

COMMON /COM  1/RST,Q1,Q2.PHI.CIN(4).DM1,DM2,TAU,VLCAP 
IF  (RST*TP,GT.75. )  GO  TO  5 
AXP  =OEXP(-RST*TP) 

RAT=AXP/ ( 1 , ODO-AXP) 

GO  TO  6 

5  AXP=0 . ODO 
RAT= 1 .ODO 
6  CONTINUE 

VLCAP=(Q1 +Q2 ) » ( 1 . ODO-AXP) 

C 

RATE  1 =  DMl'YPf 1  ,  1 )*VP( 1  ,2) 

RATE2=  DM2*YP( 1 , 2) *VP( 1 , 3) 

C 

RXN ( 1 )=-RATE1 
RXN(2)=-RATE2  -RATE  1 
RXN(3)=RATE1-RATE2 
RXN(4)=RATE2 


50 

60 


DY(1)=RST*RAT»(CIN(1)*(1.0DO/(1.0DO+PHI))-YP(1.1))+RXN(1) 
DV(2)=RAT*RST*(CIN(2) *(PHI/( 1  . ODO 'PH  I  ) )-YP( 1  , 2) ) +RXN ( 2) 

DY ( 3 ) =  -RAT'RST*YP( 1 ,3)+RXN(3) 

0V(4)=  -RAT*RST*VP( 1 ,4)'  RXN(4) 

IF( TP . LE . 6 . 0D-03 )  GO  TO  60 
CO  50  1=1,4 

IF  (YP(1 ,1) ,LE. 1 .00-15)  DY( I )=0V( I ) / 10 .ODO 
IF  (VP( 1 , I ) ,LE . 1 .00-16)  DY( I ) =DY( I )/ 100 .ODO 
IF  ( VP( 1 . I ) . LE . 1 .00-20)  DV( I ) =0V( I ) / 1 00 . ODO 
RETURN 


ENO 


V.v.v.v.V) 


APPENDIX  A. 4 — Interaction  by  Exchange  with  Mean  (IEM)  Model 
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cc - 

C  IEM  SIMULATION  PROGRAM 

cc - 

IMPLICIT  REAL*8(A-H.O-Z) 

EXTERNAL  DFUN 
REAL’S  MTC 

DIMENSION  WK( 138) ,CMAX(6) 

COMMON/  CONC/  C(8.6) 

COMMON /CON 2/  CAAVG , CBAVG . CRAVG . CSAVG 
COMMON/PARM/  MTC  .RK1.RK2 

CC  —  DEFINE  SYSTEM  PARAMETERS  - 

TAUM=2.00 
TEMP=27. 3 
U  =7.60 

RK1  =  7.3D+6 

RK2=1 . 1 1D09»OEXP(-3.957D03/(TEMP+273) ) 

C  RK2=1 .9301 2*DEXP( -5 .989003/ (TEMP+273 ) ) 

01=131.0 
Q2=  1 .31 

CAO  =  0.300D-03 
CBO  =  30.080-03 
CAOBAR  =  CAO  *Q 1 / ( Q 1 +Q2 ) 

CBOBAR  =  CB0*Q2/ (Q1+Q2) 

PHI=Q1 /Q2 

SM  =  PHI/ ( 1 . DO+PH I ) 

GAMMA  =  RK1/RK2 
MTC=0.500/TAUM 
BETA  =  CBOBAR/CAOBAR 
WRITE(6. 222) 

222  FORMAT( IX,  1 20  C  )/) 

WRITEC6.223) 

223  FORMAT  (20X.  'IEM  SIMULATION') 

WR ITE(6,500)  CAO. CBO, BETA, RK1 ,RK2, GAMMA 

500  F0RMAT(/2X, 'CAO  =  ' , 0 1 1 . 4 , 2X , ' CBO  =  ' , D 1 1 . 4 , 2X , ' BETA  =  ',D1 

*  , ' RK 1  =' .011 ,4,2X. 'RK2  =  ' , D 1 1 . 4 , 2X , ' GAMMA  =  ',011.4/) 

WRITE  (6,224)  TAUM.Q1.Q2 

224  F0RMAT(/2X, 'TAUM=  ' , F7 . 3 , 3X , ' Q 1 =  ' , F7 . 3 . '  CC/ S ' , 3X , ' 02=  ', 

*  F  7 . 3 ,  '  CC/S'/) 

WRITEC6, 222) 

CC  —  DEFINE  NUMERICAL  PARAMETERS  - 


T  = 

0.00 

N  = 

6 

MTH  = 

2 

MAXDER= 

6 

J  START  = 

0 

H  = 

0.010-3 

HMIN  = 

1  .  00-  1  1 

HMAX  = 

0.0100 

EPS 

1 

a 

o 

TGOAL  = 

0.  1 

IND 

0 

CC  -  INITIAL  CONDITIONS 


1  . 4 , 2X 


cc 


C  (  1  .  3  }  =  0.00 
C( 1 .4)=  0.00 
C(  1  ,5)=  CBO 
C(1 .6)=  0.00 

DO  1  I  =  1  .  N 
1  CMAX ( I )  =  1.00 

CC  INTEGRATION  BY  GEAR  PACKAGE 
CC 

ICOUNT  =  0 

20  ICOUNT  =  ICOUNT  +1 

5  CALL  0V0GER(DFUN,C,T,N,MTH,MAXDER, J START , H . HMI N . HMAX , EPS , CMAX , 

*  ERROR ,WK. IER) 

IF(IER.GE.33.AND. IER.LE.36)  GO  TO  60 
CAAVG  =  SM*C ( 1.1)  +  ( 1 .DO-SM)*C( 1 .4) 

CBAVG  =  SM*C(1,2)  + ( 1  . DO- SM ) »C ( 1  , 5  ) 

CRAVG  =  SM*C  ( 1,3)  +( 1 .DO-SM)*C{ 1 .6) 

IF(ICOUNT/50*50.EQ. ICOUNT)  GO  TO  10 
GO  TO  20 

10  FA  =  1.00  -CAAVG/CAOBAR 
Z  =  T*U 

CAAVG 1  =  CAAVG/CBOBAR 
CRAVG 1  =  CRAVG/CBOBAR 
CBAVG 1  =  CBAVG/CBOBAR 

CSAVG1  =  CAOBAR/CBOBAR  -CAAVG 1 -CRAVG 1 
WR I TE ( 6 , 100)  Z.T 

WR I TE ( 6 , 200 )  CAAVG 1 , CBAVG 1 , CRAVG 1 .CSAVG1 

TGOAL  =  TGOAL+I.O 

IF  (CAAV01 . GT . 1 .OD-06)GO  TO  59 

CC 1 , 1 ) =0.000 

C( 1 , 4 ) =0 . 000 

59  IFCCBAVG1 .GT.O. 10-03)  GO  TO  20 
GO  TO  30 

60  WRITEC6.300)  IER 

300  FORMAT ( 5X , 'IER  =  '.15,' -  PROGRAM  WAS  TERMINATED'/) 

100  FORMAT ( 5X , ' Z  =  ' .F9.4.2X. 'TIME3' ,D1 1 .4) 

200  FORMAT ( 5X , ' CA  =  ' , 0 1 1 . 4 , 2X , ' CB  =  ' , D 1 1 . 4 . 2X , ' CR  =  ' . 0 1 1 . 4 , 2X 

*  . 'CS  =  ' ,01 1 .4/) 

30  CONTINUE 

STOP 

END 

CC 

CC 

SUBROUTINE  DFUNCCP.TP.N.DC.PW. IND) 

IMPLICIT  REAL«8(A-H,0-Z) 

REAL*8  MTC 

DIMENSION  CP ( 8 , 6 ) ,0C(6) , PW( 6 , 6 ) 

COMMON/  CONC /  CC8.6) 

COMMON /CON 2/  CAAVG . CBAVG . CRAVG . CSAVG 
COMMON/PARM/  MTC  .RK1.RK2 


CC  -  FIRST  EDDY  - 

OC ( 1  )  =  MTC* (CAAVG-CPC  1  .  1  ) )  -R K 1  *  CP (  1  .  1  ) *  CP (  1  , 2 ) 


-  /- 


»V.>, 


£3 


DC ( 2) =MTC» ( CBAVG-CP ( 1 .2) )  -  RK  1 *CP( 1 . 1 )*CP( 1 .2)  -RK2«CP( 1 . 2) *CP( 1 ,3) 
DC ( 3) =MTC* ( CRAVG-CP (  1  ,3))+RKl*CP(1  .  1 ) »  C  P (  1  .2)  -RK2»CP( 1  , 2 ) *CP( 1  .3) 
-  SECOND  EDOV  - 

DC  C  4 ) =  MTC*(CAAVG-CP(  1 ,4) )  -RK1 *CP  ( 1 ,4)*CP{  1  .5) 

DC(5 ) =MTC* ( CBAVG -CP( l,5))-RKt«CP(l,4)*CP(1,5)  -RK2*CP(1,5)*CP(1,6) 
DC(6)=MTC*(CRAVG-CP( 1 ,6) )+RK1 *CP( 1 . 4)«CP( 1 .5)  -RK2*CP( 1 , 5)*CP( 1 .6) 
RETURN 


APPENDIX  A. 5 — four -Environment  (4E)  Model 


- - 

C  FOUR  ENVIRONMENT  MOOEL 

- - 

IMPLICIT  REAL«8  (A-H.O-Z) 

EXTERNAL  OFUN 

DIMENSION  V ( 8 , 8 ) .VMAX(8) , WK ( 200 , 1 ) . ERROR ( 8 ) , CEXI T ( 4 ) 
COMMON/COM  1 /  V1LCAP.V2LCAP.RST.Q1  .02, DM1  ,DM2,CIN(8) 
READ  (5,10)  TEMP , C80 , CAO 
READ(5 , 20 )  TAUM 

READ (5,30)  HMAX , H , HMI N , EPS .TSTOP 
READ(5 , 32 )  Q1 ,Q2 
32  F0RMAT(4X,F9.2,5X.F9.2) 

RK 1 =7 . 30D06 

RK2=1 . 1 1D09*0EXP( -3 .957003/ (TEMP+273) ) 

C  RK2=1 ,93D12*DEXP(-5.989D03/( TEMP+273 ) ) 

U=7 . 60D0 

TAU= 19 . 000*2 . 54D0/U 

RET=2.54DO*2.0DO*U/0.0100 

HIN=H 

RST=TAU/TAUM 
PHI =Q2/ ( Q 1 +Q2) 

BETA=  CBO*Q2/CAO/Q1 

CBOBAR=CBO*PHI 

CA0BAR=CA0*Q1/(Q1+Q2) 

WRITE( 6 , 130) 

WR I TE ( 6 , 55 ) 

WRITE! 6 , 65 )  Q  1  .  Q2 . RET . CAOBAR , CBOBAR , TEMP 
IF  (BETA. GT. 2.0)  GO  TO  200 
CA10=( 1 . 0+Q2/Q1 ) / BET A 
CB20=( 1 .0+01 /Q2) 

DLCA0B=1 . O/BETA 
DLCB0B=1 .0 
COBAR=CBOBAR 
GO  TO  210 

200  CA 1 0= ( 1 . 0+Q2 /Q 1 ) 

CB2Q=BETA*(  1  .0+Q1/Q2) 

0LCA0B=1 .0 
DLCBOB  =  BET  A 
COBAR=CAOBAR 
210  CONTINUE 

DM  1 =  RK1 «COBAR*TAU 
DM2=RK2*DM1 /RK 1 
TM=TAU/RST 
TRXN1=TAU/DM1 
TRXN2= 1 . 0/RK2/CA0BAR 
C 

C  INITIAL  CONCENTRATIONS 

C 

DO  96  1=1,8 
Y( 1 . I ) =0 . 0D0 
96  YMAX ( I ) = 1 .000 
Y( 1  ,  1 ) =C A  1 0 
Y ( 1 , 6 ) =CB20 


CIN(I)=V( 1 .1 ) 

CA=CA10*C0BAR 
CB=CB20* COBAR 
WRITE  (6,71)  RK1 ,RK2 

FORMAT ( 5X , ' K 1 =  ' . 09 . 2 . '  L/MOL/S'.  5X,'K2=  ',09.2,'  L/MOL/S1/) 
WR ITE(6,70)  TAUM.RST.BETA 
WRITE(6,75)  TAU. TM.TRXN1 ,TRXN2 
WRITE(6,80)  HMAX, H.HMIN.EPS 
WR I TE (6,90)  CA.CB 
WRITEC6, 130) 

WRITE  (6. 135) 

WR I TE ( 6 , 140) 

DEFINE  PARAMETERS  FOR  OVOGER 

T=HIN/100.0DO 
N=8 
MTH  =  2 
MAXDER  =  6 
J  START  =  0 
H=HIN 

TKEEP=0. IDO 
TSTEP=0 . 002 
DIST=0.0D0 

CALL  DVOGER ( DFUN , Y , T , N , MTH . MAXDER , JSTART , H , HMI N , HMAX . 

♦  EPS, YMAX, ERROR. WK. IER) 

IF(T.LT.TKEEP)  GO  TO  100 
TKEEP=TKEEP+TSTEP 

I F ( TKEEP . GE .0.1)  TSTEP=0 . 1 
IF  (RST*T.GT. 150. )  GO  TO  250 
THRI=DEXP( -RST*T) 

GO  TO  260 
THR 1=0 . 0 
CONTINUE 
DO  110  L=1  .4 

CEXIT(L)=(V1LCAP*Y( 1 ,L)»V2LCAP*Y(1 . L+4 ) +THR I « ( Q 1 » 

*  CIN(L)+Q2*CIN(L+4) ) )/(Ql+Q2) 

DLCA  =  CEXIT(  1 ) 

IF  (BETA. LE. 2.000)  DLC A  =  OLC A  *  BETA 

DLCR= ( 1 ,0D0/( 1 .0D0-RK2/RK1 ))»( (OLCA ) ** (RK2/RK 1 J--DLCA) 

IF  (BETA. LE. 2. 000)  DLCR=DLCR/BETA 
DLCA=CEXIT ( 1 ) 

OLCS=DLCAOB-DLCA-DLCR 
DLCB=OLCBOB-DLCR-2 . 000* DLC S 

CHECK1=(DLCA0B  )-CEXIT( 1 ) -CEX I T ( 3 ) -CEX I T ( 4 ) 
CHECK2=DLC80B-CEXIT( 2 ) -CEXIT( 3 ) -2 . 0D0*CEXIT(4) 

SELECT  =  2 . 0D0*CEXIT(4) / (CEXIT(3 )  +  2 .0*CEXIT(4) ) 

AT=T*TAU 
01  ST  =  T*U*  TAU 

WRITE(6, 1 20) AT . 01  ST .  (CEXIT( JK)  . JK= 1  , 4)  , DLCB , OLCR .DLCS , SELECT 
«  , CHECK  1  , CHECK2 

IF  (CEXIT(2) . LT.O. 10-3)  GO  TO  160 
GO  TO  100 


A-  AT ^ 


160 


WR I TE ( 6 ,  I  20) AT, GIST,  ( CEXIT( JK) , JK=1 , 4) ,  DLCB, DLCR , DLCS. SELECT 

*  , CHECK  1  , CHECK2 
WR I TE ( 6 , 130) 

SELECT=CEXIT(3)/CEXIT( A) 

C0NVA=  (DLCA0B-CEXIT( 1 ) ) /DLCA08 
C0NV8=  ( DLCBOB-CEXI T ( 2 ) ) /DLCBOB 
WR I TE (6, 170)  SELECT , CONVA . CONVB 
WRITEC6. 130) 

C 

C 

10  FORMAT ( 5X , F4 . 1 ,6X,09.2,8X,09.2) 

20  FORMAT ( 7X , F9 . 3 ) 

30  FORMAT ( 4 ( 5X , D 1 0 . 3 ) , 6X , F5 . 2 ) 

55  FORMAT ( 40X , '  SELECTIVITY  IN  PLUG  FLOW  REACTOR ' /35X . 

1  'BOURNE-S  KINETICS  :  FOUR  ENVIRONMENT  MODEL’/) 

65  FORMAT (10X,'Q1=',F5.1,'  CC/ S ' , SX . ' Q2= ' . F6 . 2 ,  '  CC/S'.5X, 

1  'RER=' ,F10. 1 .5X/10X, 'WT.ED  STREAM  A='.D10.3. 

2  '  MOL/L' ,5X, "WT.ED  STREAM  B=',D10.3,'  M0L/L',5X, 

3  'TEMP  (DEG  C)=" ,F6. 2/) 

70  FORMAT! 10X, ' TAUM= ' ,010.3, 5X , 

1  "RSTAU="  , DIO. 3, 5X,  'FLOW  WT.ED  ( CB/CA )  =  '  , D 1 0 . 3/ ) 

75  FORMAT! 15X,  ' TAU= '  .  F7 . 3,  '  S ' , 5X , ' TM= ' ,D10.3, '  S '  ,  5X  , 

«  ' TRXN 1  =  '  , D 1 0 . 3 , '  S'  ,5X, 'TRXN2='  ,010.3. '  S'/) 

80  FORMAT! 1  OX ,  ' HMAX  = '  , D 1 0 . 3 , 5X ,  ' H= '  . D 1 0 . 3 . 5X ,  ' HM I N= '  . D 1 0 . 3 . 

1  5X. 'EPS=' ,010.3/) 

90  FORMAT!  1 5X , ' CA=  ',09.2,  '  MOL/L ', 5X ,' CB=  ' . D9 . 2 , '  MOL/L'/) 

120  FORMAT! 2X.F6. 3,4X,F7.2, 7(4X.F6.3) ,3!2X,D9. 2) ) 

130  FORMAT! IX, 129! '-' )/) 

135  F0RMAT(32X, 'SCALED  CONCENTRATIONS' , 12X, 'MAX.  MIXED'/) 

140  FORMAT! 3X. ' TIME! S) '. 3X , ' OIST. (CM) ', 2X,  ' CA ' ,  9X,'CB',  8X, 

*  'CR'.  8X,  'CS'  ,7X,  'CBM'  ,7X,  'CRM'  ,7X,  'CSM'  ,7X,  'SELECT'  , 
tt  7X  ,  'DELI  '  ,  6X  ,  '  DEL2  '  ) 

i/O  FORMAT (  1 UA .  SELECT! V  I T\  :  CR/C3  - '  , U i u . J ,  1  OX .  ' CONVERS I  ON  OF  A 

*  =', DIO. 3.  10X,  'CONVERSION  OF  B=',D10.3/) 

STOP 

ENO 

C 

C 

SUBROUTINE  OFUN(VP,TP,M,DY,PW.INO) 

IMPLICIT  REAL*8  (A-H.O-Z) 

DIMENSION  VP (8 , 8) . PW(8 ,8 ) ,0Y(8) , RXN(8 ) 

COMMON / COM  1/  V 1 LCAP . V2LCAP , RST . 0 1  . 02 , DM  1  , DM2 . C I N ( 8 ) 

IF  (RST*TP.GT.75. )  GO  TO  5 
PSI 1 =DEXP( -2 . 0*RST*TP) 

PSI2=DEXP(-RST*TP) 

GO  TO  6 
5  PSI1=0.0 

PS  I  2  =  0 .0 
6  CONTINUE 

V1LCAP=(Q1/2.0)*(1 .0-( 1 .0-02/Q1 )*PSI 1 +02 /Q 1 *( 1 .0-2.0* 

*  PS  I  2  )  ) 

V2LCAP=(Q2/2 .0)*(1.0-(1.0-Q1/Q2)*PSI1«-Q1/Q2»(1. 0-2.0* 

*  P  S I  2  )  ) 


RATE  1 =  DM1 *VP( 1 . 1 ) *VP( 1 . 2) 
RATE2=  DM2*VP( 1 ,2)*VP( 1 ,3) 
RATE3=  DM 1*YP( 1 ,5)*YP( 1 ,6) 
R ATE4=  DM2*YP( 1 ,6)*YP( 1 ,7) 

RXN ( 1 )=-RATE1 
RXN(2)=-RATE2  -R ATE  1 
RXN(3)=RATE1 -RATE2 
RXN ( 4 ) =RATE2 
RXN(5)=-RATE3 
RXN(6)=-RATE4-RATE3 
RXN(7)=RATE3-RATE4 
RXN ( 8 ) =  RATE4 


DO  10  K=1  ,4 

DV(K)=RST* ( (Q1 /VI LCAP) * (CIN( K) -VP{ 1 ,K) ) *PSI 2+ (V2LCAP/ 
V 1 LCAP ) * ( YP ( 1 , K+4)-YP( 1 ,K) ) )+RXN(K) 

DO  20  K=5 , 0 

OV(K)=RST* ( (Q2/V2LCAP) » (CIN(K) -VP( 1  .K) ) *  PS  I  2+ ( V 1 LCA P/ 
V2LCAP ) * ( YP ( 1 , K-4)-VP( 1 . K ) ) )+RXN(K) 

IF(TP.LE.6.0D-03)  GO  TO  60 
DO  50  1=1.0 

IF  (YP(I.I).LE. 1 .00-15)  DV(I)=DV(I)/10.0D0 
IF  (YPd.D.LE.  1  .0D-16)  DY  (  I  )  =OY  (  I  )  /  1  00  .  ODO 
IF  (YP( 1 , I) .LE. 1 .00-20)  DY ( I ) =OY ( I ) / 1 00 . ODO 
RETURN 


APPENDIX  B 


SAMPLE  CALCULATIONS  FOR  THE  TURBULENCE  MEASUREMENTS 


The  phasic  output  of  the  pulsed  Doppler  Ultrasound  Velocimeter 
is  sampled  and  processed  by  the  microprocessor  and  mean  and  the  root- 


mean-square  (r.m.s.)  value  of  the  signal  are  displayed  in  mV  on  the 
microprocessor  screen.  These  values  are  used  in  the  equation  below, 
which  is  based  on  Equation  (8)  in  the  text,  to  calculate  the  mean  and 
r.m.s.  velocities.  Basically,  it  incorporates  the  conversion  of  the 
Doppler  shift,  Af^,  by  the  zero-crossing  frequency  meter  to  an  analog 
voltage: 


u (cm/ s) 


(1  KHz/0.25  volt)  •  (u  cm/s)  •  V(mV) 
_ s _ 

(f  KHz)  (2  cos  0)  (1000  mV/volt) 


(B.l) 


For  the  fully  developed  pipe  flow  experiments, 


u  =  153100  cm/s  (value  obtained  from  CRC  handbook) 
s 

e  =  43° 

f  =  104  KHz 
o 

Then, 

u(cm/s)  =  41.87  x  10  ^  V  (B.2) 


can  be  used  to  calculate  the  mean  velocities  by  substituting  the  mean 

output  voltage  (in  mV)  for  V.  For  example,  at  Re  =  8029,  the  mean 

voltage  V  =  1315  mV  and  correspondingly,  the  experimental  mean  velocity 

u  =55.1  cm/ s . 
exp 

The  calculation  is  a  little  different  for  the  r.m.s.  velocity.  If 
the  turbulence  field  is  isotropic,  the  measured  fluctuating  velocity 
would  be  the  same  regardless  of  the  angle  at  which  the  probe  observes 
the  flow.  This  requires  that  for  the  r.m.s.  velocity  computations, 
9=0.  Equivalently, 


The  colorimeter  circuit  gives  an  analog  output,  V  ,  which  is 

proportional  to  £n(IX/IX) .  IX  is  the  constant  internal  reference 

o  o 

intensity.  According  to  Equation  (14)  (see  Chapter  IV),  however,  the 
time  average  should  be  taken  only  after  a  natural  logarithm  of  the 
signal  is  taken,  i.e.,  £n(l“/IX).  Since  Equation  (16)  (i.e.,  AX  ~ 
£n(IX/lX))  has  been  used  to  analyze  the  experimental  data,  it  is 
necessary  to  estimate  the  error  introduced  by  such  an  approximation. 

If  Equation  (12)  (cf.  Chapter  IV)  is  rewritten  as 

IX 

-j  =  exp  (er.>ICr  +  GgS-Cg)  (C.l) 

and  the  concentrations  as 


CR  = 


CR 


(C.  2) 


and 


"  5S  +  cs 


(C .  3) 


(overbar  denotes  the  time-mean  value  and  prime  denotes  the  fluctuating 
component),  then  the  following  equation  is  obtained: 


o 


exp(eRJ,CR 


exp(eR£CR 


esics) 


(C.4) 


Taking  natural  logarithm  of  two  sides  of  the  above  equation,  after  they 
are  time  averaged,  results  in 


Zn(~)  =  ch  C  +  e*fcC_  +  £n  ( exp  (cX£C '  +  eX£C')]  =  AX  +  A. 


(C.  5) 
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aX  ■  £r"5r  +  £s^s 


(Equation  (14)), 


and 


a  =  £n[e.xp(mRCR  +  m^)]  , 


mR  =  e^£  and  m„  = 


(C.6) 


'R 


S  s  • 


(C.  7) 


A  series  expansion  of  the  exponential  term  in  Equation  (C.6)  leads 
to 


A  =  Hn[l  +  (mRC^  +  m^)  +  (mRC^  +  msCg)2/2 


+  higher  order  terms] 


nh 


m„ 


=  ?.»(!  +  ~f  CR2  +  “f  q2  *  VsVs1 


(C.8) 


Estimates  of  r.m.s.  concentration  fluctuations  are  obtained  from  the 

work  of  Keeler  (1964).  Keeler's  reactor  assembly  was  similar  to  the 

one  used  in  this  study  in  terms  of  the  injector,  grids  and  the  tubular 

reactor.  He  measured  concentration  fluctuation  intensity  for  an  inert 

tracer  downstream  of  grid  and  found  that  for  a  mesh  Reynolds  number 

Re  =  800  (M  =  0.64  cm),  the  intensity  decayed  from  about  10  percent 
M 

at  5  cm  to  4  percent  at  50  cm.  The  intensities  were  based  on  the  mean 
concentration,  which  was  constant  across  the  reactor.  Based  on  these 
results , 


CR  ~  °'15 


(C.9) 


m 

• 
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y.y.v] 
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'/.V 

-  4 ’»%J 

m 


.V.V 
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Cg  ~  0.15  Cgh(t)  (C.10) 

may  be  used  as  estimates  of  fluctuations  where  g(t)  and  h(t)  are  random 

functions  of  time  which  nodulate  the  respective  amplitudes  of  C'  and 

K 

C^.  Also, 

|  g ( t)  |  <.1 

|  h  ( t )  |  <  1 

and 

g(t)  =  h(t)  =  0. 

The  largest  possible  value  of  each  term  in  Equation  (C.8)  can 

now  be  estimated  for  the  experimental  concentration  range,  and  the 

magnitude  of  A  calculated. 

For  l  =  0.001  m(gap  of  the  fiber-optic  probe), 

£  =  2438  m^/mol  (A  =  520  nm) 

K 

£g  =  888  m“/mol 

and  typical  concentrations 

C  =  C  =  0.07  mol/m\ 

K  b 

A  is  0.3  percent  of  given  by  Equation  (14). 

Similarly, 

2 

for  e  =  3216  m^/mol  (Kozicki's  value  at  A  =  520  nm) 

R 

Eg  =  1998  m“/mol 
l  =  0.001  m 

and 

C  =  C_  =  0.15  mol/m^ 

R  j 

A  is  0.8  percent  of  a\ 


The  microprocessor  (DATA  6000)  provides  the  one-sided  spectral- 

density  of  the  PUDV  signal  (G(f))  as  a  function  of  frequency  (f).  The 

2 

units  of  G(f)  are  mV  /Hz  and  that  of  frequency  are  Hz.  The  spectral 
data  needs  to  be  normalized  for  comparison  with  the  data  of  Garbini 
et  al.  (1982).  These  researchers  have  compared  their  ultrasonic  data 
with  their  own  hot-film  anemometer  data.  The  following  normalization 
(Equation  (3.11)  of  the  above  reference)  was  used  for  their  spectral 
data,  which  are  plotted  in  Figure  15: 

Efl (k)  S_(w)U 

-| - -  (D.l 

u  D  TTk/u  D 
rms  D  rms 


where  U  is  the  local  mean  velocity,  k^  is  the  magnitude  of  the  ultra- 

2 

sonic  wave  vector  (k  =  4irf  /u  ) ,  u  is  the  mean-square  value  of  the 

u  os  rms 

velocity  fluctuations,  D  is  the  pipe  diameter,  SD(to)  is  the  two-sided 

spectral-density  of  the  detected  velocity  signal,  with  units  of 
2 

(rad/s)  /(rad/s),  and  Eg(k)  is  the  one-dimensional  turbulence  spectrum. 

By  definition  (refer  to  Equations  (2.3)-(2.4)  of  the  reference), 
the  one-  and  two-sided  spectral  densities  are  related  by 


SD(W) 


G(io) 

2 


(D.  2) 


The  conversion  of  G(f)  to  G(cd)  is  obtained  as  follows.  The  zero¬ 
crossing  meter  of  the  PUDV  gives  0.25  volt  as  output  voltage  for  every 
Doppler  shift  of  1  kHz.  Also, 


o)(rad/s)=  2Tr(rad)  f(Hz) 


(D.  3) 
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Therefore, 


(rad/s) ' 
(rad/s) 


,mV3s  ,  1  kHz  103  Hz/kHz 

G(f) ("hz)  x  (o^r^n x  — r  : } 

10  mV/volt 


x  =  32  TTG(f) 


Substituting  Equation  (D.4)  into  Equation  (D.2),  the  following  is 
obtained: 


SD(u)  =  16ir  G(f ) 


(D.4) 


(D.  5) 


Equation  (D.5)  can  now  be  used  in  conjunction  with  Equation  (D.l)  to 
get  normalized  spectral  density  data.  Note  that  the  abscissa  of 
Figure  15  is  the  normalized  wave  number,  which  is  related  to  frequency 


(D.  6) 


For  the  fully  developed  pipe  flow  in  this  study: 

D  =  1.82  cm  (pipe  diameter),  f  =  10  MHz  (carrier  frequency 

of  the  PUDV) ,  and  u  =  153100  cm/s  (velocity  of  sound  in 

s 

water  at  11°C)  . 

The  spectral  data  displayed  in  Figure  15  (circles)  correspond  to  the 
numbers  listed  in  Table  D.l,  and  the  associated  V  =  1081  mV  and  V 

rms 

72  mV. 
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Case  1 :  Honeycomb  Turbulence 

Range  =  0.3487  volt 
(1.12  cm) 


Temperature 
U  (or  u  ) 


=  23-26°C 

=  29. 93x10" 3  V  (or  V 


Distance 

x' 

(cm) 


V  (mV) 


PUDV  Output 


U 

(cm/s) 


Rise  velocity  of 
bubbles  under  no 
flow  conditions 


7 

0.90 

0. 

6 

1.29 

0. 

3 

1.08 

0. 

6 

0.93 

0. 

3 

1.32 

0. 

3 

1.17 

0. 

0 

1.92 

0. 

2 

1.74 

0. 

1 

2.75 

0. 

1 

2.78 

0. 

5 

4.25 

0. 

4 

4.28 

0. 

6 

5.93 

0. 

5 

5.90 

0. 

6 

6.58 

0. 

5 

6.54 

0. 

6 

6.82 

0. 

l !4 

6.85 

0. 

>4 

6.97 

0. 

23 

6.97 

0. 

21 

7.06 

0. 

12 

7.03 

0. 

11 

7.18 

0. 

12 

7.09 

0. 

24 

7.27 

0, 

24 

7.27 

0. 

23 

7.45 

0. 

23 

7.39 

0. 

23 

7.42 

0, 

24 

7.42 

0 

L7 

0.24 

0 

L7 

0.24 

0 

Downstream  distance  from  the  honeycomb. 


Not  corrected  for  noise;  corrected  value  = 

V  -  (V  )  .  ,  .  • 

rms  rms  rise  velocity 


rise  velocity 
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Case  2:  (Honeycomb  +  Injector)  Turbulence 


Range  =  0.3487  volt 
(1.12  cm) 


Temperature 

U  (or  u  ) 
rms 


25-26°C 

29.93xl0“3  V  (or  V  ) 


Rise  velocity  of 
bubbles  under  no 
flow  conditions 


9 

10 


17 

17 


0.27 

0.30 


0.51 

0.51 


+Axial  distance  downstream  of  the  injector  needle  tips. 
'‘Not  corrected  for  noise. 
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Case  4b:  (Honeycomb  +  Injector  +  Screen  (M^))  Turbulence 
(Screen  upstream  of  the  injector) 


Range  =  0.4750  volts 
(1.52  cm) 


Temperature  =  11°C 


x  =  39.0  cm 

Radial 
Distance 
Y  (cm) 

0.75 


U  (or  u  )  =  30.56  V  (or  V  )  cm/s 
rms  rms 


Distance"1" 

X 

(cm) 

V  (mV) 

V  (mV) 

rms 

j. 

u" 

(cm/s) 

u  * 

rms 

(cm/ s) 

13.6 

235 

25 

7.18 

0.76 

233 

23 

7.12 

0.70 

233 

25 

7.12 

0.76 

18.6 

237 

25 

7.24 

0.76 

242 

26 

7.40 

0.  79 

237 

30 

7.24 

0.92 

23.7 

245 

27 

7.49 

0.83 

245 

26 

7.49 

0.79 

244 

25 

7.46 

0.76 

28.8 

256 

25 

7.82 

0.76 

252 

27 

7.70 

0.83 

254 

28 

7.76 

0.86 

33.9 

257 

26 

7.85 

0.79 

256 

27 

7.82 

0.83 

262 

25 

8.01 

0.76 

39.0 

262 

26 

8.01 

0.79 

263 

26 

8.04 

0.79 

266 

26 

8.13 

0.79 

44.0 

264 

31 

8.07 

0.95 

268 

27 

8.19 

0.83 

267 

26 

8.16 

0.  79 

266 

27 

8.12 

0.83 

259 

28 

7.92 

0.86 

261 

25 

7.98 

0.  76 

258 

28 

7.88 

0.86 

263 

28 

8.04 

0.86 

256 

28 

7.82 

0.86 

Axial  distance  downstream  of  the  injector  needle  tips. 

J- 

Not  corrected  for  noise. 
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l-Naphthol  is  a  primary  standard  and  can  be  directly  used  to 

prepare  a  standard  solution  of  it.  As  shown  in  Figure  F.l,  it  has  a 

characteristic  peak  in  its  absorption  spectrum  at  A  =  245  nm  (Na^CO^/ 

3 

NaHC03  buffer,  pH  =  10.00  ±  0.002,  T  =  298  K,  I  =  40  g  ion/m  ).  A 
typical  calibration  data-set  is  listed  below. 


Calibration  of  l-Naphthol  Solution 
(C. .  =  0.1998  mol/m3) 


(mol/m  ) 


Absorbance 
at  245  nm 
A 


* \“w  ' 

A.-..-..-' 


»  "  »  "V  * 

r  v  ■>  V 
- 

fvV.Vx; 


"  %  lN 
_  *  m  %  “  * 

>:# 

$ 

.V« 

»  A  ■.  > 


A  linear  regression  results  in  the  following  equation: 


A  =  23.0918  CA  +  0.0009  (F.l 

A 

This  equation  is  then  used  to  determine  the  concentration  of  a  sample 
taken  from  the  feed  tank  A  for  a  mixing-reaction  run.  Typically,  the 
absorbance  of  such  samples,  at  16:1  dilution  with  the  buffer,  was 
found  in  the  range  0.40-0.42,  which  corresponded  to  C^q  =  0.294  - 
0.309  mol/m3. 
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Experiment 

Colorimeter  Output 

Date: 

No.  9 

VX  (mV) 

March  9,  1985 

Reference 

3.7 

8.8 
14.0 
19.0 
24.0 
29.0 

Reference 


Exit  1 
Reference 
Exit  2 
Reference 
Exit  3 
Reference 


(mV/mo 1/m  ) 


(mV/mol/m  ) 


Injector  clogged,  hence  purged  with  high  flow  of  B 
+0btained  by  regression  with  MINITAB. 


(mol/m  ) 


(mV/mol/m 


kA 

ks 

(mV/mol/m^) 


1369  1224 


654  836 


916 


819 


800 


857 


903 


168 


Experiment 

No.  12 


Colorimeter  Output 
7  (mV) 


Date : 

April  2,  1985 


s 


«  _  a  _  *•  b  a 

VvJ\' 

(-.%v 

•-W 

•Jv3 

A  *■ 


X 

A  (nm) 

_  + 

CR 

*s+ 

(cm) 

490 

520 

545 

450 

470 

(mol/m 

3) 

Reference 

781 

425 

230 

992 

981 

14.6 

913 

57  6* 

382* 

1085* 

1083 

0.089 

0.038 

19.7 

926 

598 

386 

1101 

1099* 

0.090 

0.053 

24.8 

943 

613 

403“ 

1108 

1105 

0.102 

0.054 

29.9 

946 

621 

411 

1116* 

1119 

0.104 

0.060 

34.9 

949 

620 

422* 

1120 

1133 

0.103 

0.065 

40.0 

951 

631* 

438 

1122* 

0.108 

0.068 

Reference 

787 

430 

241 

994 

Exit  1 

1154 

742 

524 

1176 

1041 

Reference 

883 

435 

240 

992 

823 

0.180 

0.075 

Exit  2 

1151 

738 

521 

1170 

1039 

Reference 

879 

436 

240 

995 

821 

kA 

K 

(mV/mol/m  ) 

1119 

1369 

1224 

654 

836 

kA 

ks 

3 

(mV/mol/m  ) 

916 

819 

800 

857 

903 

Injector  clogged,  hence  purged  with  high  flow 
^Obtained  by  regression  with  MINITAB. 
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Experiment 

Colorimeter  Output 

Date: 

i. 

a 

No .  16 

^  (mV) 

April  16,  1985 

S' 

o’1 

A 

Ref  erence 

14.6 

19.7 

24.8 

29.9 

34.9 
40.0 

Reference 


Exit  1 
Reference 
Exit  2 
Reference 
Exit  3 
Reference 


(mV/mol/m  ) 


(mV/mol/m  ) 


(mol/m  ) 


Ll 

I 


m 


Obtained  by  regression  with  MINITAB 


1322  1166 


612  790 


(mV/mol/m  ) 


605 


547 


543 


546 


602 


Experiment 
No.  22 


Colorimeter  Output 
(mV) 


Date : 

May  4,  1985 


X 

X  (nm) 

—  + 

CR 

C  + 

(cm) 

490 

520 

545 

450 

470 

(mol/m 

3) 

Reference 

919 

560 

306 

1042 

870 

5.1 

1006 

663 

399 

1087 

935 

0.087 

0.000 

10.2 

1038 

688 

437 

1117 

973 

0.0773 

0.041 

15.6 

1079 

725 

464 

1144 

994 

0.093 

0.061 

20.7 

1093 

746 

484 

1163 

1007 

0.099 

0.077 

25.8 

1110 

770 

503 

1167 

1022 

0.123 

0.064 

30.9 

1128 

770 

518 

1183 

1035 

0.102 

0.113 

Reference 

929 

559 

311 

1044 

884 

Exit  1 

1357 

1001 

715 

1305 

1168 

Reference 

966 

569 

317 

1054 

866 

Exit  2 

1350 

1002 

717 

1294 

1170 

Reference 

965 

567 

316 

1049 

869 

0.257 

0.125 

Exit  3 

1347 

997 

712 

1289 

1167 

Reference 

965 

567 

317 

1040 

864 

4 

3 

(mV/mol/m  ) 

1112 

1360 

1193 

650 

830 

kx 

ks 

(mV/mol/m  ) 

744 

716 

704 

657 

723 

+Obtained  by  regression  with  MINITAB. 
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Exper imenc 

Colorimeter  Output 

Date 

: 

No .  25 

vA  (mV) 

May  9 , 

1985 

X 

A  (nm) 

—  + 

CR 

(cm) 

490 

520 

545 

450 

470 

(mol/m 

,3> 

Reference 

897 

464 

252 

1019 

848 

4.1 

1083 

678 

453 

1142 

997 

0.129 

0.06 

9.2 

1122 

697 

469 

1155 

1007 

0.140 

0.05 

14.3 

1126 

720 

495 

1184 

1036 

0.129 

0.11 

16.8 

1142 

721 

507 

1176 

1041 

0.127 

0.11 

19.3 

1151 

730 

506 

1184 

1044 

0.131 

0.11 

24.5 

1156 

741 

511 

1189 

1060 

0.139 

0.12 

Reference 

894 

474 

255 

1024 

858 

Exit  1 

1252 

869 

625 

1241 

1131 

Reference 

875 

444 

232 

994 

823 

Exit  2 

1251 

867 

627 

1242 

1127 

0.239 

0.14 

Reference 

877 

442 

232 

991 

825 

Exit  3 

1255 

867 

627 

1243 

1130 

Reference 

876 

444 

231 

995 

825 

.A 

\ 

(mV/mol/m^) 

1112 

1360 

1193 

650 

830 

3 

(mV/mol/m  ) 

744 

716 

704 

657 

723 

+Obtained  by 

regression 

with 

MINITAB. 
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Experiment 
No.  26 


Reference 


Colorimeter  Output 
VA  (mV) 


\  (nm) 
545 


1360  1193 


Date: 

May  13,  1985 


(mol/m'"’) 


0.127 

0.124 

0.128 

0.151 

0.146 

0.143 


0.053 

0.083 

0.108 

0.084 

0.094 

0.105 


0.218  0.143 


Reference 


Exit  1 
Reference 
Exit  2 
Reference 


(mV/mol/m^) 


(mV/mol/m  ) 


Obtained  by  regression  with  MINITAB. 

One  bubble  on  the  screen. 

Two  bubbles  on  the  screen;  effectively  10  percent  screen  blockage 


v.v.i 

CvS- 


8$ 

m 


m 

m 

mi 


VAyivtv5!rovi'iV.,i' 


Experiment 
No.  27 


Date: 

May  14,  1985 


Colorimeter  Output 
VX  (mV) 


X 

X  (nm) 

_  + 

CR 

<S+ 

(cm) 

490 

520 

545 

450 

470 

(mol/ 

m3) 

Reference 

875 

389 

199 

975 

815 

4.1 

1046* 

591 

397 

1090 

955 

0.115 

0.057 

9.2 

1075 

619 

416 

1049 

968 

0.152 

0.036 

14.3 

639 

435 

1120 

992 

0.144 

0.074 

16.8 

649 

453 

1131 

1007 

0.132 

0.106 

19.3 

1099 

649 

454 

1133 

1012 

0.127 

0.115 

24.5 

1107 

660 

454 

1137 

1011 

0.148 

0.096 

Reference 

877 

399 

083 

819 

Exit  1 

1224 

780 

560 

1212 

1092 

Reference 

884 

399 

203 

979 

817 

Exit  2 

1230 

787 

565 

1212 

1096 

Reference 

885 

397 

203 

977 

815 

0.212 

0.147 

Exit  3 

1235 

793 

579 

1211 

1104 

Reference 

884 

399 

202 

977 

817 

(mV/mol/m3) 

1112 

1360 

1193 

650 

830 

ks 

744 

716 

704 

657 

723 

(mV/mol/m3) 


Experiment 
No.  30 


Colorimeter  Output 
V^(mV) 


Date: 

May  21,  1985 


X 

X(nm) 

_  + 

CR 

s+ 

(cm) 

490 

520 

545 

450 

470 

(mol/m^) 

Reference 

547 

375 

178 

934 

770 

4.1 

761 

590 

357 

1050 

925 

0.117 

0.072 

9.2 

784 

602 

388 

1061 

931 

0.141 

0.054 

14.3 

770 

610 

395 

1082 

944 

0.059 

0.189 

16.8 

788 

604 

390 

1083 

949 

0.098 

0.140 

19.3 

613 

403 

1079 

950 

0.098 

0.140 

24.5 

613 

399 

1075 

942 

0.112 

0.112 

Reference 

394 

184 

937 

111 

Exit  1 

1117 

828 

596 

1192 

1085 

Reference 

725 

380 

171 

936 

773 

0.216 

0.215 

Exit  2 

1140 

851 

617 

1214 

1111 

Reference 

731 

384 

177 

936 

767 

,A 

“r 

(mV/mol/m^) 

1148 

1384 

1156 

672 

853 

kA 

ks 

[mV/mol/m^) 

732 

733 

688 

570 

655 

Obtained  by  regression  with  MINITAB. 


209 

444 


Reference 


0.152 

0.173 

0.188 

0.174 

0.189 

0.163 


0.053 

0.096 

0.104 

0.132 

0.119 

0.170 


Reference 


Exit  1 
Reference 
Exit  2 
Reference 
Exit  3 
Reference 


(mV/mol/m  ) 


(mV/mol/m  ) 


Obtained  by  regression  with  MINITAB. 
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ExperimenC 
No.  33 


Colorimeter  Output 
VA  (mV) 


Date: 

May  26,  1985 


•  W  ! 
S.VJ ■  1 N 

\  v" 
»  Vv 

.S£vv 


X 

X  (mi) 

r  +  c  + 

lr 

(cm) 

490 

520 

545 

450 

470 

(mol/m^) 

183 


Experiment  Colorimeter  Output  Date: 


No.  34 

VA  (.mV) 

June  14, 

1985 

—  *+■ 

r  + 

X 

X  (nm) 

CR 

cs 

(cm) 

490 

520 

545 

450 

470 

3, 

(mol/m 

i  ) 

Reference 

848 

531 

244 

955 

790 

4.0 

911 

607 

318 

997 

839 

0.057 

0.000 

9.1 

957 

650 

345 

1024 

875 

0.063 

0.047 

14.1 

982 

687 

383 

1046 

905 

0.082 

0.064 

16.  7 

697 

3  92 

1054 

910 

0.087 

0.072 

19.2 

699 

406 

1061 

918 

0.084 

0.  090 

24.3 

720 

422 

1076 

937 

0.086 

0.116 

Reference 

855 

540 

243 

966 

797 

Exit  1 

1329 

972 

749 

1284 

1182 

Reference 

859 

450 

257 

970 

802 

Exit  2 

1332 

973 

752 

1288 

1184 

Reference 

858 

453 

264 

967 

802 

0.258 

0.276 

Fxi»-  3 

1332 

977 

755 

1285 

1186 

Reference 

858 

449 

263 

971 

799 

>4 

1117 

1342 

1186 

645 

818 

(mY/mol/m^) 

-A 

ks 

670 

641 

608 

543 

622 

(mV/mol/m^) 

Obtained  by 

regression 

with 

MINITAB. 

Experiment 

Colorimeter  Output 

Date : 

No.  35 

(mV) 

June  17,  1985 

Experimenc 
No.  36 


Date : 

June  18,  1985 


n  i 


Experiment 
No.  3  7 

Colorimeter  Output 
MX  (mV) 

Date : 

June  19,  1985 

X 

X(nm) 

(cm) 

490 

520 

545 

450 

470 

(mol/m^) 

Reference 

787 

534 

159 

913 

742 

16.1 

859 

615 

224 

956 

797 

0.061 

21.2 

862 

608 

240 

962 

803 

0.070 

26.3 

869 

62/ 

249 

967 

810 

0.077 

32.4 

872 

632 

259 

973 

817 

0.082 

37.5 

878 

645 

256 

973 

815 

0.084 

41.5 

883 

635 

262 

981 

819 

0.088 

Reference 

788 

537 

158 

915 

744 

Exit  1 

1102 

534 

303 

1002 

858 

Reference 

971 

382 

162 

919 

754 

Exit  2 

1100 

535 

305 

1003 

857 

0.119 

Reference 

969 

381 

162 

919 

754 

Exit  3 

1100 

534 

304 

1004 

859 

Reference 

971 

381 

162 

919 

754 

kR 

1117 

1342 

1186 

645 

818 

(mV/mol/m^ 

) 

^Obtained 

through 

approximate 

analvsis. 

k>"; 

iSO 

■'>>>> 


E 


